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PREFACE 

The purpose of this series of technical reports is to provide current documentation of the SIMBIOS Project activities, 
NASA Research Announcement (NRA) research status, satellite data processing, data product validation and field 
calibration. This documentation is necessary to ensure that critical information is related to the scientific community 
and NASA management. This critical information includes the technical difficulties and challenges of combining ocean 
color data from an array of independent satellite systems to form consistent and accurate global bio-optical time series 
products. This technical report is not meant to substitute for scientific literature. Instead, it will provide a ready and 
responsive vehicle for the multitude of technical reports issued by an operational project. 
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Chapter 1 

SIMBIOS Background 

Charles R. McClain 
NASA Goddard Space Flight Center 
Greenbelt , Maryland 


The Sensor Intercomparison and Merger for 
Biological and Interdisciplinary Ocean Studies 
(SIMBIOS) program was conceived in 1994 as a 
result of a National Aeronautics and Space 
Administration (NASA) management review of the 
agency's strategy for monitoring the bio-optical 
properties of the global ocean through ocean color 
remote sensing from space. At that time, the NASA 
ocean color flight manifest included two data buy 
missions, the Sea-viewing Wide Field-of-view Sensor 
(SeaWiFS) and Earth Observing System (EOS) Color, 
and two sensors, Moderate Resolution Imaging 
Spectroradiometer (MODIS) and Multi-angle Imaging 
Spectro-Radiometer (MISR), scheduled for flight on 
the EOS-AM (1998) and -PM (2001) satellites. 

Considerable effort was spent by Dr. McClain 
(Project Scientist) and Dr. Kirk (Study Manager) on 
examining mission scenarios for EOS Color because 
of the slips being encountered with SeaWiFS. 
However, with the delay of SeaWiFS and an uncertain 
launch schedule, it was not clear that EOS Color was 
needed to fill a potential gap between SeaWiFS and 
MODIS, especially when five additional ocean color 
systems with similar global capabilities [Ocean Color 
and Temperature Sensor (OCTS), Japan; Global 
Imager (GLI), Japan; Polarization Detecting 
Environmental Radiometer- 1 (POLDER- 1) and -2, 
France; and Medium Resolution Imaging 
Spectrometer (MERIS), European Space Agency] and 
several other non-global missions by Argentina, 
Germany, Taiwan, India, Korea, the U S. Navy, and 
the People’s Republic of China, were planned for 
launch during the late 1990s. 

The review led to a decision that the international 
assemblage of ocean color satellite systems provided 
ample redundancy to assure continuous global 
coverage, with no need for the EOS Color mission. At 
the same time, it was noted that non-trivial technical 
difficulties attended the challenge (and opportunity) of 
combining ocean color data from this array of 
independent satellite systems to form consistent and 
accurate global bio-optical time series products 
(Figure 1). Thus, it was announced at the October, 
1994 EOS Interdisciplinary Working Group meeting 
that some of the resources budgeted for EOS Color 
should be redirected into an intercalibration and 
validation program. NASA Goddard Space Flight 


Center (GSFC) was directed to develop a plan for 
submission to NASA Headquarters (HQ) by May 

1995. As a result of the directive from NASA/HQ, the 
ocean color group lead by Dr. McClain at 
NASA/GSFC organized an international 
organizational meeting in February 1995 at the 
University of Miami Rosenstiel School for Marine and 
Atmospheric Sciences. The objective was develop a 
conceptual plan for a comparison program. The plan 
(Ocean Color Multisensor Data Evaluation and 
Utilization Plan, under SIMBIOS Documents, 
http://simbios.gsfc.nasa.gov/) which outlined NASA’s 
contribution to the international effort was completed, 
submitted, externally reviewed, and revised in 1995. 

Based on the final approved plan, a NASA 
Research Announcement (NRA) was released in July 

1996, and the SIMBIOS Project Office was 
established at the NASA (GSFC) in January 1997 (co- 
located with the SeaWiFS Project). The initial 
SIMBIOS Program was scoped for five years (1997- 
2001) and included separate support for a science 
team (NRA selections) and the Project Office. 

Dr. Mueller (San Diego State University) acted 
as an interim project manager at NASA/GSFC under a 
one-year assignment to assist in getting the project 
office organized and the science team contracts 
executed. His assistance in this capacity was essential 
as the SIMBIOS Project was beginning just as the 
SeaWiFS Project was preparing for launch. 

In parallel with the NASA SIMBIOS Program 
planning, the international effort was being organized. 
The initial meeting was held in Vittoria, British 
Columbia during September, 1995. As a result of the 
recommendations from that meeting, the International 
Ocean Colour-Coordinating Group (IOCCG) was 
formed (http://www.ioccg.org/). The IOCCG 
presently operates under the auspices of the Scientific 
Committee on Oceanic Research (SCOR) and 
chairmanship of Dr. Platt (Bedford Institute of 
Oceanography). The IOCCG meets one or two times 
per year and is generating a series of special reports on 
topics essential to the coordination of the international 
ocean color community (e g., IOCCG, 1998). 

During calendar year 1997, Dr. Mueller’s tenure 
as SIMBIOS Project Manager, the primary 
accomplishments included the following: 
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• The science team contracts or interagency 
agreements were negotiated (3 -year contracts) 
and first year funds were distributed. Several 
additional investigations, primarily for 
atmospheric correction studies, were solicited and 
funded. 

• The initial Project Office staffing actions were 
completed. 

• The SIMBIOS processing systems were 

purchased and installed. 

• An instrument pool was established. Submersible 
instruments for field studies were purchased 
through some of the science team contracts. The 
Project also purchased 12 MicroTops hand-held 
sun photometers, 2 PREDE (Japanese built) sun 
photometers, and a micro-pulse lidar as part of the 
instrument pool. Twelve CIMEL (French built) 
sun photometers were also purchased to augment 
the Aerosol Robotic Network (AERONET) for 
deployment at coastal sites. 

• A satellite overflight prediction service was 
established to assist in coordinating field studies 
in synchronizing data collection with satellite 
passes (http://simbios.gsfc.nasa.gov/, and http: // 
simbios.gsfc.nasa.gov/~schedule/Predictions/ 
Current_Cruises.html). 

• The first SIMBIOS science team meeting was 
held in August at Solomons Island, Maryland. 

During the second year of the SIMBIOS Project, 
Dr. McClain assumed project management of both the 
SeaWiFS and SIMBIOS as both Dr. Cleave and Dr. 
Mueller stepped down in their roles as project 
managers of these two projects, respectively. In the 
second year, many of the original objectives of 
SIMBIOS began to come to fruition, some of which 
are discussed in more detail in subsequent chapters. 
Project accomplishments included: 

• The first SIMBIOS intercalibration round-robin 
experiment (Riley and Bailey, 1998) was 
completed. 

• The CIMEL sun photometers were re-engineered 
by Meridian Engineering of Hanover, Maryland 
to improve their survivability in the field. 
Confirmed (delivered or negotiated) CIMEL sites 
include Lanai (Hawaii) with a backup in 
Honolulu, Ascension Island, Bahraine, Turkey 
(Black Sea), Tahiti, Wallops Island (Virginia), 
and South Korea. Additional sites being 
considered are Perth, Australia and Iceland. 
Delivery of equipment for calibrating polarized 
bands was received from the University of Lille 
as part of the collaboration with the AERONET 
group at NASA/GSFC [Brent Holben, principal 
investigator (P I.)]. An automated procedure for 
extracting data from the AERONET archive for 
satellite match-up comparisons was completed. 


• The micropulse lidar environmental container was 
re-engineered by Meridian Engineering to 
improve temperature stability. 

• A second copy of the SeaWiFS Transfer 
Radiometer (SRX), the SXR-2, was built by 
Reyer Corporation of New Market, Maryland for 
use in the calibration round-robin. 

• A Marine Optical Spectroradiometer (MOS) 
ground station subsystem was purchased from 
Antrix Corporation Limited (Bangalore, India) for 
installation at Wallops Flight Facility. 

• Data processing (Level -0 to Level-3) software for 
OCTS and MOS were completed. This activity 
included evaluations of the navigation, band 
coregistration, calibration, destripping, 
atmospheric correction, and bio-optical 
algorithms. Match-up comparisons with field data 
and OCTS were also completed using match-up 
subscenes provided by National Space 
Development Agency of Japan (NASD A) Earth 
Observation Research Center (EORC). 

• SeaWiFS Bio-Optical Archive and Storage 
System (SeaBASS) holdings were greatly 
expanded as the science team field activities got 
underway. 

• The SIMBIOS website was substantially 
enhanced including monthly reports posted by the 
science team investigators. 

• Science team contract performance was reviewed 
and the contracts were revised as necessary and 
renewed. 

• The second science team meeting was held at 
Scripps Institution of Oceanography (SIO) in 
September 1998. 

In 1999, the SIMBIOS Project will finalize the 
processing of it holdings of OCTS data, i.e., all the 
OCTS data collected at Wallops Flight Facility (east 
coast of North America, Gulf of Mexico, and the 
Caribbean Sea), which should be completed by early 
spring. The OCTS processing software will be 
incorporated into SeaWiFS Data Analysis System 
(SeaDAS) Version 3.3. The next ocean color 
instrument data set to be evaluated will be ADEOS- 
1/POLDER in cooperation with Centre National 
d’Edudes Spatialle (CNES) and the POLDER science 
team. The SIMBIOS Project and science team will 
assist the ROCSAT Ocean Color Imager (OCI) 
(launched in January) and MODIS (launch scheduled 
for July) instrument teams in their evaluations, 
primarily by providing match-up data for calibration 
and validation analyses. The SIMBIOS Project will 
also assist NASAHQ in drafting the second 
SIMBIOS NRA which is scheduled for release 
sometime this summer. Finally, the third science team 
meeting will be held in September in the Washington, 
D C. area. 
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Figure 1 . Ocean Color Missions 
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Chapter 2 

SIMBIOS Program and Project Goals 

Giulietta S, Fargion 
SAIC General Sciences Corporation 
Beltsville, Maryland 


2.1 PROGRAM GOALS 

The purpose of the SIMBIOS project is to 
develop a methodology and operational capability to 
combine data products from various ocean color 
missions in a manner that ensures the best possible 
global coverage and best exploits the complementary 
missions of the sensors. NRA 96-MTPE-04 
specifically defines the scientific goals in support of 
this undertaking. “Specifically, the objectives are (1) 
to quantify the relative accuracies of the products from 
international ocean color missions, (2) to improve the 
level of confidence and compatibility among the 
products, and (3) to generate merged, improved level- 
3 products” (NASA, 1996). A major goal of the 
Program is to foster collaborations with the various 
space agencies and science working groups to assist 
each other in achieving these objectives. Presently, 
SIMBIOS is collaborating with Japan on OCTS and 
GLI, with India and Germany on MOS and just 
started a relationship with France on POLDER. 

2.2 PROJECT OBJECTIVES 

The SIMBIOS Program consists of the SIMBIOS 
Science Team and the SIMBIOS Project. The 
SIMBIOS Science Team, as initially defined by the 
NRA selections, was expanded to include additional 
investigations for atmospheric correction algorithm 
validation. The SIMBIOS Project has been established 
to provide support and coordination for the SIMBIOS 
Program such as administration, project 
documentation, complete data processing of 
complementary ocean color missions, support for each 
of the individual elements, and oversight of data 
processing systems supported by the Project. 
“Specifically, the SIMBIOS Project is responsible for 
the following: 

• Assistance in the collection, processing, archiving 
and documentation of the calibration (including 
instrument round-robins, pre- and post-launch 
sensor calibration) and bio-optical data sets; 

• Maintenance of a data archive and storage system 
for tracking instrument calibrations and validating 
bio-optical algorithms or merger methodologies 
with in situ data; 


• Planning and execution of instrument and 
analysis round-robins which will include 
coincident intercomparisons of hardware and 
software performance; 

• Oversight of the regular revision of the ocean 
optics protocols which will establish the 
methodologies and standards for calibrating 
instruments, collecting data, and producing final 
results from approved analysis procedures. 

• Assistance in data acquisition during field studies; 

• Management of a centralized facility for 
providing portable field sources which will be 
used for iracking the" calibration of instruments 
while they are used in the field; 

• Assistance in product validation and quality 
control assessment; 

• Technical monitoring of the prototyping, 
deployment, or production of new technology' 
developments; 

• Assistance in the development, evaluation, and 

implementation of product merger schemes and 
algorithms; and § 

• Administer and oversee the regular meeting of 
working subgroups concerned with specific issues 
important to the success of the SIMBIOS Project” 
(NASA, 1996). 

The SIMBIOS Program is augmented by the 
participation of the NASA-supported MODIS Oceans 
Team and SeaWiFS Calibration and Validation 
Program. SIMBIOS calibration activities are 
accomplished using the calibration activities of the 
various independent participating missions, 
coordinating activities of the SIMBIOS Project, and 
complementary activities of individual investigators 
selected by the NRAs. The SIMBIOS Project goal is 
to help science users with a long time series of 
calibrated radiances extending across the boundaries 
of individual missions. Close cooperation with the 
SeaWiFS Project is essential to the success of the 
SIMBIOSProject. The SIMBIOS Project is co-located 
with the SeaWiFS Project at GSFC, and many of 
SIMBIO S sta ff are shared with SeaWiFS to guarantee 
a close cooperation to quantify the relative accuracy of 
the data products from different ocean color sensors 
and the improve level of confidence and compatibility 
among the products. 
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Table 1, SIMBIOS Project 


SIMBIOS Project Address 

NASA/Goddard Space Flight Center 

SIMBIOS Project 

Code 970.2, Building 28 

Greenbelt, MD 20771 

Telephone: (301)286-9676 

FAX: (301) 286-0268 

E-mail: project@simbios.gsfc.nasa.gov 

Personnel of SIMBIOS Project 

Sean Bailey 

FutureTech Corporation 

Robert Barnes * 

SAIC General Sciences Corporation 

Patty Clow * 

NASA 

Giulietta Fargion 

SAIC General Sciences Corporation 

Gene Feldman * 

NASA 

Ian Durham * 

SAIC General Sciences Corporation 

Bryan Franz * 

SAIC General Sciences Corporation 

Lynne Hoppel * 

NASA 

Alice Isaacman 

SAIC General Sciences Corporation 

Charles McClain * 

NASA 

Fred Patt * 

SAIC General Sciences Corporation 

Christophe Pietras 

SAIC General Sciences Corporation 

Tom Riley 

NASA 

India Robinson * 

NASA 

Kellie Ruebens * 

SAIC General Sciences Corporation 

Brian Schieber * 

SAIC General Sciences Corporation 

Lin wood Smith * 

FutureTech Corporation 

Paul Smith * 

SAIC General Sciences Corporation 

Sung Lee * 

SAIC General Sciences Corporation 

Tamara Tucker * 

SAIC General Sciences Corporation 

Grey Valenti * 

SAIC General Sciences Corporation 

Menghua Wang 

Univ. of Maryland Baltimore 

Bill Woodford * 

FutureTech Corporation 


Note: total civil servant staffs 2.5 man year 


* shared with SeaWiFS Project 

2.3 PROJECT ORGANIZATION 

The SIMBIOS Project organization structure is 
scoped to support the primary activities: 1) sensor 
engineering and calibration, 2) satellite data 
processing; 3) data product validation, 4) SIMBIOS 
Science Team, and 5) SIMBIOS Project Office. The 
SIMBIOS Project organization chart is shown in 
Figure 1 and illustrates the individual subtasks under 
each primary activity with staff associations. The 
SIMBIOS Project key personnel and their areas of 
responsibility at GSFC are also given in Table 1. The 
satellite data processing elements requires close 
coordination between the data product validation and 
the sensor engineering and calibration. A complete 
overview of the primary activities 1 to 4 will be 
described in the following chapter. The Project Office 
administrative structure is as follows: system support, 
administrative support, technical support, science 
tpam support, procurement support, and resources 


support. At the present time, only SeaWiFS and 
SIMBIOS participants and collaborators associated 
with other relevant programs, and science teams who 
actively contribute data or algorithms, are given 
access to SIMBIOS Project resources, such as the 
SeaWiFS Bio-optical Archive and Storage System 
(SeaBASS) and Field Instrument Pool. This limitation 
is to protect the interests of the researchers currently 
submitting to the SeaBASS archive. 

The SIMBIOS Project policy is to have all 
documents easily available to die SIMBIOS Pis and to 
the larger ocean color scientific community. Access to 
the information is gained through the SIMBIOS 
website. The SIMBIOS website was organized to 
serve as the main information resource to access the 
Project activities. Project Office and Science Team 
The website is organized under five main topics: 
News and Information, Support Services and 
Schedule, Project Status, Instrument Pool and 
Contacts. Briefly, under Project Status are found 
monthly progress updates from the Project and the 
Science Team. The Project update covers processing 
and calibration accomplishments for the various 
missions, the instrument pool and real-time cruise 
support, and other key events or Project issues. The 
Science Team oversight is divided into two categories: 
funded Pis and international collaboration Under 
Status of Contracted Pis, the abstracts of the funded 
research, monthly progress reports, documents in 
process, and other key events can be found. 

Under News and Information, is an archive of 
most of the produced SIMBIOS documentation and 
calendar of events. In this section, a user can easily 
find out which conferences or symposia the Project 
staff is attending, documents on SIMBIOS workshops 
and team meetings, and special activities, such as the 
SIMBIOS participation in Indian Ocean Experiment 
(INDOEX). An online bibliography includes 
SIMBIOS technical memoranda or papers presented 
and published, with online abstracts. SeaWiFS 
scheduling and the SeaBASS database are accessible 
online. For further information see http://simbios.gsfc. 
nasa.gov. All sections are updated as needed. 

In addition, the SIMBIOS Project has an internal 
web interface to closely follow the contractual status 
of the funded Pis. Contractual issues such as data 
delivery, contract evaluation performance, contract 
modifications, requested extensions, procurement 
support performance, and other items are tracked here. 
In addition, the SIMBIOS Project holds bimonthly 
staff meetings, weekly review on technical papers or 
project issues (such as calibration, in situ matchups, 
etc ), and key staff attends the SeaWiFS meetings 
Monthly SIMBIOS Project progress reports are 
generated and posted to the SIMBIOS website. The 
Project also gives a formal oral presentation to the 
Earth Science Directorate (Code 900) management. 
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Chapter 3 

SIMBIOS: NRA Contracts 

Giulietta S. Fargion 
SAIC General Sciences Corporation 
Beltsville, Maryland 


3.1 PROJECT RESPONSIBILITY 

Contract administration involves ensuring 
compliance with contract terms and conditions during 
the performance of the contract. The performance of 
SIMBIOS Project contracts are monitored in several 
ways by several individuals in the Project office. 
These individuals include a Contracting Officer 
specialist, a Resource/Finance person, the Assistant 
Project Manager with technical background and 
experience with projects and their costs, and the 
Program Manager. The SIMBIOS Project contract 
personnel and their areas of responsibility at GSFC are 
given in Table 1. The most important step of contract 
administration is to review the requirements and 
specific obligations set forth in the contract. This role 
is performed primarily by the Assistant Project 
Manager. In matters related to contract laws and 
interpretation, assistance is given by the Contracting 
Officer. In matters related to the status of payments or 
provisions relating to payments, assistance is given by 
the Resource/Finance person. While SIMBIOS 
contractors have historically addressed questions 
about contract administration to the Assistant Project 
Manager, they are entirely free to address questions 
directly to the appropriate personnel. 

A wide selection of contract types is available for 
use between the Project office and contractors. This 
selection provides a wide selection of services and a 
large degree of flexibility of contract conditions (FAR, 
1998). GSFC Procurement Office selected a Research 
and Development Service, firm-fixed price (FFP) 
contract-type for the NASA Research Announcement 
(NRA) selections. The FFP contract provides for a 
price that is not subject to any adjustment by reason of 
the cost experience of the contractor in the 
performance of the contract. Use of the FFP contract 
imposes a minimum administrative burden on the 
contracting parties. Contracts are very different from 
the usual scientific research grants in three main 
aspects: 1) a work statement concerning the end 
objectives of the research (i.e., the proposal), 2) the 
technical data to be delivered under the contract with 
time delivery specifications, and 3) scientific and 
technical reports, consistent with the objectives of the 
effort involved, as a permanent record of the work 
accomplished under the contract. 


3.2 CONTRACTS OVERVIEW 

The initial SIMBIOS program was scoped for five 
years (1997-2001) and includes support for the NRA 
selections (i.e., science team) and the SIMBIOS 
Project. The Interim SIMBIOS Office (ISPO) 
successfully negotiated specific work statements for 
the 21 NRA selections and issued contracts starting 
from late summer 1997. SIMBIOS funded contracts, 
PI names and delivery requirements are summarized 
in Table 2. 

NASA Procurement requires formal evaluations 
for all contracts at the end of each contract year These 
evaluations are to go into a database and will be 
shared with the Pis’ institutions. The performance 
evaluation has four categories: quality, time, other and 
price (price is not relevant here because contracts are 
fixed cost). SIMBIOS Project Office procedure is to 
coordinate an inside panel to perform an across-the- 
board evaluation of all funded contracts. Under quality 
will be considered 1) data quality and completeness, 
2) ancillary information provided on the data 
(metadata), 3) the data's usefulness in relation to 
SIMBIOS goals, i.e., calibration, validation, and 
algorithm development, and 4) quality of technical 
reports. Time is a mixed bag, but will be viewed in 
respect to data and documentation (monthly, year-end 
reports, and special topic publications) and delivery 
times. Under “other” will be considered 1) scientific 
publications and scientific achievements, 2) science 
team collaboration and involvement, and 3) other 
significant events occurring during the contract period 
evaluate. 


Table 1. Contract key personnel 


SIMBIOS NRA Organization 

Contracting Officer: 

Lynne Hoppel 

Resource/Financial Officer: 

Patty Clow 

Assistant Project Manager: 

Giulietta Fargion 

Program Manager: 

Charles McClain 
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3.3 SCIENCE TEAM 

NASA SIMBIOS Science Team Principal 
Investigators are composed of those selected under the 
NRAs, some members of the MODIS Ocean Teams, 
and certain members the SeaWiFS Project. The 
Science Team can be grouped under three working 
areas: 1) Ocean Bio-optical and Sensor 

Characterization Studies, 2) Data Merger Studies, and 
3) Atmospheric Correction Studies. There are many 
more U S. and international co-investigators and 


collaborators actively participating in the NASA 
component of the international SIMBIOS program. 

The funded Principal Investigators and 
collaborators are shown in Table 3. SIMBIOS Science 
Team meeting is held each year, information and 
documentation’s produced in the two previous 
meetings (1997 and 1998) are posted on the SIMBIOS 
website. 

Figure 1 shows the global distribution of the 
NRA-selected SIMBIOS field studies, and following 
chapters from 5 to 27 describe the funded research 
topics, field studies and the first year research results. 


Table 2. SIMBIOS contracts and government interagency agreements * 


Principal Investigator 

Contract 

Number 

Begin Date 

Monthly Report 

End-year TM 

Cruise Report 
Metadata 

Special Topic TM 

Field Data Within 
3 or 6 Months 

Algorithm 

Development 

Science Team 
Meeting 

William Balch 

97268 

9/23/97 

V 

✓ 

V 


V 


✓ 

J. Brock / C. Brown 

* 

7/21/97 

V 

/ 

V 


V 


V 

D. Capone /E. Carpenter 

97131 

9/11/97 

/ 

✓ 

V 


V 


V 

Ken Carder 

97137 

9/23/97 

V 


V 


V 


V 

Francisco Chavez 

Hrilitl 


■a 

m 



V 


V 

Glenn Cota 



D 


V 


V 


V 

Tom Dickey 



n 

✓ 

V 


■a 


V 

Dave Eslinger 

97133 


D 

✓ 

V 


■a 


V 

Robert Frouin 

97135 


D 

V 

✓ 


V 


V 

Robert Green 

* 


D 

V 

V 


V 


V 

Watson Gregg 



D 

V 

V 



V 

V 

Rick Miller 

■ 


D 

D 

V 


V 


V 

Greg Mitchell 

97130 


V 

n 



V 

V 

V 

Jim Mueller 

97126 


V 

✓ 

H 


mm 


V 

Frank Miiller-Karger 

97128 


✓ 

V 

■a 


m 


V 

Dave Siegel 



✓ 

V 

✓ 

D 

✓ 


V 

Ron Zaneveld 



V 

V 

✓ 

V 

V 


✓ 

Piotr Flatau 

97139 


✓ 

V 

✓ 


V 

V 

V 

Mark Miller 

HI 

HHB 

✓ 

V 

✓ 

✓ 

V 


V 

John Porter 



✓ 

V 

✓ 


V 


V 

Knut Stamnes 


9/15/97 



✓ 

✓ 

V 

V 

V 


* Interagency agreement 
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Table 3. Funded Principal Investigators (PI) and collaborators *. Collaborators are funded via EOS 
(MODI S Ocean Team) and/or SeaWiFS Project. 

I. Ocean Bio-optical and Sensor Characterization Studies 


K. Arrigo * (NASA) 

W. Balch (Bigelow Lab.) 

J. Brock (NOAA) 

C. Brown (NOAA) 

D. Capone (U. of Maryland) 

E. Carpenter (State U. of New York) 

K. Carder (U. South Florida) 

F. Chavez (Monterey Bay Aquarium Res. Inst.) 

D. Clark * (NOAA) 

G. Cota (Old Dominion U.) 

T. Dickey (U. of Califomia/Santa Barbara) 

M. He (Peoples Republic of China, Ocean U. of Qingdao) 

W. Esaias * (NASA) 

D. Eslinger (U. of Alaska) 

R. Evans *(U. of Miami) 

R. Frouin (Scripps Inst. Oceanography) 

R. Green (Jet Propulsion Lab.) 

F. Hoge *(NASA) 

S. Hooker *(NASA) 

O. Kopelevich (Soviet Union, P.P. Shirshov Inst, of Oceanology) 

G. Korotaev (Ukraine, Marine Hydrophysical Inst.) 

H. Li (Taiwan, Nat. Taiwan Ocean U.) 

R. Miller (NASA) 

G. Mitchell (Scripps Inst, of Oceanography) 

J. Mueller (San Diego State U.) 

H. Li (Taiwan, Nat. Taiwan Ocean U.) 

R. Miller (NASA) 

G. Mitchell (Scripps Inst, of Oceanography) 

J. Mueller (San Diego State U.) 

F. Miiller-Karger (U. South Florida) 

S. Saitoh (Japan, Hokkaido U.) 

D. Siegel (U. of Califomia/Santa Barbara) 

R. Zaneveld (Oregon State U.) 


IL Data Merger Studies 


W. Gregg (NASA) 

D. Siegel (U. of Califomia/Santa Barbara) 


HI Atmospheric Correction Studies 


P. Flateau (Scripps Inst, of Oceanography) 
H. Gordon * (U. of Miami) 

M. Wang (UMD/Baltimore) 

M. Miller (Brookhaven National Lab.) 

J. Porter (U. of Hawaii) 

K. Stamnes (U. of Alaska) 
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Chapter 4 


SIMBIOS Project Activities (1997-1998) 

Menghua Wang 

University of Maryland Baltimore County 
at Goddard Space Flight Center, Greenbelt, Maryland 

Brian Franz, Alice Isaacman, Christopher Pietras, Brian Schieber and Paul Smith 

SAIC General Sciences Corporation 
Beltsville, Maryland 

Tom Riley and Gene Feldman 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 


Sean Bailey 

FutureTech Corporation 
Greenbelt, Maryland 


The project accomplishments during this first 
period (1997-1998) under (a) satellite data processing, 
(b) data product validation, and (c) sensor engineering 
and calibration activities are described below. 

4.1 SATELLITE DATA 
PROCESSING 

One of the primary goals of the SIMBIOS project 
is to develop methods for meaningful comparison and 
possible merging of data products from multiple ocean 
color missions. Direct comparison of such products is 
complicated by differences in sensor characteristics 
and processing algorithms, as well as spatial and 
temporal coverage. 

4.1.1 MOS 

The German Modular Optoelectronic Scanner 
(MOS) (Zimmermann and Neumann, 1997) is an 
imaging pushbroom Charge-Coupled Device (CCD) 
spectrometer that was launched in a sun-synchronous 
polar orbit in the spring of 1996 on the Indian Remote 
Sensing (IRS-P3) satellite. 

MOS is a technology demonstrator instrument 
with limited geographic coverage capabilities. With 
the successful launch of NASA’s SeaWiFS on August 
1 of 1997, there are now two ocean color missions in 
concurrent operation. Table 1 provides characteristics 
of MOS compared with SeaWiFS. Therefore, we have 
an unprecedented opportunity to compare ocean color 


data from two sensors in simultaneous operation on 
two different satellite platforms. We have extensively 
studied the ocean color measurements compared 
between MOS and SeaWiFS and developed a 
vicarious calibration approach in which the MOS 
spectral bands can be recalibrated from the SeaWiFS 
measurements, which are considered as “truth,” 
thereby allowing remotely retrieved ocean color 
results from the two sensors to be meaningfully 
compared (Wang and Franz, 1998 and 1999). The 
vicarious calibration method has also been 
successfully applied to Japan’s OCTS data, and will 
also be tested with the French POLDER data. In here, 
we give some summaries about these studies. 

The MOS calibrated and geo-located data (LIB 
data) are provided by the MOS project at Institute of 
Space Sensor Technology, Deutshe Forshungsanstalt 
Luft-und Raumfarhrt (DLR), Berlin, Germany The 
MOS LIB data file is prepared as size of 384 x 384. 
In the following sub-sections, we give brief discuss of 
our efforts in comparing the MOS results with 
SeaWiFS. 

Atmospheric Correction and De striping 

It is well known that the atmospheric correction, 
which removes more than 90% of sensor-measured 
signals contributed from atmosphere in the visible, is 
the key procedure in ocean color imagery data 
processing. We have implemented the SeaWiFS 
atmospheric correction algorithm (Gordon and Wang, 
1994b) into the MOS imagery data processing using 
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the method outlined by Wang (1999). Therefore, a 
consistent atmospheric correction method can be 
applied to both SeaWiFS and MOS. 

The MOS radiance image has along-track stripes 
due to variations in the relative response of the in- 
dividual detectors on the MOS CCD array (total of 
384 CCD detectors). Therefore, we have developed a 
simple destriping algorithm and applied it to the MOS 
radiance imageries. The MOS destriping procedure 
can be briefly outlined as follows. First, for each scan 
(along the detector array) and a given spectral band, fit 
the radiance to a least-square cubic polynomial along 
the scan and compute relative gain at each detector. 
Next, for each detector select the median gain over all 
scans in the scene to derive the nominal gain factor for 
that detector. Finally, the MOS radiance image can be 
recomputed with the destriping correction using the 
median gain factor. 

A complicated stray light and striping correction 
algorithm is now available from the MOS project. 
Some study and tests are needed to see the efficacy of 
this algorithm. We applied the atmospheric correction 
to both MOS and SeaWiFS for co-located images and 
compared the retrieved ocean and atmospheric optical 
properties. Two MOS-SeaWiFS co-located images 
acquired on January 29, 1998 in the Atlantic ocean 
(latitude 27° and longitude -32°) and February 28, 
1998 in the Mediterranean Sea Oatitude 38° and 
longitude 3°) were first tested. These two scenes, 
acquired one month apart, differ significantly in their 
ocean and atmospheric optical properties. In compar- 
ing the MOS retrieved ocean and atmospheric optical 
results with that of SeaWiFS, however, we found that: 

• the MOS retrieved aerosol optical thickness at the 
NIR band was usually a factor of 2-3 times higher 
than that of SeaWiFS; 

• the MOS retrieved e (7,8) which characterizes the 
spectral variation of aerosol optical properties is 
unreasonably low; and 

• the MOS retrieved normalized water-leaving 
reflectances [p w 00]w in ^ visible 
significantly different from those of SeaWiFS. 

Since we are applying an identical atmospheric 
correction process to the two sets of measurements, 
the large discrepancy in the retrieved results between 
the two sensors can probably be interpreted as a 
difference in sensor calibrations. It is therefore 
necessary to recalibrate one sensor to the other, to 
allow for meaningful comparisons of the retrieved 
ocean optical properties. 

A Vicarious Intercalibration for MOS 

The vicarious intercalibration procedure between 
MOS and SeaWiFS can be described as follows. We 
assume that the gain of the MOS 868 nm band is 


unchanged because of differences in the orbits of 
MOS and SeaWiFS, thereby using the aerosol 
concentration from the MOS measurements, and only 
accept that the aerosol model determined by SeaWiFS 
is still valid. Next, by using the SeaWiFS retrieved 
aerosol models we can theoretically predict the 
atmospheric effects in the MOS imagery. The 
whitecap radiance contribution can be estimated in the 
same way as SeaWiFS (Gordon and Wang, 1994a). 
Finally, using the SeaWiFS retrieved normalize water- 
leaving reflectance, [ A* (A)] n > the water-leaving and 
total radiance at the TOA in the MOS imagery can be 
computed, and the gain coefficients for the MOS 
bands can be derived. To reduce the variation of the 
derived gain coefficients with various scans, multiple 
scans within the MOS scene can be used to obtain 
coefficient data and derive a best fit for the MOS 384 
detectors. Table 2 provides the derived MOS re- 
calibration gain coefficients fitted with the least- 
square cubic polynomial. These coefficients were 
derived from two MOS data set acquired on January 
29 and February 28, 1998. We found that the derived 
gain coefficients for the MOS bands 1-6 have very 
similar values in the two different cases, indicating 
that they are nearly independent of temporal and 
spatial variations. The derived gain coefficients for 
band 7, however, are different in the two cases. It 
appears that the MOS 750 nm band performance is 
related to the atmospheric optical conditions and its 
gain adjustment is in opposite to other bands (gain 
coefficient > 1). 

Since the derived MOS 750 nm band re- 
calibration gain coefficients depend on the 
atmospheric optical properties, we have modified the 
atmospheric correction algorithm such that the 
correction can also be operated using the MOS 685 
and 868 nm bands. Therefore, a consistent re- 
calibration gain coefficients for the MOS bands 1-6 
and band 8 can be applied. 

Results and Discussions 

We applied the derived MOS gain coefficients as 
in Table 2 to the MOS measured-radiance at the TOA, 
and retrieved ocean and atmospheric optical properties 
for comparison with results from the SeaWiFS 
measurements. Figure 1(a) and 1(b) provide examples 
of the histogram (%) for the retrieved normalize 
water-leaving reflectances (%) for bands 2 from MOS 
in comparison with SeaWiFS for the case of January 
29 and February 28, 1998. There are four cases in 
each figure: (i) results from the SeaWiFS 

measurements with the bands 7 and 8 used in the 
atmospheric corrections, (ii) results from the MOS re- 
calibrated radiances with the bands 7 and 8 used in the 
corrections, (iii) same as in (ii) except that the MOS 
bands 6 and 8 were used in the corrections, and (iv) 
results from the MOS original radiance data with the 
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bands 7 and 8 used in the corrections. Note that for 
cases in which the MOS bands 7 and 8 were used for 
the atmospheric corrections, two different calibration 
gain coefficients were applied for the MOS 750 nm 
band for cases of January 29 and February 28, 1998. 
However, when the MOS bands 6 and 8 were used in 
the corrections, the MOS band 7 reflectance data were 
simply not used, thereby allowing a consistent set of 
recalibration gain coefficients for the MOS bands 1-6 
and 8 to be applied for both cases. Figures 1(a) and 
1(b) show that the vicarious calibration improves the 
agreement significantly. 

To further test the efficacy of the vicarious re- 
calibration approach, we have applied the MOS re- 
calibration gain coefficients, which were derived from 


January 29 and February 28, 1998 data, to a MOS 
image acquired on September 24, 1997 (4-5 months 
before) at a location of about latitude 45° and 
longitude 13° in the Adriatic Sea, and compared the 
results to those obtained from a co-located SeaWiFS 
image. Again, MOS results were improved 
significantly with vicarious re-calibration (Wang and 
Franz, 1999). 

We conclude from these studies that, with the 
vicarious calibration approach, the retrieved results 
from different sensors can be meaningfully compared 
and possibly merged. More importantly, with the 
same procedure, one can re-calibrate satellite sensors 
using in situ ocean and atmospheric optical property 
measurements. 


Table 1. Characteristics of MOS with SeaWiFS. 


Instrument 

MOS 

SeaWiFS 

Platform 

I RS-P3 

Orb View- 2 

Launch date 

March 21, 96 

August 1 , 97 

Altitude (km) 

817 

705 

Equatorial crossing time 

10:30 AM 

12:00 Noon 

Resolution at nadir (km) 

0.5 

1.1 

Scan swath (km) 

200 

2800 

Time for one orbit (minutes) 

101 

99 

Spectral range (nm) 

408-1010 

412-865 

Instrument calibration 

Lamp 

Solar & Lunar 


Table 2. The derived MOS gain coefficients as G(> l,/) = (A)i" , for / = 1-384 

n= 0 


MOS 

/- (nm) 

Co<?0 

Ci(X) 

C 2 W 

CsW 

408 

0.9029 

3.5X10" 4 

-9.0xl0- 7 

1.4xl0- 9 

443 

0.8453 

3.8X10- 4 

-7.0X10- 7 

6.5xl0- 10 

485 

0.8097 

3.8X1CT 4 

-5.3xl0- 7 

2.1xl0- 10 

520 

0.8693 

1.7X10- 4 

-4.9xl0- 8 

2.8xl0- 10 

570 

0.8701 

1. 8X10- 4 

2.2xl0- 7 

-4.6xl0' 10 

685 

0.9287 

7.6x10-4 

-2.6X10- 6 

3.5xl0' 9 

750t 

1.3208 

-3,5x10-4 

-2.9xl0" 6 

7.4xl0' 9 

750t 

1.2287 

5.1x10-4 

-5.3X10- 6 

-3.2xl0- 10 

868 

1.0000 

0.0 

0.0 

0.0 


t From case of Jan. 29, 1998. t From case of Feb. 28, 1998. 
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MOS and SeaWiFS (Band 2) 
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Figure 1 . The histogram (%) of the MOS retrieved normalized water-leaving reflectances (%) with and 

without recalibrations in comparison with the SeaWiFS measurements for 443 nm band for case of (a) January 
29, 1998 (Atlantic Ocean) and (b) February 28, 1998 in the Mediterranean Sea. 
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4.1.2 OCTS 

The OCTS is an optical radiometer which flew on 
the Japanese Advanced Earth Observing Satellite 
(ADEOS) from August 1996 to June 1997, collecting 
10-months of global ocean color data (Kawamura, 
1998). During the ADEOS mission lifetime, 
approximately 450 GB of real-time, 700m-resolution 
OCTS data was collected by the SeaWiFS project 
through NOAA ground stations at Wallops, Virginia 
and Fairbanks, Alaska. The archive consists of 337 
scenes of the U. S. East Coast and 1311 scenes over 
Alaska. 

Since no standard software was publicly 
available, the coastal U S. data set was originally 
processed and distributed in near-real time using 
NASA software developed under rapid prototyping 
conditions. Over the past year, the SIMBIOS project 
has been engaged in a comprehensive program to 
enhance the algorithms and software used to process 
the OCTS data, with specific attention given to 
improving the geolocation and image registration, 
atmospheric correction algorithms, and vicarious 
calibration. In addition, through its close working 
relationship with the OCTS team at the Japanese 
Space Agency (NASDA) Earth Observation Research 
Center (EORC), the project has collected an archive of 
OCTS Level-IB and Level-2 NASDA products 
corresponding to times and locations where in situ 
data is available. This match-up data set has been 
used to validate the Level-2, Version-3 products from 
NASDA (Shimada, 1998b), as well as the Level-2 
products generated using SIMBIOS project software 
and algorithms. 

Processing to Level-0 

The U.S. Coastal data set must first be converted 
from a raw 10-bit format to a 16-bit Level-0 format. 
The original processing made use of a NASA-defined 
Level-0 format, which was inconsistent with the 
NASDA Level -0 format (NASDA, 1994). To 
eliminate this inconsistency, the SIMBIOS software 
was enhanced to generate NASDA format, including 
headers containing the definitive spacecraft 
ephemeredes and spacecraft clock timing data. The 
entire raw-data archive was then reprocessed to 
NASDA Level-0 format, and the data set was again 
made available to the public. It can be accessed via the 
"File Request and Staging" function of the SeaWiFS 
data processing system at http://falefa.gsfc.nasa.gov 
/~seawifsd/sdpsdoc/html/main. html . 

Processing to Level-IB 

In general, the processing of Level-0 data to 
Level- IB includes conversion of raw sensor counts to 
physical units and the assignment of geolocation 


information to each observation (pixel, detector, and 
band). Several factors in the design of the OCTS 
instrument complicate this process. Each OCTS band 
is divided into 10 individual detectors, where each 
detector is associated with a separate scan line Each 
detector has a slightly different and apparently non- 
linear responsively, which gives rise to horizontal 
striping in the OCTS imagery. The basis for the 
SIMBIOS calibration is the NASDA Version 3 
preflight calibration and relative detector calibration 
(Shimada, 1998a). The relative detector calibration 
substantially reduces the striping effects, but low-level 
artifacts are still visible and they can become 
significant following removal of the large, 
comparatively smooth atmospheric signal. 

The design of the OCTS scan mirror and focal 
plane results in severe misalignment of the line-of-site 
for different spectral bands viewing at essentially the 
same time. To perform atmospheric correction or 
derive chlorophyll concentration at a given location 
requires that the bands be co-registered to that 
location. 

The SIMBIOS project developed a registration 
technique in which an idealized focal-plane geometry 
is defined, and all bands are re-sampled using a 
nearest neighbor approach to match the idealized 
ground track. This method is similar in concept to the 
NASDA approach (NASDA, 1994), but it will result 
in slightly different re-sampling. Furthermore, while 
the registration process can correct for differences in 
geolocation between bands, the re-sampled 
observations come from different lines-of-sight with 
unique sensor-to-ground and sun-to-ground 
geometries. 

To ensure accurate atmospheric correction, the 
true solar and sensor view angles per band are saved 
in the Level-IB file for use in Level-2 processing 
Before the registration process can begin, however, 
the observations from each band must be accurately 
navigated. The geolocation algorithm used for 
SIMBIOS OCTS processing is a modified version of 
the method used for SeaWiFS navigation (Patt and 
Gregg, 1994). This technique has the potential to yield 
an exact solution to the geolocation of each scan line, 
when the spacecraft position and sensor geometry are 
accurately known. 

Some sources of error include: inexact knowledge 
of the instrument design and installed, post-launch 
sensor and mirror alignments (NASDA, 1994), 
uncertainties in orbit location or spacecraft attitude, 
and variations in the scan-rate or sensor tilt 

To improve geolocation accuracy, the SIMBIOS 
project derived time-dependent attitude and tilt 
adjustments using an automated navigation 
assessment technique based on island targets (Patt et 
al., 1997). The final navigation is still being assessed, 
but preliminary indications show that it is accurate to 
approximately one kilometer. 
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Processing to Level-2 

Within the SIMBIOS project, Level-IB to Level- 
2 processing includes application of the vicarious 
calibration, removal of the atmospheric signal to 
retrieve the water-leaving radiances, and generation of 
derived products such as chlorophyll concentration. It 
is most appropriate to apply the vicarious calibration 
at this stage, as it must account for deficiencies in the 
atmospheric correction algorithm, as well as for 
changes in the preflight calibration. Since this same 
approach was used by NASD A, the Level- IB products 
from the Japanese Space Agency and the SIMBIOS 
project are calibration equivalent, though the 
atmospheric correction algorithms and file formats 
differ. 

The SIMBIOS OCTS atmospheric correction 
algorithm is based on the SeaWiFS algorithm of 
Gordon and Wang (1994a). The SIMBIOS approach 
has been to identify those few parts of the algorithm 
which are sensor or band-pass specific, and develop a 
software package which can process data from 
multiple ocean color remote sensors. The sensitivity of 
the algorithm to differences in spectral bands was 
carefully assessed by Wang (Wang, 1999), wherein he 
showed by simulation that these differences can be 
accurately accounted for through exact calculation of 
the Rayleigh reflectances and minor modifications to 
the diffuse transmittance calculation. The multi-sensor 
Level-IB to Level-2 code (MSI 12) has been used 
successfully to process and intercalibrate MOS and 
SeaWiFS data (Wang and Franz, 1999), and it is now 
being employed to process the entire SIMBIOS OCTS 
archive. 

When processing OCTS data, the MSI 12 software 
is able to make use of the band-dependent solar and 
view angles. If these angles are provided in the Level- 
IB input file, they will be used in the atmospheric 
correction process when computing the Rayleigh 
radiances, the air-mass calculations that affect 
estimates of white-cap radiance and diffuse 
transmittance, and the normalization of the water- 
leaving radiances. 

The atmospheric correction process includes 
masking for land, clouds, and saturation and masking 
for stray-light in the vicinity of those bright sources. 
To mitigate the effects of residual striping or other 
forms of systematic or random noise, the MSI 12 code 
includes options for various types of statistical, spatial 
filtering. For OCTS processing, a 5x5 -pixel median 
filter is applied to the Rayleigh-subtracted radiances in 
the near-IR bands, based on the assumption that the 
spatial scale of aerosol variability is expected to be 
greater than a few kilometers. This smoothing 
effectively reduces the tendency for the atmospheric 
correction algorithm to select different aerosol types 
from pixel-to-pixel. Additional smoothing with a 3x3- 


pixel median on the Rayleigh-subtracted radiances is 
applied to the visible bands, reducing the effect of 
residual striping at the cost of degraded resolution of 
the in-water spatial structure. 

For purposes of vicarious calibration, the 
SIMBIOS project is fortunate to have time and space- 
coincident in situ water-leaving radiance spectra from 
the Marine Optical Buoy (MOBY) buoy (Clark et al., 
1997) to match a series of OCTS over-flights. This 
same buoy is the basis for the SeaWiFS calibration, 
providing a convenient bridge to connect the two 
mission life-spans. Unfortunately, there were some 
problems with the early deployment of MOBY, 
limiting the number of usable calibration samples to 
seven measurements spanning the period November 
1996 to February 1997. In addition, MOBY 
measurements are not sufficient to calibrate the near- 
IR bands at 765 nm and 865 nm, where water-leaving 
radiances are very close to zero. It was therefore 
assumed that the 865 nm vicarious calibration was 
equal to the value suggested in the NASDA Version-4 
calibration (Fukushima, pers. comm ), and the 765 nm 
band was adjusted to force the atmospheric correction 
algorithm to select, on average, a maritime aerosol 
type over the MOBY site. This assumption is 
consistent with the typical aerosol-type retrievals 
derived from SeaWiFS data in the vicinity of MOBY. 
The entire vicarious calibration process can be 
performed in a single processing step in which the 
OCTS match-up data set is corrected to retrieve the 
normalized water-leaving radiances, the retrieved 
normalized water-leaving radiances are replaced with 
the in situ values, and the atmospheric correction 
algorithm is reversed to derive the theoretical top-of- 
atmosphere radiances. These theoretical radiances are 
divided by the observed radiances to derive correction 
factors for each in situ match-up point, and those 
ratios are averaged within each band to derive the 
vicarious calibration coefficients. 

The final calibration coefficients are listed in the 
Table 3, with the equivalent NASDA Version-4 
calibration shown for comparison. Considering that 
the two projects use different atmospheric correction 
algorithms ( Fukushima et al., 1998) and different in 
situ measurements for calibration, the results are in 
good agreement. 


Table 3. 


Band 

NASDA 

SIMBIOS 

1 

1.14 

1.13 

2 

1.03 

1.01 

3 

0.9394 

0.94 

4 

1.00 

1.00 

5 

1.04 

1.03 

6 

LOO 

0.99 

7 

1.02 

0.91 

8 

0.89 

0.89 


1 1 ! 
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The largest difference is in the 765 nm band, 
which NASDA does not use for atmospheric 
correction (Shimada, 1998). Since SIMBIOS does 
make use of the 765 nm band for aerosol model 
selection, a correction for oxygen absorption was 
necessarily performed, which accounts for the 
discrepancy. 

Validation Analyses 

An extensive set of in situ data taken during the 
time of the OCTS mission was made available to the 
SIMBIOS project for match-up analysis from the 
SeaBASS database. Over 600 data points taken on 
134 separate days exist in this set. 

OCTS Satellite data used in the match-up analysis 
was obtained from two different sources and two 
different processing methods were used. Files in the 
U.S. Coastal archive which matched in situ data points 
were processed from Level 0 to Level 2 using the 
SIMBIOS algorithms and calibration described above. 
Level IB data files supplied to the SIMBIOS Project 
by NASDA/EORC were processed to Level 2 using 
the SIMBIOS atmospheric correction algorithms and 
vicarious calibration. Level 2 data supplied by 
NASDA/EORC were also matched to in situ data , but 
these data were processed using Version 3 of the 
NASDA calibration and so results obtained are of 
historic interest only. 

Normalized water-leaving radiance (nLw) values 
were extracted from the Level 2 records and 
acceptable pixels within a 1.05-km radius were 
weighted by the inverse of their distance to the match- 
up site and then averaged to derive nLw values for 
comparison to in situ data. This technique usually 
resulted in obtaining data from 9 pixels in a square 
centered about the in situ data point, although the 
actual number can vary due to the complex scan 
geometry of the OCTS instrument. Match-up values 
were also calculated for an area within a 2-km radius 
if there were not enough acceptable pixels within the 
1.05-km limit. 

The in situ water-leaving radiances were 
normalized using theoretical estimates of surface 
irradiance (E,). The E s values were calculated from 
the Solar radiance model of Neckel and Labs (Neckel 
and Labs, 1984), weighted by the OCTS band passes, 
and transmitted to the surface using measurements of 
the local ozone concentration, knowledge of the 
Rayleigh optical thickness, and estimates of the 
aerosol transmittance. The resulting nLw values were 
then compared to the weighted average nLw values 
from OCTS. 

Where possible, the match-ups of satellite to in 
situ data were judged acceptable by the same criteria 
used to judge SeaWiFS match-ups to in situ data. 
Specifically, match-ups were rejected if: 


• The time difference between satellite and in situ 
measurements was more than 4 hours (6 hours for 
chlorophyll). 

• Difference in solar zenith angle was more than 
15° 

• In situ data were from a duplicate cast and had a 
lower K490 value than other casts. 

• Wavelengths of in situ data differed from OCTS 
wavelengths by more than 5 nm. 

• Satellite data coefficient of determination (the 
ratio of standard deviation to average) was greater 
than 0.5 

• The usable number of satellite pixels was less 
than 5, after tests for clouds and sunglint. 

For the data processed from NASDA Level IB to 
Level 2 using SIMBIOS methods, 13 match-up points 
were deemed acceptable. An additional 10 match-up 
points were obtained from the Wallops Level 0 
dataset. Data collected over the MOBY site were 
excluded from this analysis, since they were used to 
derive the vicarious calibration. 


Table 4. 


Band 

Number of 
Match-ups 

nLw ratio 

OCTS retrieval / in situ 

i 

18 

1.20 

2 

18 

1.06 

3 

19 

1.01 

4 

13 

0.998 

5 

9 

1.04 


The Table 4 summarizes the match-up results 
between in situ water-leaving radiance data (excluding 
MOBY) and OCTS retrievals processed using the 
SIMBIOS atmospheric correction algorithms. 

Sargasso Sea Validation Study 

Due to the small number of usable match-up in 
situ data in existence for the OCTS mission, additional 
validation of the SIMBIOS OCTS calibration was 
sought. A time series of averaged normalized water- 
leaving radiances derived from 27 scenes centered in 
the Sargasso Sea was generated and compared to 
nominal water-leaving radiances derived from MOBY 
in situ data and to a time series from one year later in 
the Sargasso Sea area created from SeaWiFS data. 

The Sargasso Sea area was chosen for the OCTS 
data because it yielded the largest number of dear- 
water scenes within the OCTS data set acquired by 
the Wallops receiving station. In addition, six dear- 
water scenes centered over Bermuda were generated 
from the same archive and compared to in situ data 
from that area. 

The SeaWiFS Sargasso time-series data was 
taken from the standard 8-day time -binned dataset, for 
a 3° square box centered at Latitude 25.5° and 
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Longitude -7L5° from 57 8-day periods between 
November 1997 and December 1998. 

The SeaWiFS time binning procedure averages 
all acceptable values over each time period on a pixel- 
by-pixel basis. The nominal dear-water radiances 
used for comparison were derived from MOBY in situ 
data tuned to the OCTS and SeaWiFS wavelengths 
respectively. In situ chlorophyll and water leaving 
radiances for the Bermuda site were taken from the 
Bermuda Atlantic Time Series (BATS) conducted by 
the Center for Remote Sensing and Environmental 
Optics at the University of California, Santa Barbara 
(Garverand Siegel, 1997). 

For each scene in the Sargasso Sea and Bermuda 
area, normalized water-leaving radiances were 
averaged, excluding values where chlorophyll was not 
within acceptable limits, i.e. data not between .0005 
pg/L and .30 pg/L, as well as data of unacceptable 
quality, e g., clouds, etc. 

The averaged Sargasso Sea nLw values were 
compared to nominal MOBY-derived OCTS nLw 
values and to SeaWiFS values from 1997-1998 while 
the Bermuda area nLw values were compared to 
nominal MOBY-derived nLw values and to 
contemporary in situ data from the BATS experiment. 
For both Sargasso Sea and Bermuda area data, 
chlorophyll values were computed from the average 
nLw using both the OC2 algorithm of O'Reilly et al. 
(1998) adjusted to the OCTS spectral bands and the 
standard NASDA OCTS-C algorithm (Kishino et al., 
1998). In addition, pigment concentrations were 
computed using the Normalized Difference Pigment 
Index (Frouin, personal comm ). 

The agreement in the general behavior of the 
OCTS data with the SeaWiFS data from one year later 
suggests that the variability seen in the Sargasso Sea 
nLw values is due to seasonal variation rather than 
instrument instability. The mean nLw values from 
OCTS agree reasonably well with MOBY-derived 
nominal nLw values, suggesting that the SIMBIOS 
calibration is working well. 

The three chlorophyll algorithms used on OCTS 
data agree in the average but the distribution of the 
chlorophyll derived by the OC2 algorithm is 
somewhat broader than that obtained using the 
NASDA chlorophyll or NDPI pigment algorithm. The 
results of this study were presented at the Alps ’99 
Conference in Meribel, France in January, 1999. 

4.2 PRODUCT VALIDATION 

The Project is comparing and validating 
SeaWiFS satellite measurements with in situ 
atmospheric and bio-optical measurements. 


4.2.2 SeaWiFS AOT Comparisons 

The theoretical basis of the SeaWiFS aerosol 
retrieval algorithm are described in both Gordon and 
Wang, (1994b) and Wang et al., (1999). In this 
section, we briefly outline our efforts in comparing 
and validating the SeaWiFS aerosol optical products 
with the in situ measurements mainly from the data of 
the Aerosol Robotic Network (AERONET) (Holben et 
al., 1998). The majority of the instruments within the 
AERONET network are in continental locations The 
SIMBIOS Project augumented this network with 
coastal and island stations. Some other in situ 
measurements from calibration and validation 
campaigns within the SIMBIOS project are also 
analyzed. The primary objectives for these 
comparisons are to: 

• validate the SeaWiFS aerosol optical thickness 

(AOT) products, and 

• determine the validity of the suite of aerosol 

models currently used by SeaWiFS for 

atmospheric correction. 

In Situ Data Acquisition 

The ground-based measurements utilized for the 
aerosol optical thickness matchup analyses come from 
two primary sources: automated CIMEL sun/sky 
scanning radiometer managed as part of the 
AERONET network, and hand-held MicroTops II 
sunphotometers. A select group of the ground stations 
from the AERONET was chosen. These instruments 
were located at either coastal or island stations and 
were operational for a reasonable length of time after 
SeaWiFS went into operation. Table 1 provides the 
AERONET station name, location (latitude and 
longitude), and corresponding responsible AERONET 
principle investigator (PI), the hand-held MicroTops 
II sunphotometer data, on the other hand, were 
collected by various investigators in field campaigns 
associated with the SIMBIOS Project. Data from the 
MicroTops instruments are reprocessed from raw 
voltages using code adapted from the AERONET 
standard CIMEL processing code. This ensures that 
the T a data derived from the MicroTops 

measurements will be comparable to the data provided 
by the AERONET. 

SeaWiFS Data Acquisition 

The SeaWiFS aerosol optical thickness data were 
obtained by spatially co-locating a 5x5 pixel grid box 
around the pixel containing the ground-base 
measurement station, thereby providing a maximum 
25 SeaWiFS retrievals in each matchup. The 
SeaWiFS operational code has been modified to 
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output, at a pixel by pixel level, the aerosol optical 
thickness at wavelength 865 nm, values of retrieved 
two aerosol models, as well as the model partition 
ratio r a value. Therefore, aerosol optical thicknesses 
at all the SeaWiFS wavelengths can be calculated 
(Wangetal., 1999). 

Data Analyses 

SeaWiFS data were obtained by averaging over a 
5x5 pixel grid box spatially, whereas the T a data from 
the CIMEL measurements were derived by a time- 
weighted averaging during the SeaWiFS overpass (± 1 
hour). Usually, the CIMEL instruments routinely take 
one measurement every 15 minutes. Therefore, for a 
given SeaWiFS file there may be as many as eight 
AERONET measurements that qualify as a match for 
the 2-hour time window. The number of hand-held 
MicroTops II measurements that match a given 
SeaWiFS file, however, varies greatly since the 
measurement protocol for these instruments is not yet 
well-defined. In general, there should be a minimum 
of three MicroTops measurements per matched 
SeaWiFS file. The ground-based measurements are 
averaged after weighting by the time difference 
between the in situ measurement and the SeaWiFS 
overpass. Once averaged, the ground based 
measurements are compared with the SeaWiFS 
derived values on a band by band basis for each 
ground station. Due to a limited number of 
MicroTops data, these data are grouped together as 
one “station.” 


Preliminary Results 

We compared the SeaWiFS derived aerosol 
optical thicknesses with those from the ground in situ 
measurements. Figure 2 provides an overall 
comparison results of T a (X) between SeaWiFS and 
CIMEL measurements at wavelength 865 nm (870 nm 
for CIMEL). The CIMEL measurements were from 
the AERONET stations listed in Table 5. The dotted 
lines in Figure 2 is the 1:1 line. Though the 
comparison results varying both in time and location, 
it appears that SeaWiFS has tendency of 
overestimating T a {X) with respect to the in situ 

measurements. Note that, however, since the 
SeaWiFS band 8 has not been absolutely calibrated on 
orbit, any error in calibration may contribute to the 
error in the r a (865) evaluations. On the other hand, 
some in situ data are suspected to be erroneous due to 
instrument calibration. Obviously, more studies are 
needed to understand all of these. We want to 
emphasize that all results are preliminary. 

Similarly, the in situ MicroTops II data, which 
were from the various SIMBIOS calibration and 
validation campaigns, have been compared with the 
SeaWiFS measurements. Although this work is still 
in the initial phase, some results are promising. Table 
6 shows three sample comparison results from three 
field experiments. In these three examples, the 
SeaWiFS results were almost all underestimated as 
compared with the MicroTops II measurements (Ratio 
< 1), though the three results usually agreed 
reasonable well. 
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Figure 2. The retrival SeaWiFS aerosol optical thickness r a (/ 1) compared with the ground in situ 
measurements from AERONET at 865 nm (870 nm for CIMEL Data). 
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Table 5. AERONET sites utilized for the aerosol matchup analyses. 


AERONET Station 

Latitude 

Longitude 

AERONET PI 

Bahrain 

26.32 

50.50 

Charles McClain* 

Bermuda 

32.37 

-64.70 

Brent Holben 

Dry Tortugas 

24.60 

-82.80 

Ken Voss/Howard Gordon 

Kaashidhoo 

4.97 

73.47 

Brent Holben 

Lanai 

20.83 

-156.99 

Charles McClain* 

San Nicolas Island 

33.26 

-119.49 

Robert Frouin 


* SIMBIOS Project Office 


Table 6. Three samples of MicroTops II data compared with SeaWiFS 


Aerosol Optical Thickness T a (X) 

X (nm) 

Jason Project 

GOCAL97-98 

JULNAN98 



G. Feldman 

J. Mueller 

B. Schieber and 




A. Subramaniam 


SeaWiFS Ratio 

SeaWiFS Ratio 

SeaWiFS 

Ratio 

440 


0.1635 0.8678 

0.1775 

0.9168 



0.1469 0.9767 

0.1538 

0.9031 

670 

0.0395 0.8857 

0.1076 1.1811 

0.1013 

0.9217 

865 

0.0304 0.9102 

0.0774 0.8113 

0.0641 

0.8456 


4.2.3 SeaBASS INTERFACE 

The SeaWiFS Bio-optical Archive and Storage 
System (SeaBASS) was created by the SeaWiFS 
Project to provide a data archive for in situ products 
used in scientific analysis. This system has been 
expanded to additionally contain data sets of interest 
to the SIMBIOS Project. 

The in situ data in SeaBASS include 
measurements of water-leaving radiance and other 
related optical and pigment measurements, from ships, 
moorings and drifters. Various methods are deployed 
in the collection of SeaBASS data, including the use 
of standard profiles, floating radiometers and above- 
water measurement devices. In addition to the optical 
and pigment products listed, measurements of total 
suspended matter (TSM), chromophoric dissolved 
organic matter (CDOM), and other typically non- 
plankton derived optical components are recorded in 
SeaBASS. 

The products are used by the SIMBIOS Project in 
validation of SeaWiFS and other (OCTS, POLDER, 
etc.) postlaunch imagery and development of new 
operational chlorophyll algorithms. Additional uses 
under SIMBIOS will likely include field data 
calibration comparisons, algorithm workshops, and 


time series studies. A detailed description of the 
original SeaBASS design and scope is provided in the 
SeaWiFS Technical Report Series, specifically 
volume 20. A current description of the SeaBASS 
system is available via the World Wide Web at 
http : / / seabass . gsfc . nasa . gov . 

Data Design 

The SeaBASS system contains data from well 
over 150 separate experiments with thousands of files 
requiring over 22 megabytes of storage space. In 
addition, the historical optics database holds over 
300,000 records related to phytoplankton pigments 
and is expected to grow as new pigment data is 
collected. Since the SeaBASS system is continuously 
growing with new data, a method for efficiently 
ingesting and storing field data was needed It was 
deemed important that the data ingest be as 
straightforward and effortless as possible on the part 
of the contributing investigators, while still offering a 
useful and correct data archive for analysis efforts. 
The following aspects were considered to be the most 
important in the design of the system: 

• Simple data format 

• Global portability 
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• Data update simplicity 

• Web accessible data holdings 

• Support for existing Web standards 

To accomplish the above goals, the SeaBASS 
system was designed to support standard ASCII files 
which can be managed from any computer platform. 
The files use a simple descriptive header format which 
contains descriptor=value pairs to explain the specifics 
(“metadata”) of a data file. For example: 

/beginheader 

/affiliations = University _of_Maryland 
/investigators = John_Creamer 
/contact=iohnc@umd . edu 
/parameters=Ed,Lu,Kd,CHL 
/north_latitude=59. 10[DEG] 

/south_latitude= 58.65 [DEG] 

/east_longitude=-13. 1 1 [DEG] 

/west_longitude= -13.50 [DEG] 

! Comments: data file to be updated 

/end_header@ 

The header or metadata section is intended to 
fully describe the data in a file. Each header contains 
information on time and location of the in situ 
measurements, the investigators involved, values 
measured, comments, and other related descriptive 
information. Frequently, SeaBASS files will contain 
references to other files with additional descriptive 
details (e.g., README files). 

Following the header information, the actual 
measured optics or pigment data is listed in columnar 
format. Typical data files include optical profiles, 
summary pigment files, or along track measurements. 
Specific information and examples of the formats can 
be referenced at the SeaBASS web site. 

Following an series of redesign modifications in 
the first year of the SIMBIOS effort, the SeaBASS file 
formats have not changed significantly. However, 
additional changes may need to be made in the near 
future as new instrumentation and collection methods 
continue to populate the archive, some of which may 
warrant variations in the current data archive approach 
(i.e., hyper-spectral data). 

Format Checking 

Early in the development of the SeaBASS web 
site it was realized that each data provider can have 
very different views on data format or presentation. In 
order to establish a standard format for SeaBASS 
incoming data files which met the dual goals of 
simplicity and portability, feedback software was 
developed. 

The primary component of this software is the 
FCHECK script. FCHECK is written in the PERL 


scripting language with connections to look-up tables 
and UNIX mail handling utilities. Data providers 
from any platform can test their data for compatibility 
with the SeaBASS format by simply e-mailing their 
file to fcheck@seabass.gsfc.nasa.gov. The FCHECK 
program will receive the incoming mail file from the 
system mail handler, parse the data, and thoroughly 
compare it to the established SeaBASS format. 
FCHECK will then report back to the data provider on 
the file’s compliance. This file testing function 
happens automatically with no intervention required, 
but the SeaBASS administrator is sent a copy of all 
incoming data files, for review. The FCHECK 
program has been used by many SeaBASS data 
providers as a quick format checker and has saved 
many hours of processing time for both the SeaBASS 
administrator and contributing investigators. 

Once the data provider is sure that their files will 
work with the SeaBASS archive, they can send the 
data, and related documents, to the SeaBASS archive 
via the file transfer protocol (FTP). From there, the 
SeaBASS administrator will move the files into their 
appropriate position in the archive and will initiate a 
database update to register the new data to the 
SeaBASS online system. 

Data Security Issues 

A concern of some contributors to SeaBASS 
involves outside user accessibility. Some investigators 
have been concerned that when data is quickly 
submitted to the archive, others may use the data 
without proper notification or acknowledgement. To 
address these concerns, and to afford continued rapid 
submission of data sets, the SeaBASS web server has 
been configured as a password protected system 
Typically, only SeaWiFS or SIMBIOS authorized 
investigators who also submit data to the archive are 
allowed access. Others may be allowed access on a 
case-by-case basis. Additionally, the web server and 
SeaBASS software log all user activity. This 
information is available to contributing investigators. 
It is expected that all users who use data from the 
SeaBASS archive will offer collaboration or 
authorship to the data provider. 

Specific Matchup Format 

To expedite comparisons between in situ data and 
SeaWiFS (or other satellite) products, a specific file 
format has been developed. This file format, known as 
the “matchup” format, follows a similar logic to that 
of the other SeaBASS in situ files, that being a 
simple, easy to implement format in standard ASCII 

However, unlike the general format for SeaBASS 
data, the matchup format expects a strict adherence in 
the positions and units of data fields. This requirement 
is to allow for automated comparisons between in situ 
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by phytoplankton) water characterization. Additional 
investigators are also supported to develop new 
algorithms or scientific approaches in accordance with 
the goals of the SIMBIOS Project. As described 
above, the data and information provided are used for 
both internal analysis and satellite validation and to 
populate the SeaBASS archive for future activities and 
community use. Table 7 summarizes the data sets 
which have been received from investigators during 
the first year of funding. As shown in the table, the 
majority of funded principal investigators provide 
field data for the archive. However, some additional 
investigators were funded solely to develop algorithms 
or analysis techniques. 

Future Plans 

The SIMBIOS project has created an environment 
by which in situ data may be archived, accessed and 
analyzed to improve satellite retrievals through 
scientific data validation. In its first year of operation, 
the Project has implemented software and format 
standards which should greatly aid in future satellite 
efforts. Additionally, the Project’s focus on thorough 
data documentation should aid in both current and 
future research activities. Subsequent years will focus 
on improving the usefulness of in situ data by 
continuing to improve the data ingest and data 
analysis approach. 


Table 7. SIMBIOS supported data submitted to SeaBASS 


Principle 

Investigator 

University 
or Laboratory 

Optics 

Data 

Pigment 

Data 

Documentation 

Provided 

mmm 

William Balch 

Bigelow Labs 

V 

V 

v 


J. Brock / C. Brown 

USGS/NOAA 

V 

V 

V 


Douglas Capone 

U. MD 

V 

V 

V 


Ken Carder 

U. South Florida 

V 

V 

V 


Francisco Chavez 

Mont. Bay ARI 

V 


V 


Glenn Cota 

ODU 

V 

V 

V 


Tom Dickey 

UC SantaBarbara 

V 


V 


Dave Eslinger 

U. AK, NOAA/CSC 

V 

V 

V 


Piotr Flatau 

Scripps, UCSD 



V 

V 

Robert Frouin 

UCSD / Scripps 

V 


V 


Robert Green 

JPL 

V 


V 


Watson Gregg 

GSFC 



V 

V 

Mark Miller 

BNL 

V 


V 


Rick Miller 

Stennis Space Ctr. 

V 




Greg Mitchell 

Scripps, UCSD 

V 

V 

V 

V 

Jim Mueller 

San Diego St. U. 

V 

V 

V 


Frank Mtiller-Karger 

U. South Florida 

V 

V 

V 


John Porter 

U. Hawaii 

V 


V 


Dave Siegel 

UC Santa Barbara 

V 

V 

V 

V 

Knut Stamnes 

U. Alaska 

V 


V 

V 

Ron Zaneveld 

osu 

V 

V 

V 
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and satellite data. The format is also used in 
interactions with the SeaWiFS/SIMBIOS satellite data 
archive. The matchup format is structured in standard 
ASCII columnar format with the following fields in 
each data record: 

Year, Month, Day, Hour, Minute, Second, Latitude, 
Longitude, 1^412, Lw 443, Lw 490, Lw 510, Lw 555, 
U 670, E.412, E 9 443, E s 490, E s 510, E* 555, E 9 670, 
Kd490, Chi a 

Where the time and location values identify when 
and where the measurement was taken. Water-leaving 
radiance (Lw) and downwelling irradiance (E s ) values 
may be extrapolated from below the surface or 
measured above the surface, depending upon the 
instrumentation used. The vertical attenuation of 
irradiance at 490 nm (Kd 490) as well as in situ 
measured values of chlorophyll a (Chi a) are also 
recorded if available. Missing data values are 
represented by a placeholder. As in the case of other 
SeaBASS submitted data, the matchup file also 
contains a metadata section to describe the enclosed in 
situ data. 

Data Contributions 

The SIMBIOS project funds numerous 
investigators in order to obtain in situ optical and 
pigment data for Case I (optically characterized solely 
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4.2.4 Support Services 

In an effort to improve the quality and quantity of 
calibration and validation data sets, the SIMBIOS 
Project offers several support services to field 
investigators. These services include; scheduling of 
on-board Local Area Coverage (LAC) recording for 
SeaWiFS; overflight predictions for operational 
sensors (currently SeaWiFS, MOS-B and OCI); near 
real time SeaWiFS imagery for cruise locations; and 
optical instrumentation from a pool of investigator- 
and project-owned instruments. These services may 
be requested via the World Wide Web at 
http://simbios.gsfc . nasa. gov. 

In return for these services, the SIMBIOS Project 
requests that the field investigators provide in situ 
validation data to the Project's bio-optical archive, 
SeaBASS. Since October of 1997, when these 
services were initially offered, the SIMBIOS Project 
has supported 81 cruises (Table 8). 

SeaWiFS Scheduling 

Since much of the worlds oceans are not covered 
by a SeaWiFS High Resolution Picture Transmission 
(HRPT) station, high-resolution data may be recorded 
onboard the SeaWiFS sensor. As a service to the 
science community, the SIMBIOS project in 
conjunction with the SeaWiFS project can schedule 
SeaWiFS onboard LAC for cruises that occur outside 
HRPT coverage. SeaWiFS has the ability to record a 
maximum of 10 minutes of high-resolution data per 
downlink. Typically, a 30-second interval is allotted 
for LAC target, which corresponds to 180 scan lines 
or approximately 200 km along track at nadir. 
Detailed information on LAC scheduling is available 
on the SIMBIOS web site. 


Overflight Predictions 

For calibration and validation purposes, in situ 
measurements should be made as close to the sensor 
overflight time as is possible. To aid investigators in 
determining when sampling should occur, the 
SIMBIOS project offers overflight predictions for all 
operational ocean color remote sensors. Currently, 
the sensors supported are SeaWiFS and MOS-B 
With the successful launch of ROCS AT- 1 on 26 
January 1999, OCI is being added to this list. Detailed 
information on overflight predictions is available on 
the SIMBIOS web site. 

Near-Real Time SeaWiFS Imagery 

In addition to providing predictions for satellite 
overflight times, the SIMBIOS project offers near real 
time imagery of the operations SeaWiFS products in 
JPEG format to cruises at sea. These images provide 
field investigators with additional information with 
which they may maximize in situ sampling of 
transient oceanographic features. The default 
specifications for the images provided include: 

• Available LAC, HRPT, and Global Area 
Coverage (GAC); 

• Chlorophyll-a and pseudo-true color images; 

• 2-degree box about a designated location or the 
entire designated region; 

• Image Width: 600 pixels; and 

• Minimum percent valid chlorophyll pixels: 5% . 

Images may be customized to best accommodate 
individual investigators needs. Detailed information 
on near real time imagery is available on the 
SIMBIOS web site. 


Table 8. Table of SIMBIOS supported cruises with services provided. 


Cruise Location 

Begin Date 

End Date 

0 

5 

1 

2 

1 

c 

o 

Over-flight 

Prediction 

Near Real Time 
Imagery 

Instrument Pool 

Principal Investigator 

Chesapeake Bay 

09/28/1997 

11/23/1997 

V 




Curt Davis 

Black Sea 

10/03/1997 

10/13/1997 



g 

■ 

Oleg Kopelevich 

Alaska 

10/10/1997 

10/17/1997 


■ 

| 

H 

.cj r -] ii i i*T»i 

Gulf Of Maine 

11/10/1997 

11/21/1997 


U 

u 

■ 


Taiwan 

11/23/1997 

11/25/1997 

•j 

V 



Hslen-wen Li 

Equatorial Pacific 

09/27/1997 

12/16/1997 

/ 

/ 



Francisco Chavez 

Ross Sea 

10/01/1997 


D 


D 



Equatorial Atlantic 

11/14/1997 

12/12/1997 


✓ 

✓ 


C. Menkes 

HOTS 

12/01/1997 

12/21/1997 

•j 




John Porter 
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Begin Date 


12 / 07/1997 


10 / 19/1997 


01 / 09/1998 


01 / 31/1998 


01 / 26/1998 


01 / 25/1998 


02 / 10/1998 


03 / 02/1998 






02 / 15/1 998 


01 / 12/1998 

02 / 11/1998 


02 / 12/1998 


02 / 07/1998 


03 / 27/1 998 


03 / 16/1998 
99 



East China Sea 


Bahamas 


HOTS 


Sicily Channel 
alifornia Coast 



MMEM 


m 


Equato 

rial Pacific 

Subtroi 

Dical North Pacific 

HOTS 


Atlantic 

) Ocean - Baltic Sea 

Gulf Of 

Cadiz 


Baltic Sea 

Atlantic Meridional Transect 
Gulf Of Maine 
Equatorial Pacific 
Massachusetts Bay 
Labrador Sea 
East China Sea 
South China Sea 


Gulf Of Mexico 


Arctic Ocean 
Resolute Bay 


Equatorial Atlantic 
Long Island Sound 
Puget Sound 


Puget Sound 
Arctic Ocean 
Puget Sound 
Taiwan Straits 


02 / 28/1998 


03 / 16/1998 


04 / 1/1998 


04 / 14/1998 


03 / 27/1998 


04 / 02/1998 


04 / 03/1998 


04 / 17/1998 


03 / 24/1998 


04 / 6/1998 


04 / 6/1998 


04 / 26/1998 


04 / 21/1998 


998 

998 




05 / 14/1 998 


06 / 01/1 998 


06 / 6/1998 


07 / 05/1 998 


06 / 22/1 998 


06 / 27/1998 


06 / 28/1998 


06 / 27/1998 



Instrument Pool 

Principal Investigator 


Marlon Lewis 

V 

Greg Mitchell 


John Porter 


Marlon Lewis 

V 

Dennis Clark 


Eurico D'Sa 

V 

Greg Mitchell 


Marcel Wernand 


Jim Mueller 


ohn Morrison 


03 / 27/1998 


04 / 07/1998 




imiai 


V Robert Steward 


M. Ribera d'Alcala' 


04 / 24/1998 




V Sonia Gallegos 

V V Gene Feldman 


04 / 30/1998 


99 

05 / 01/1998 


05 / 2/1998 


05 / 06/1998 


mpgBEESM 



05 / 19/1998 


06 / 03/1998 


06 / 18/1 998 


06 / 15/1998 


06 / 16/1998 


07 / 10/1 998 


07 / 10/1 998 


arcel Wemand 


Stan Hooker 


Peter Strutton 


effrey Polovina 


ohn Porter 


Roland Doerffer 




Peter Land 


Stan Hooker 




— amai 


998 1 08 / 25/1 998 
99 



07 / 20/1998 
06 / 04/1998 06 / 25/1998 V 
998 " 08 / 22/1 998 7 “ 

08 / 05 / 1998 “ 08 / 07/1 998 
07 / 24/1998 08 / 22/1 998 
08 / 10/1998 08 / 20/1 998 


V V Peter Strutton 


jit Subramaniam 


lenn Cota 


Gwo-Ching Gong 


Hsien-Wen Li 


Dariusz Stramski 


Glenn Cota 


C. Menkes 


S. Papanikolaou 


V Ron Zaneveld 


/ Ron Zaneveld 


Sergey Shirshov 
M.J. Perry 
^Ming-Xia HE 
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Cruise Location 

Begin Date 

End Date 

On-board LAC 

Over-flight 

Prediction 

Near Real Time 
Imagery 

Instrument Pool 

Principal Investigator 

HOTS 

08/08/1998 

08/12/1998 


✓ 



John Porter 

Gulf of Maine 

08/14/1998 

08/18/1998 



✓ 


Bruce Bowler 

SE Bering Sea 

08/15/1998 

09/09/1998 



V 


Stephan Zeeman 

South West U.S. 

08/17/1 998 

09/25/1998 

v 

✓ 

y/ 


Edward Zalewski 

South China Sea 

09/11/1 998 

09/27/1998 

s/ 

/ 


/ 

Rick Miller 

California Coast 

09/13/1998 

10/02/1998 


/ 


V 

Gregg Mitchell 

North Pacific 

08/30/1998 

09/21/1998 

y/ 

V 


V 

Ron Zaneveld 

HOTS 

09/26/1998 

09/30/1998 

V 

V 



John Porter 

North Pacific 

10/01/1998 

10/12/1998 

y/ 

y/ 

✓ 


Carrie Leonard 

Atlantic Meridional Transect 

9/14/1998 

10/30/1998 

y/ 

y/ 



Stan Hooker 

Barents Sea 

9/15/1998 

10/10/1998 

/ 

/ 



Oleg Kopelevich ! 

U.S. East Coast 

10/5/1998 

10/13/1998 


/ 



Stephan Howden 

Gulf Of Maine 

9/10/1998 

10/31/1998 

V 

V 

✓ 

V 

Bruce Bowler 

HOTS 

10/18/1998 

10/20/1 998 


y/ 



John Porter 

East China Sea 

10/15/1998 

10/31/1 998 

/ 

/ 

y/ 


Gwo-Ching Gong 

Equatorial Pacific 

10/23/1 998 

11/09/1 998 

v 

y/ 

yj 


Peter Strutton 

HOTS 

11/10/1 998 

11/12/1998 


/ 


~ 

John Porter 

South Atlantic Bight 

10/27/1998 

11/24/1998 


/ 

✓ 


Ajit Subramaniam 

Gulf of California 

11/25/1998 

12/09/1998 

✓ 

V 

V 

y/ 

Andrew Barnard 

Gulf of Maine 

12/02/1998 

12/18/1998 


~T 

y/ 


Heidi Sosik 

HOTS 

12/13/1998 

12/16/1998 

/ 

“7“ 



John Porter 

Indian Ocean 

01/04/1999 

02/23/1999 

/ 


✓ 


J. Le Fevre 

Indian Ocean - INDOEX 

01/14/1999 

04/01/1999 

/ 


/ 


Piotr Flatau 

Southern Ocean 

01/01/1 999 

04/01/1999 

/ 


yj 


Steve Groom 

California Coast 

01/09/1999 

01/31/1999 


y 

v' 


Greg Mitchell 


Field Instrument Pool 

The SIMBIOS project provided funding to 
several of the science team members for the purchase 
of in situ ocean optical instrumentation. The funding 
was provided with the stipulation that these 
instruments would be made available for three years to 
an instrument pool to be maintained by the Project 
Office. The Project augmented this instrument pool 
with atmospheric instrumentation. Table 10 
summarizes the instruments available. 

4.3 SENSOR ENGINEERING AND 
CALIBRATION 

AERONET is a network of ground-based 
automated sun photometers owned by national 
agencies and universities. AERONET data provides 
globally distributed, near-real time observations of 


aerosol spectral optical depths, aerosol size 
distributions, and precipitable water. These data allow 
for algorithm validation of satellite aerosol retrievals 
as well as characterization of aerosol properties. Such 
validation and characterization are critical for 
atmospheric correction of ocean color sensors. The 
majority of the instruments within the AERONET 
network are in continental locations. The CIMEL Sun 
Photometer manufactured in France is the instrument 
used for AERONET. 

The StMBIOS Project augmented his instrument 
pool with atmospheric instruments by purchasing 12 
MicroTops hand-held sun photometers, 2 PREDE 
(Japanese built) sun photometers, one micro-pulse 
lidar and the 12 CE3 18 CIMEL sun photometers. 

Specifically, the 12 CIMEL instruments augment 
the AERONET network with coastal and island 
stations. These CIMEL instruments better withstand 
the corrosive marine environment after undergoing a 
robust re-engineering. These “hardened” CIMEL 
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instruments are placed, with other instruments, at sites 
near marine environments. Confirmed (delivered or 
negotiated) CIMEL sites include Lanai-Hawaii (with a 
backup in Honolulu), Ascension Island, Bahraine, 
Tahiti, Wallops Island (Virginia), South Korea, and 
Turkey (Black Sea). 

CIMEL Re-engineering 

In its normal configuration, the CIMEL Sun 
Photometer is not suitable for use in a corrosive 
marine environment. Meridian Engineering of 
Hanover (Maiyland) under Tom Riley’s guidance 
hardened these instruments for placement in locations 
where they may be subject to salt spray. The 
modifications include: 

• Replacement of the base housing — The new 
base housing is welded and has more weather 
resistant connectors. 

• Replace motor housing — The new housing is 
welded. 

• Replace instrument housing — Again the new 
housing is welded. 

• Cable and connectors — Connectors were 
replaced with military style connectors. Cables 
replaced with ones rated for this environment. 
The cables were then sealed with both 3M 2242 
Linerless Electrical Rubber Splicing Tape and 
then 3M Scotch Super 88 Vinyl Electrical Tape. 

• Hardened moisture sensor — The moisture sensor 
cable and connector were replaced. 

The following devices were added to allow 
remote operation: 

• Electronic boxes — The two electrical boxes 
(Control Box and Power Box) are now sealed 
plastic NEMA boxes and are painted white. 

• Photovoltaic panel mount — The two 
photovoltaic panels and the moisture sensor are 
mounted on a fiberglass panel that is held to the 
power box with stainless steel hardware. 

• Sun Shield — The Control Box also has a fiber 
glass sun shield held one inch above the box and 
has stainless steel hardware. 

• Mounting Plate — A heavy galvanized steel 
mounting plate is now available. 

Sun Photometer Calibration: 

Non-Polarized Channels 

The calibration consists of using a sun photometer 
calibrated in high altitude conditions (Mauna Loa, 
Hawaii) and transferring the calibration to the other 
sun photometers by taking simultaneous sun direct 
measurements from ground during a clear day. The 


top of atmosphere (TOA) signals (VO) of the non- 
calibrated sun photometer are deduced from those of 
the CIMEL sun photometer. This method is currently 
applied to all the CIMEL sun photometers involved in 
the AERONET network (Holben et al., 1998). 

The processing steps for analyzing the calibration 
measurements performed by a sun photometer are 
described below and summarized in Diagram 1 : 

• Download the raw signal of the uncalibrated sun 
photometer [Vu (i)] and of the calibrated sun 
photometer [Vr (i)]. 

• Select the simultaneous measurements (time 
difference lower than 36 seconds). 

• Compute for each channel, the TOA signals for 
the un-calibrated sun photometer using the TOA 
signals of the calibrated sun photometer (obtained 
at Mauna Loa). 

• The TOA signals, computed for each 
simultaneous event, are averaged 

The comparison with the aerosol optical thickness 
(AOT) retrieved from the calibrated sun photometer 
allows for validation of the protocol. AOT are 
retrieved using the same steps as in Diagram 2 except 
for the Ozone retrieval. Total Ozone Mapping 
Spectrometer (TOMS) measurements are used for the 
computation of ozone contribution. 

Currently, eight MicroTops, one SIMBAD and 
two PREDE sun photometers have been successfully 
transfer calibrated according to the inter-calibration 
method. The TOA signals of the MicroTops #3773, 
SIMBAD #972396 and PREDE #PS090063 are 
summarized during several transfers of calibration 
between August and December 1998 in Table 1 
Means and standard deviation are presented for each 
channel of each instrument. The standard deviation of 
the TOA signal is overall lower than 2%. The results 
are especially good for the MicroTops and SIMBAD 
sun photometers, considering the long period used, the 
design (i.e., hand held) and affordable cost The 
accuracy of the CIMEL sun photometer TOA signals 
(calibrated at Mauna Loa and used for our calibration 
transfer) is 0.5 % for the channels between 440 and 
1020nm (Holben et al., 1998), The accuracy of the 
calibration transfer for the PREDE PS090063 is only 
presented on Table II. Another PREDE (#PS090064) 
instrument, used between October and December 
1998, is still under study. 

Sun Photometer Calibration: 

Polarized Channels 

The polarized version of the CIMEL sun 
photometers is composed of three channels at 870nm 
mounted on the filter wheel. Three polarizing sheets 
are located in front of the filters and are respectively 
adjusted according to 0, -60 and +60 degrees. 
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The polarization rate of the light is deduced from three 
measurements in these three channels. 

The calibration of the polarized version of the 
CIMEL sun photometer is performed using a 
calibration device designed by the Laboratoire 
d’Optique Atmospherique. The calibration device is 
located ahead of a non-polarized source and allows the 
polarization of the light in the output of the device. 
The reading of the angular position of the two 
SFllwindows mounted inside the device allows for 
the computation of the polarization rate of the light 
The calibration consists of plotting the computed 
polarization rate versus the measured one and 
determining the slope and the intercept of the trend 
line. 

The calibration have been performed on six 
polarized instruments belonged to the AERONET 
group. The results are summarized in Table 12. The 
same calibration is scheduled for the SIMBIOS 
polarized CIMEL sun photometers. The calibration is 
performed for a polarization rate ranging between 0 to 
65%, because of the dimensions of the device. The 
linear squared fit determined for the calibration of the 
six CIMEL sun photometers ranges between 0.9993 to 
1 . 

Sun Photometer Data Processing 

The data processing is applied to all SIMBIOS 
sun photometers (i.e., CIMEL, PREDE, MicroTops II 
and SIMBAD). The data processing principle is 
described below and summarized in Diagram 2: 

1. Download the raw signal V, corresponding to the 
measurements performed in the channel i. 

2. Compute the ratio between the raw signal and the 
top of Atmosphere (TOA) signal V,o provided by 
the calibration of the sun photometer, V/V, 0 . 

3. Compute the ratio between the averaged Sun- 
Earth distance and the Sun-Earth distance of the 
day, (do/d). 

4. Compute the air mass m, according to the day, the 
date, the time and the location when/where the 
measurements were taken. 

5. Compute the total optical thickness in the channel 

i, Ttot(i) 

6. According to the location, compute and remove 
from the total optical thickness the Rayleigh and 
the ozone optical depth for each channel, r R (i) 
and x 0 (i). The same equation and Dobson table is 
used for the data processing of all the sun 
photometers. 

7. Compute the aerosol optical thickness (AOT), 
t a(0- 

Several AOT measurements have been performed 
on the roof of a building at GSFC (from August to 
December 1998), from cruises during several 


campaigns and from ground, allowed performance 
comparisons with SeaWiFS sensor (Wang et al., 
1999). 

Calibration Round Robin 

NASA personnel carried out the first SeaWiFS- 
SIMBIOS Intercalibration Round-Robin Experiment 
(SIRREX-6) from August 1997 to February 1998 
SIRREX-6 was performed in a completely different 
manner from SIRREX-1 through SIRREX-5. In those 
tests, laboratory references were brought to a central 
location and tested against one another. In SIRREX-6, 
the same four common field instruments (eg., 
Satlantic in-water radiometers) were taken to nine 
separate laboratories and tested using the laboratories’ 
standards and procedures. Two of the sensors were 
seven-channel radiance heads and two were seven- 
channel irradiance heads. The calibration and data 
reductions procedures used at each site followed that 
laboratories’ normal procedures. The reference lamps 
normally used for the calibration of these types of 
instruments by the various laboratories were also used 
for this experiment. NASA personnel processed the 
data to produce calibration parameters from the 
various laboratories for comparison. These tests 
showed an overall agreement at better than the +/-2% 
level. Test equipment, calibration procedures, data 
reduction, and specific handling procedures have been 
published in the Riley and Bailey (1998) 

SXR-II 

The first SeaWiFS Transfer Radiometer (SXR) 
was built for the SeaWiFS Project to verify' and 
compare measurements of spectral radiance at six 
discrete wave-lengths in the visible and near infrared 
Johnson et al. (1998). In addition SXR is used to 
compare these sources to standards of spectral 
radiance maintained at the National Institute of 
Standards (NIST). SIMBIOS project had a second 
copy of the SeaWiFS Transfer Radiometer (SXR-II) 
built for use in the calibration round robin. This unit 
will supplement the first unit and is designed for 
easier travel. The SXR-II is currently being 
characterized and calibrated. It will be used in future 
SIRREX round robin experiments. 

Additional SQMs 

NASA developed a portable optical calibration 
source called the SeaWiFS Quality Monitor (SQM) 
with NIST. Last year, NASA provided limited 
financial support for the development of the SQM as a 
commercial product. Two companies elected to take 
advantage of the opportunity, Satlantic Instruments 
and Yankee Environmental Systems, Inc. The 
Satlantic version was evaluated by Dr. Stan Hooker. 


27 



SIMBIOS Project Annual Report 


The Yankee Environmental Systems, Optical 
Calibration Source prototype was evaluated at GSFC 
by Tom Riley. Their computer interface is a 
significant improvement over the earlier models, A 
number of small problems were uncovered and 
reported back to the manufacter who has taken 
corrective action. 

4.4 SIMBIOS COMPUTING 
RESOURCES 

The SIMBIOS computing facility contains many 
different types of equipment (Diagram 3). The main 
goal of the facility is to provide support for the 
following areas: 1) routine bulk data processing, 2) 
large scale analysis and match-up between sensors, 3) 
focused analysis on data subsets, and 4) storage of 
large sensor data sets. 

To support these efforts, the facility has grown to 
include two large server systems, a 6 terabyte tape 
library, four desktop workstations, and numerous 
personal computers and printers. Co-location of the 
facility with the SeaWiFS Project has also made 
available other resources such as high-quality color 
printers, network equipment, and other supporting 
equipment. Large scale tasks such as bulk data 
processing and large scale analysis tasks are 
performed by two Silicon Graphics (SGI) Origin 2000 
servers. To prepare for increased workload, the 


systems were upgraded to 6 processors each, and 
approximately 3 gigabytes (GB) of RAM. The influx 
of large data sets for additional sensors requires 
increasing amounts of space, especially when 
performing comparative analysis. The Data Processing 
server has been increased to 98 GB of disk, and the 
Data Analysis server has over 372 GB of total disk 
space. Long term storage of data sets is handled by 
the 6 terabyte (TB) tape library attached to the Data 
Analysis server. 

Silicon Graphics 02 workstations are used for 
software development and more focused analysis 
tasks. The sizing of the workstations varies with the 
task requirements. One workstation has 320 
megabytes (MB) of RAM and 28 GB of disk space, a 
second workstation has 640 MB of RAM combined 
with 72 GB of disk space. The third workstation has 
640 MB of RAM and 50 GB of disk space. 
Acquisition of a fourth workstation is planned for this 
fiscal year. 

SIMBIOS is connected to the Internet via a 100 
megabit link, making network transfer of large data 
sets more feasible. Internal connections are also 100 
megabit for rapid transfer of data between the servers 
and workstations. Data ingest for SIMBIOS is also 
done via removable media. The Project can read a 
wide variety of removable media including CD-ROM, 
8 millimeter tape, 4 millimeter (DAT) tape, and DLT 
media up through DLT 7000. 


Table 10, SIMBIOS Pool Instruments 


Instrument 


Description 

Manufacturer 

Custodian(s) 

MicroTops II 

Sunphotometer 

w/GPS 

12 

5 channel handheld sun photometer 
w/Garmin GPS-38 

Solar Light 
Company 

SIMBIOS Project 

HI STAR Package 

1 

2 hyper-spectral absorption/beam attenuation 
meters mounted on a cage with a SeaBird 
CTD 

WETLabs, Inc 

Zaneveld/Pegau 

Hydroscat 6 

3 

Backscattering meter 

HOBILabs 

Siegel 

Carder 

Mitchell 

Pure Water System 

3 

Water purification system for calibration of 
WETLabs AC-9 and HISTAR absorption and 
attenuation meters 

Barnstead 

Zaneveld/Pegau 

AC-9 

3 

Absorption/beam attenuation meter 

WETLabs, Inc. 

Cota 

Capone 

Muller-Karger 

SIMBAD 

2 

5 channel radiometer and sun photometer 


Frouin 

SeaWiFS 
Multichannel 
Profiling Radiometer 

2 

Free-fall profiling radiometer measuring L u (z), 
E d (z), E, , U(1 m) with cables, deck unit and 
PC with data acquisition and processing 
software 

Satlantic, Inc. 

Capone 

Chavez 

Micropulse LIDAR 

1 _j 

Continuous operation LIDAR system 

SESI 


Prede Sun 
Photometer 

2 

Automated ship-board sun photometer 

Prede 

SIMBIOS Project 

CE 318 
CIMEL Sun 
Photometer 

12 

j 

Marine-hardened, automated sun photometer 

Cimel 

SIMBIOS Project 
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Table 11. Top of Atmosphere signals (V 0 ) of MicroTops, Simbad and Prede sun photometers determined by 

transfer calibration from a calibrated C1MEL at GSFC between August and December 1998. 


MicroTops #3773 

440nm 

500nm 

675nm 

870nm 

940nm 

08/20/98 - Cimel # 37 

1236.95i0.78 

980.63il.37 

1214.27i0.98 

821.80il.12 

1419.70il.22 

08/21/98 - Cimel # 37 

1243.64il.92 

■ 

1220.24i0.86 

825.59i0.76 


10/16/98 - Cimel # 27 

1216.46±0.91 


1187.49i0.77 

822.38i0.58 

1400.10i0.55 

1 1/24/98 - Cimel # 27 

1219.75i0.45 

976.43i0.46 

1191.42i0.34 

821.64±0.61 

1401.36i0.75 

MEAN 

1229.20 

979.14 

1203.35 

822.85 

1410.43 

STD % 

1.07 

0.79 

1.35 

0.22 

0.79 

Simbad #972306 






08/20/98 - Cimel # 37 

3 93 564. 40i 1.89 

4 82 820. 40i 1.07 

403245. 16i0.35 

420247.88i0.84 

307582.92i0.59 

08/21/98 - Cimel # 37 

393525. 05±3. 13 

481855. 71il.02 

407705. 35i0.41 

423453. 93il. 03 

305498.45i0.56 

10/16/98 - Cimel # 27 

381894. 63il. 25 

469676. 76i0.99 

393682.49il.82 

409544. 15i0.56 

306692.23i0.98 

1 1/24/98 - Cimel # 27 

392974. 50±0.89 

477633. 99i0.97 



307275. 49i0. 17 

12/14/98 -Cimel# 101 

386234. 51±0. 95 

471653.55i0.37 

393092.41±0.21 


310487.83i0.22 

MEAN 

389638.62 

476728.08 

401421.16 

415116.63 


STD % 

1.36 

1.24 

1.91 

1.54 


Prede #PS090063 






08/20/98 - Cimel # 37 

1.52E-04±3.38 

2.86E-04il.43 

3.59E-04±2,22 

2.81E-04i0.36 


08/21/98 -Cimel #37 

1.48E-04±1.13 

2.84E-04il.45 

3.57E-04il.75 

2.81E-04i0.45 


MEAN 

1.50E-04 

2.82E-04 

3.58E-04 

2.81E-04 

2.88E-04 

STD % 

1.48 

0.51 

0.27 

0.11 

1.87 


Table 12: Slope (5p), intercept (p 0 ) and linear squared fit applied to the polarization rate given by the device 

and plotted versus the polarization rate measured by six polarized CIMEL sun photometers. 


Instrument # 

__ 5 p 

Po 


92 

1.1088 

-0.1051 


25 

1.2256 

-0.0523 

■iTOlrirM 

111 

1.0791 

0.0476 

1 

45 

1.0608x 

0.2354 


43 

1.1935X 

-0.7816 

1 

14 

1.1665 

-11.179 

0.9993 
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Diagram 1. Sun photometer transfer calibration 



Diagram 1 Sun photometer data processing 


Raw signals (in V) or (in A) in the channel i: Vi 
Date, Time, Longitude, Latitude 


Calibration Input 
Top of Atmosphere signal ViO 


Earth-Sun distance computation 
according to the averaged distance 
dO/d 

(Iqbal M., 1983) 

alrmass computation 
m 

(Karten et al, 1989) 

i _ 1 



Total optical thickness computation 
x tot(l)=Ln [(Vi rVi 0)/(d 0/d) A 2]/m 


1 - 

Rayleigh optical thickness computation 
xR(i) 

(Penndorf R., 1957) 

Ozone optical thickness computation 
xO{i) 

Lookup table by 5degrees of latitude 
(London et al, 1976) 


J 


Rayleigh and Ozone correction 
tA (i)=Ttot(i)-TR{i)-TO(i) 


Aerosol optical thickness computation 
tAO) 
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Diagram 3. SIMBIOS Analysis/Data Processing System 
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Chapter 5 

Validation of Surface Bio-Optical Properties 
in the Gulf of Maine as a Means for Improving 
Satellite Primary Production Estimates 

William M. Balch 

Bigelow Laboratory for Ocean Sciences 
West Boothbay Harbor, Maine 


5.1 INTRODUCTION 

There is a strong need to have sea-truth optical 
data in order to evaluate remotely-sensed ocean color. 
Our SIMBIOS contract is to use the M/S Scotian 
Prince ferry as a ship of opportunity, running between 
Portland, Maine and Yarmouth (NV). While the Gulf 
of Maine is cloudy and foggy more than it is clear ; the 
climatology shows on average, about 1 in 4 days is 
completely clear and clear days are slightly more fre- 
quent in die late summer and early fall months. A 
ship of opportunity program (where one has choice of 
the days of sampling) provides much better flexibility 
to sample during clear periods with good satellite cov- 
erage. Measurements include continuous, surface, 
along-track fluorescence, light scattering, absorption, 
beam attenuation, above- water remote sensing reflec- 
tance, calcite-dependent light scattering, temperature, 
and salinity. Expendable Bathythermograph (XBT) 
drops allow acquisition of vertical temperature infor- 
mation, useful for defining isopycnal slope, which 
affects primary production. These data are compara- 
ble to a previous program from early 1982, where a 
Ship Of Opportunity Program (SOOP) was run on the 
truck ferry, M/V Marine Evangeline, which ran along 
the same transect (Boyd, 1985). These surface data 
were combined with satellite-derived sea surface tem- 
perature fields to examine the Maine coastal current 
(Bisagni et al., 1995). Unfortunately, this program 
stopped in 1982. The ongoing SIMBIOS results will 
dovetail nicely with the previous work (which also 
had CZCS coverage) for looking at any long-term 
changes in the Gulf of Maine hydrography or bio- 
optics. 

The Gulf of Maine is an ideal site for gathering 
sea-truth optical data. It is a semi-enclosed basin with 
an average depth of about 150m. Isopycnal slope (and 
productivity) appears to be strongly topographically 
driven, with the water column becoming isothermal at 
the 60m isobath around Georges Bank. Integrated 
productivity in the region varies seasonally by a factor 

of 45x (0.1 to 4.4 gC m* 2 d' 1 ; O'Reilly, 1987), greater 
than in the California Current (15x variability; 0. 2-3.0 


gC m" 2 cT 1 ). The Marine resources Monitoring As- 
sessment and Prediction (MARMAP) database along 
the eastern seaboard of the U S has been the equiva- 
lent to the California Cooperative Oceanic Fisheries 
Investigation (CalCOFI) database (West Coast) in 
terms of its sheer size and length of sampling How- 
ever, MARMAP does not maintain pigment or inher- 
ent optical property data. This SIMBIOS project at- 
tempts to provide this unique data for the region. 

5.2 RESEARCH ACTIVITIES 

We participated in several field cruises to field 
test equipment, and collect data. We have completed 
the following field activities: 

• June '98-Second set of sea trials of sampling and 
underway system in the Gulf of Maine and Geor- 
ges Bank (lid cruise). 

• September-October '98-Eleven cruises aboard 
M/S Scotia Prince completed. 

Real time enumeration of coccolithophores is es- 
sential to this work (as well as cell counts done after 
the cruise). We have an Olympus BH2 microscope 
with epi-fluorescence and polarization optics which 
allows us to quantify birefringent CaC03 coccoliths, 
empty coccospheres, and fluorescing, plated coccoli- 
thophore cells on board the ship. Due to limited time 
at each station, only about two samples typically can 
be enumerated. However, counts are recorded on 
video, which provides an important record for post- 
cruise comparison to the preserved counts, as well as 
for any retrospective analyses. We take preserved cell 
count samples at all stations (preserving with 4% 
buffered formalin) and, after the cruise, measure the 
concentration of detached coccoliths and plated coc- 
colithophores (Utermbhl, 1931 and 1958). 

The technique of Fernandez et al. (1993) is used 
to measure CaC03 concentrations. Briefly, samples 
are filtered onto 0.4(im polycarbonate filters, and 
rinsing first with filtered sea water, then borate buffer 
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(pH=8) to remove seawater calcium chloride. Filters 
are placed in trace metal free centrifuge tubes with 5 
ml 0.5% Optima grade Nitric acid. Next, the Ca con- 
centration is measured using graphite furnace atomic 
absorption spectrometry. The sensitivity of the tech- 
nique, after correction to the volume of seawater fil- 
tered, is about 2ng Cal" 1 . 

Particulate organic carbon and nitrogen are meas- 
ured according to Sharp (1974; 1991). One liter sam- 
ple bottles will be emptied via tygon tubing, through 
in-line pre-combusted glass fiber filters. Low vacuum 
will be maintained during filtration. These filters will 
be placed in plastic petri dishes using clean forceps, 
and stored at -20°C. Prior to running the samples, they 
will be fumed over HC1 to remove any carbonates, 
then dried in a dessicator, and shipped on dry ice to 
the Bermuda Biological Laboratory for analysis. 
Chlorophyll is measured fluorometrically according to 
Strickland and Parsons (1972) and modified according 
to Joint Global Ocean Flux Study (JGOFS) (1996) 
which involves filtering 200 ml seawater with Gelman 

GFF filters, extracting overnight at 4°C in 10 ml of 
90% acetone, and measuring fluorescence with a 
Turner 1 1 1 fluorometer. 

We locate high densities of coccolithophores in 
the field with a flow-through system which estimates 
suspended CaC 03 using an optical technique. A 
Wyatt Technologies laser-light scattering photometer 
(equipped with a flow-through cell) is used to measure 
optical volume scattering at 18 angles. Integration of 
this signal in the backward direction allows calcula- 
tion of backscattering in real time. A pump injects 
weak acid every 4 minutes to dissolve the CaC 03 , and 
then backscattering is re-measured. The difference 
between the total bb and the acidified bb is called the 
“acid-labile bb,” and is calibrated to suspended 
CaC03 concentrations. The acidified b^ represents 
backscattering due to organic matter (which can be 
calibrated to CHN measurements). This value can be 
converted to POC concentration using our POC- 
specific backscattering coefficients (Balch et al., 
1999). The system monitors chlorophyll fluorescence, 
pH, temperature, salinity and 18-channel volume 
scattering (with and without CaC 03 ). An AC-9 pro- 
vides estimates of absorption and attenuation, which 
by difference, also provides scattering. A Global Po- 
sitioning System is interfaced, as well. The flow- 
through system has been field-tested over carbonate 
banks off of Florida, which provided data over 3 or- 
ders of magnitude in acid-labile b^. Such high con- 
centrations of suspended CaC 03 are found in the most 
dense coccolithophore blooms that we have ever vis- 
ited. During the Arabian Sea expedition, as well as 5 
recent Gulf of Maine cruises, the instrument ran virtu- 
ally flawlessly over thousands of kilometers, provid- 
ing a wealth of data on the distribution of particulate 


inorganic and organic carbon (Balch et al., 1999). The 
final cruise of this year was completed in October; the 
raw data have been worked up and distributed to the 
NASA SIMBIOS Project. Discrete samples for sus- 
pended calcite, POC, PON, and coccolith counts were 
also taken, and will be processed over the coming 
months. 

5.3 RESEARCH RESULTS 

We completed 33 days at sea in 1998, with 11 
trips aboard the M/S Scotia Prince ferry (Table 1). 
Two of 11 trips were cloudy (82% were sunny) In 
the Gulf of Maine, where 28% of the time it is sunny, 
we have demonstrated that we take full advantage of 
ferry as a ship of opportunity. For the year, we have 
processed about 147 stations in the Gulf of Maine and 
Georges Bank, and provided continuous underway 
data along 5000 miles of ship track for SeaWiFS cali- 
bration. Below are summarized the main sampling 
activities: 

• R/V Delaware 10-21 November 1997 R/V Dela- 
ware, Gulf of Maine and Georges Bank Under- 
way data analysis completed and submitted to Se- 
aBASS database. Data collected: underway fluo- 
rescence, backscattering and calcite-dependent 
backscattering (510 nm; Wyatt light scattering 
photometer, absolutely calibrated with a glass 
standard), absorption, scattering (AC-9; SeaWiFS 
bands; recently calibrated), temperature, salinity 
Discrete measurements at 27 stations of sus- 
pended calcite, POC/PON, chlorophyll, coccolith 
and coccolithophore concentration. Discrete cell 
counts and analytical analyses are still being 
processed. 

• R/V Albatross 3-12 June 1998, Gulf of Maine and 
Georges Bank. Data collected: underway fluores- 
cence, backscattering and calcite-dependent back- 
scattering (510 nm), absorption, scattering (AC-9: 
SeaWiFS bands), temperature, salinity, L w and Ed 
at SeaWiFS wavelengths measured with Satlantic 
SAS system. Discrete measurements at 28 sta- 
tions of suspended calcite, POC/PON, chloro- 
phyll, coccolith and coccolithophore concentra- 
tion. 

• M/S Scotia Prince on a total of 11 cruises (Table 
1). Data collected: underway fluorescence, back- 
scattering and calcite-dependent backscattering 
(510 nm; Wyatt light scattering photometer, ab- 
solutely calibrated with a glass standard), absorp- 
tion, scattering (AC-9; SeaWiFS bands; recently 
calibrated), temperature, salinity. L w and E d at 
SeaWiFS wavelengths measured with Satlantic 
SAS system. Microtops sun photometer measur- 
ments performed for overpasses. Discrete meas- 
urements at stations of suspended calcite, 
POC/PON, chlorophyll, coccolith and coccolitho- 
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phore concentration. Discrete cell counts and 
analytical analyses being processed. 

5.4 WORK PLAN 

Our contract calls for 20 ferry trips next year. We 
anticipate no problem completing these trips since our 


laboratory container and sample arm are already com- 
pleted and ready. 

Moreover, we anticipate no problems in timely 
data work-up and submission as our processing soft- 
ware is now streamlined for post-cruise calibrations. 


Table 1. 


Summary of cruises and in situ data collected. 


Delaware, Gulf of Maine and 
George Bank 


[Gulf of Maine and George Bank 


Portland, Maine and 
[Yarmouth (NV) 


Cruise Location 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Portland, Maine and 
Yarmouth (NV) 


Cruise 

Name 


lOd 


lid 


Ferry 1 


Ferry 2 


Ferry 3 
Ferry 4 


Ferry 5 


Ferry 6 


Ferry 7 


Ferry 8 


Ferry 9 


Ferry 10 


Ferry 11 


Date 


10-21 
Nov 1997 


3-1 2 

June 1998 


10-11 

Sept. 1998 
16-17 
Sept. 1998 


17-18 
Sept. 1998 


28-29 
Sept. 1998 


3-4 

Oct. 1998 


4-5 

Oct. 1998 


11-12 
Oct. 1998 


12-13 
Oct. 1998 


16-17 
Oct. 1998 


17-18 
Oct. 1998 


18-19 
Oct. 1998 
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Chapter 6 

OCTS and SeaWiFS Bio-Optical Algorithm and 
Product Validation and Intercomparison in 
U.S. Coastal Waters 

John C. Brock 

USGS Center for Coastal Geology and Regional Marine Studies 
St. Petersburg, Florida 

Christopher W. Brown 
NOAA/NESDIS, Camp Springs, Maryland 


6.1 INTRODUCTION 

The launch of the National Space Development 
Agency of Japan (NASDA) Ocean Color and Tem- 
perature Sensor (OCTS) in August 1996, and the 
launch of Orbital Science Corporation’s (OSC) 
SeaWiFS in August 1997 signaled the beginning of a 
new era for ocean color research and application. The 
ocean color data provided by these satellites promises 
the ability to remotely evaluate and monitor 1) water 
quality, 2) transport of sediments and adhered pollut- 
ants, 3) primary production, upon which commercial 
shellfish and finfish populations depend for food, and 
4) harmful algal blooms which pose a threat to public 
health and economies of affected areas. Accordingly, 
several US government agencies have recently ex- 
pressed interest in using optical remote sensing in US 
coastal waters for these purposes. 

The goal of this project is to evaluate the per- 
formance of standard bio-optical algorithms and vali- 
date satellite ocean color products derived from OCTS 
and SeaWiFS data in Case I and Case n coastal US 
waters. To accomplish this goal, in situ bio-optical 
observations were collected and analyzed. This 
document briefly describes the methods used in sam- 
ple collection and analysis, and presents preliminary 
results obtained during the first year of the project. As 
a result of the failure of the Japanese satellite on June 
30, 1997, and the consequent termination of the OCTS 
data stream, we have been unable to evaluate the per- 
formance of OCTS algorithms and data products. 

6.2 RESEARCH ACTIVITIES 

In-situ bio-optical data were collected during 
eight cruises in diverse US Case I and Case II waters 
along the Atlantic coast and in the Great Lakes since 
OCTS became operational in November 1996 (Table 
1). Optical instruments were deployed on cruises to 


measure surface spectral downwelling irradiance, in- 
water spectral downwelling irradiance and upwelling 
radiance. Deployment of other instruments and the 
collection of bottle samples enabled additional meas- 
urements that included temperature, chlorophyll fluo- 
rescence, light scattering, quantum scalar irradiance, 
total suspended solids concentration, chlorophyll and 
pigment concentrations from fluorometric and High 
Pressure Liquid Chromatography (HPLC) techniques, 
and absorption coefficients of colored dissolved or- 
ganic matter and particles. Although sampling 
strategies and instrument packages varied between 
cruises, a Biospherical Instruments Profiling Reflec- 
tance Radiometer (PRR) cage was typically deployed 
off the stem of the boat in conjunction with a refer- 
ence surface unit with matching channels. The PRR 
cage contains a split PRR600s that measures seven 
channels of downwelling irradiance, seven channels of 
upwelling radiance, depth, tilt, roll, and temperature. 
A reference surface unit PRR610 that measures seven 
matched channels of surface downwelling irradiance 
was also used. PRR600s channels 1 to 6 are narrow 
band (10-nanometer [nm], full width half maximum 
[TWHM]) centered at 380 nm, 412 nm, 443 nm, 490 
nm, 510 nm, and 555 nm, while channel 7 on the 
downwelling sensor and PRR610 measures broad 
band Photosynthetically Available Radiation (PAR) 
(400 to 700 nm). In addition, the PRR cage contains a 
10-centimeter pathlength, 660 nm Light Emitting Di- 
ode (LED) -based SeaTech transmissometer, a Bio- 
spherical Instruments Quantum Scalar Profiling sensor 
(QSP200), a SeaTech light scattering sensor (LSS), 
and a WETLabs Wetstar chlorophyll fluorometer. 
Water is drawn through the fluorometer using a Sea- 
Bird pump running at 2,000 revolutions per minute. 
The data from these instruments are multiplexed 
through the PRR600s such that each record contains a 
depth and the parameters from every instrument. 

In some cases, a Hydro-Optics, Biology, and In- 
strumentation Laboratories, Inc. (HOBI Labs) Hydro- 
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Scat-6 spectral backscattering sensor was lowered by 
hand following the deployment of the PRR cage. 
Typically, an along-track system was used to measure 
the position (latitude, longitude), time, course and 
speed of the vessel, temperature, and salinity. In situ 
temperature, salinity, and density were also measured 
at some stations with a Conductivity-Temperature- 
Depth (CTD) instrument Water samples for chloro- 
phyll biomass, particulate, and dissolved absorption 
were obtained from a separate cast using a Niskin ar- 
ray. Discrete water samples were collected following 
the PRR cast, from just below the sea surface using a 
Niskin bottle. 

Bio-Optical Data Processing 

The PRR optical data was processed using the 
Bermuda Bio-Optics Project (BBOP) processing soft- 
ware (Siegel et al, 1995c). The data were separated 
into upcast and downcast profiles and then binned to 
0.5 meter bins. Spectral attenuation coefficients were 
calculated for the optical channels over a five point 
moving window. Subsurface downwelling irradiance 
and upwelling radiance were extrapolated to just be- 
low the surface using data from the top 3 meters. The 
HydroScat-6 data were processed using software pro- 
vided by HOBI Labs, Inc. The data were despiked, in 
two passes using a difference threshold, and a moving 
average was calculated for these channels. These data 
were also separated into upcast and downcast profiles 
and then binned to 0.5-m bins. 

For fluorometric chlorophyll a and HPLC pig- 
ment determinations, water samples are filtered 
through glass fiber. The chlorophyll a samples are 
cold extracted in 10 ml of 90% acetone (10% water) 
for 24 hours in the dark, and the biomass is deter- 
mined fluorometrically with a Turner Designs fluo- 
rometer using the method of Yentsch and Menzel 
(1963). Filtered HPLC pigment samples are cut into 
small pieces, ground, and extracted in a 1.5 ml micro- 
centrifuge tube with 1.5 ml 90% acetone and are 
placed in a freezer overnight. They are centrifuged at 
0°C for 15 minutes, 0.5 ml is filtered through a Nal- 
gene nylon syringe filter, the samples are diluted to 
60% acetone with the addition of 0.25 ml water, and 
0.5 ml is injected into the HPLC. A Hewlett-Packard 
1050 Series HPLC with a Phenomonex Sphericlone 
ODS (2) reverse-phase column (250 mm x 4.5 mm 
with 5 pm particles) is used with a ternary gradient to 
separate and identify the photosynthetic pigments 
(Wright et al, 1991). On cruises that included along- 
track measurements, water from about one meter be- 
low the sea surface was pumped through a hose into a 
bucket in which a Hydrolab Datasonde 3 Multiprobe 
logger was immersed. The Datasonde was used to 
measure temperature, specific conductivity, and salin- 
ity. Data from many of these cruises are made avail- 
able via cruise reports (Culver et al. 1998, Subrama- 


niam et al. 1997a, Subramaniam et al. 1997b, Subra- 
maniam et al. 1997c, Subramaniam et al. 1998), and 
also through a website at the NOAA Coastal Services 
Center (http ://www. esc. noaa. gov/crs/cruises/) . 

6.3 RESEARCH RESULTS 

Our group has conducted eight cruises since the 
beginning of October 1996, including six cruises in 
the South Atlantic Bight (SAB), one in the Gulf of 
Maine, and one in the vicinity of Nantucket Shoals. 
Coverage of the SAB has been extensive due to pro- 
ductive collaboration with the NOAA Southeast Fish- 
eries Science Center. 

The SAB consists of a variety of environments 
including near-coastal and continental shelf regimes, 
the Gulf Stream, and the Sargasso Sea. Smooth tran- 
sitions from low temperature waters with high pig- 
ment concentrations and attenuation coefficients, that 
are typical of nearshore environments, to the higher 
temperature waters with low pigment concentrations 
and attenuation coefficients, that are typical of off- 
shore environments, are often disrupted by Gulf 
Stream eddies that appear on the continental shelf 
The variability in the biological and optical character- 
istics of these regimes complicates temporal or spatial 
analyses of changes in, for example, phytoplankton 
species composition, primaiy production rates, col- 
ored dissolved organic matter (CDOM) concentration, 
suspended sediments, or temperature. 

Coastal waters of the SAB can be classified as 
optical Case II waters. Data from the MAY970B 
cruise (Subramaniam et al,, 1998), on May 5, 1997 
and May 8, 1997, indicate that the waters of Onslow 
Bay and especially Pamlico Sound were sediment- 
dominated. Data from the May cruise suggest that the 
OC2 algorithm (O’Reilly et al., 1998) over-estimated 
chlorophyll and pigment concentrations in these wa- 
ters by a factor of at least 1.5 and up to 15, Chloro- 
phyll a and total suspended solid s concentrations in 
waters of the SAB during the APR98SAB cruise in 
April 1998 were highly variable. Chlorophyll a con- 
centrations tended to be higher nearshore, and total 
suspended solids concentration showed no discernible 
geographic pattern, with high and low concentrations 
being measured at both nearshore and offshore sta- 
tions. High attenuation coefficients for the blue 
wavelengths indicated the presence of high concentra- 
tion of colored dissolved organic material near the 
coasts of South Carolina, Georgia, and Florida how- 
ever, these high coefficients were not observed near 
North Carolina. Overall, the OC2 version 2 algorithm 
significantly over-estimated the chlorophyll a con- 
centration by a factor of 4 and as much as a factor of 
50, however, CZCS pigment concentration was esti- 
mated more successfully, with ratios of measured to 
estimated pigment ranging from 1 to 5. 
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Table 1. Summary of cruises. 


Cruise Location 

Date 


Ship 

Cruise Name 

Southern California Bight 

1/29/96 to 2/18/96 

None 

Jordan 

JAN96CAL 

Mid-Atlantic Bight 

3/1/96 to 3/10/96 


Endeavour 


South Atlantic Bight ! 

4/22/96 to 4/26/96 

None 

Ferrel 

APR96FER 

New York Bight 

5/13/96 to 5/16/96 

None 



Gulf of Maine 

5/18/96 to 5/24/96 

None 


mxern ■ 

Lake Erie 

8/26/96 to 8/30/96 

None 


witi.-i.iti*™ 

Beaufort, North Carolina 

3/13/97 to 3/13/97 

None 



Coastal North Carolina 

mSSESEMMKk 

None 


ixnmnrm 

Gulf of Maine 

6/3/97 to 6/10/97 

None 

isisSBl!®!! 


South Atlantic Bight 

9/6/97 to 9/24/97 

None 


SEP97SAB 

South Atlantic Bight 

11/4/97 to 11/6/97 

None 

R/V Palmetto 

NOV97SAR 

South Atlantic Bight 

4/5/98 to 4/27/98 

NOAA Twin Otter 


APR98SAB 

Nantucket Shoals 

||I|||I|||I 1 1 

None 


JUL98NAN 

South Atlantic Bight 

10/28/98 to 11/25/98 

None 


NOV98SAB 
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Figure 1. OC2 version 2 vs. fluorometric chlorophyll. 
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The results of cruise NOV97SAR (Culver et al, 
1998), conducted from November 3 to 5, 1997, show 
that the OC2 algorithm performs reasonably well at 
stations well off-shore, although, in general the algo- 
rithm tends to over-predict chlorophyll concentrations 
in coastal waters. Also, we found good agreement 
between our in situ optical measurements of normal- 
ized water-leaving radiance and that measured by the 
SeaWiFS sensor. We interpret this to mean that our in 
situ measurements are of good quality, and that our 
goal of validating satellite chlorophyll algorithms us- 
ing these measurements is reasonable and tractable. 

The consistent overestimation of chlorophyll a 
concentration in the SAB by the standard algorithms is 
likely due to the influence of sediments and CDOM 
on ocean color. Preliminary analysis suggests that 
sediment may be the dominant factor near North 
Carolina, and CDOM may be the dominant factor for 
the southern waters. 

When these components do not covary with chlo- 
rophyll a , the algorithms will overestimate chlorophyll 
a concentration. The agreement between our in situ 
optical measurements of normalized water-leaving 
radiance and that measured by the SeaWiFS sensor 
indicate that a chlorophyll algorithm for these waters 
is tractable, however, the coefficients will require ad* 
justment to compensate for the effect of sediment and 
CDOM on the radiance measurements. Overall, our 
results show that it is necessary to establish empirical 
bio-optical relationships for separate coastal regions, 
such as the SAB. 

There is a reasonable relationship between the 
OC2 algorithm-derived chlorophyll concentration and 
the fluorometrically measured chlorophyll concentra- 
tion in the SAB (Figure 1). Small changes in the 
power equation coefficients should yield a better re- 
gional algorithm for the SAB. The relationship be- 


tween HPLC-derived chlorophyll a and the OC2 algo- 
rithm is not as robust and requires further study 

6.4 WORK PLAN 

The project plans to complete three additional 
cruises in the SAB during October 1998 through Sep- 
tember 1999 in order to collect a more complete data 
set using our established sampling protocol. Given 
that analyses thus far have demonstrated that the OC2 
Version 2 algorithm estimates of chlorophyll a con- 
centration for SAB waters exceed acceptable error 
limits, our project will strive to develop a region- 
specific algorithm for SeaWiFS for these waters. 
During the recent NOV98SAB cruise, in addition to 
the deployment of the NOAA CSC PRR600, a Satlan- 
tic free falling profiler (SPMR), the SIMBAD, and an 
above water hyperspectral Spectrex instrument were 
used to obtain water leaving radiance data. This coin- 
cident multiple instrument data set will be used to 
compare the efficacy of each instrument in coastal 
Case II waters. In response to the observation that 
erroneous negative radiances occur in SeaWiFS data 
at 412 nm band, especially in late spring/early sum- 
mer, sunphotometer data will be used to determine the 
suitability of the standard atmospheric correction algo- 
rithms for the SAB. 
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Chapter 7 

Validation of Ocean Color Satellite Data Products 
in Under-Sampled Marine Areas 

Douglas G. Capone and Ajit Subramaniam 
Chesapeake Biological Laboratory 
Solomons, Maryland 

Edward J. Carpenter 
SUNY at Stony Brook, New York 


7.1 INTRODUCTION 

The marine diazotroph, Trichodesmium is a 
planktonic cyanobacterium which occurs throughout 
the tropical and sub-tropical oligotrophic seas. It can 
form large blooms and is thought to be responsible for 
the major fraction of N 2 fixation in the pelagic zone of 
world’s oceans. The total rate of N 2 fixation by this 
organism has been thought to be about 10 Tg annu- 
ally. However, estimates of Trichodesmium' s contri- 
bution to total marine N 2 fixation are relatively crude 
and other independent estimates suggest a higher rate 
of marine N 2 fixation. The annual rate was derived by 
scaling an average trichome N 2 fixation rate with bio- 
mass abundance and distribution data largely taken 
from historical plankton surveys. However, there are 
unique difficulties in quantifying the biomass of a 
buoyant, colonial alga such as Trichodesmium and the 
quantitative accuracy of many of these surveys are 
suspect. One of the most appealing prospects of using 
satellites in biological oceanography has been the po- 
tential to derive flux estimates from inventories of 
standing stock calculated from changes in optical 
properties of the ocean. Using algorithms to estimate 
global distributions of Trichodesmium biomass, we 
can calculate nitrogen fixed by this organism, and 
hence model its contribution to new production in the 
world oceans. Towards this goal, we proposed to use 
a previously developed hyperspectral optical model 
for the Trichodesmium to construct sensor specific 
biomass algorithms. 

The goal of our proposed work was to develop 
sensor specific algorithms to quantify the biomass of 
Trichodesmium in the world oceans. The specific ob- 
jectives of our proposal were to: 

• Use a hyperspectral optical model for Trichodes- 
mium to derive sensor specific Trichodesmium 

chlorophyll algorithms. 


• Constrain and validate these algorithms using 
field measurements. Also use the field measure- 
ments to validate standard chlorophyll algorithms 

• Estimate the error in using standard algorithms to 
estimate chlorophyll in a Trichodesmium bloom 
Develop flags for switching algorithms if neces- 
sary. 

We were required by our contract to acquire and 
process data from at least one cruise per year. The 
data included measurements of vertical profiles of 
spectral downwelling irradiance, upwelling radiance, 
spectral attenuation and remote-sensing reflectance 
calculated from the profile data, photosynthetically 
active radiation, primary productivity, and concentra- 
tions of chlorophyll a and other pigments We were 
required to purchase and make available to the instru- 
ment pool, an AC-9 instrument and a turnkey free-fall 
profiling radiometer system. Here we outline our pro- 
gress during the first year of our SIMBIOS project to 
develop an algorithm for detection and quantification 
of Trichodesmium biomass for the SeaWiFS sensor 
We also describe our field program for validating the 
OC2 algorithm and obtaining in situ bio-optical data 
in under-sampled marine areas. 

7.2 RESEARCH RESULTS 

The work detailed below is of an effort in prog- 
ress and intended to provide an interim update to the 
SIMBIOS project. Our initial approach was to derive 
a Trichodesmium chlorophyll specific algorithm from 
an optical model, and then constrain this algorithm 
based on field data. Our optical model detailed in 
Subramaniam et al. (1999) showed that at low con- 
centrations, Trichodesmium does not have singular 
enough optical properties to be uniquely identified 
(Figure 1). But as the concentration of Trichodesmium 
in the water column increases, its optical properties 
dominate and it is possible to identify and quantify 
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this organism (Figure 1). Therefore, we initially de- 
cided to develop a branching algorithm, where at low 
Trichodesmium concentrations the algorithm produced 
total chlorophyll (i.e. not Trichodesmium specific), but 
at higher Trichodesmium concentrations the algorithm 
resulted in Trichodesmium specific chlorophyll. We 
formulated the algorithm as (R555-R490/R510) where 
R555, R490, R510 are the SeaWiFS remote sensing re- 
flectance corresponding to bands 3, 4 and 5 respec- 
tively. We parameterized this formulation by fitting a 
power equation to the results of the previously devel- 
oped remote-sensing reflectance based hyperspectral 
optical model for Trichodesmium. However, we were 
not satisfied by the results of this formulation and 
have been working on other approaches and band 
combinations. 

Although the branching algorithm approach is 
very attractive for producing global maps of Tricho- 
desmium abundance based on a sound theoretical 
model, it is essential to characterize the chlorophyll 
concentration threshold where the algorithm results in 
Trichodesmium specific chlorophyll concentration. 
The optical model suggested that this threshold was 
about 0.5 mg/m 3 Trichodesmium Chi a. While we did 
not encounter such high concentrations of Trichodes- 
mium during the Trichonesia cruise, preliminary 
analysis of field data from a cruise in the South Atlan- 
tic Bight shows that this threshold is probably much 
higher - around 3 mg/m 3 . We need more field data 
with Trichodesmium cell counts to characterize this 
threshold. In the meantime we have developed a flag 
to study global occurrence of Trichodesmium blooms. 
This flag is based on the combination of high reflec- 
tance at 555 nm, and absorption at 490 nm. We are 
comparing this to climatological data on Trichodes- 
mium occurrence - in regions and seasons where 
Trichodesmium is known to occur historically. We 
think we are very close to completing an initial 
Trichodesmium chlorophyll algorithm and will com- 
municate it to the SIMBIOS office as soon as we are 
satisfied with it. 

Our 1998 field campaign was a 35-day research 
cruise in the South Western Pacific ocean. We left 
Lyttleton, New Zealand on the 24 March 1998 and 
arrived in Suva, Fiji on the 30 April 1998. We had 
requested cruise time in January-February when the 
organism is most abundant in this region. However, 
due to ship availability, our cruise was pushed to 
April. This, in combination with the effects of the EI- 
Nino, resulted in our encountering lower Trichodes- 
mium populations than we had anticipated. After our 
port call in New Caledonia on 30 March, we had to 
head north to avoid Cyclone Zumann. Thus we were 
under overcast skies during the portion of the cruise 
with some of the highest chlorophyll concentrations. 


In total, we occupied 5 1 stations and made 33 optics 
deployments including 2 diel deployments for Inher- 
ent Optical Property measurements. The in situ optical 
data, chlorophyll concentrations, and carbon and 
nitrogen fixation rates have been submitted to 
SeaBASS. 

We acquired a Wetlabs AC-9 in January 1998. 
This instrument was field tested during a SIMBIOS 
cruise (CALCOFI 9802) and recalibrated after that 
cruise. This instrument was used on the Trichonesia 
cruise and is now available at CBL. 

A Satlantic SeaWiFS Profiling Multispectral Ra- 
diometer (SPMR) with a matched SeaWiFS Mul- 
tispectral Surface Reference (SMSR) was procured in 
March 1998. The SPMR is a 13 channel free falling 
profiler that measures in situ downwelling irradiance 
(Ed') and upwelling radiance (Lu") at 340, 380, 412, 
443, 490, 510, 520, 555, 565, 619, 665, 670, and 683 
nm. The SMSR is a floating reference that measures 
downwelling irradiance (Ed 4 ) just above water and 
upwelling radiance (Lu‘) 70 cm below the surface. 
The SPMR includes a 300-m power telemetry cable, 2 
axis tilt sensors, and external temperature and con- 
ductivity sensors. The SMSR includes a 100-m power 
telemetry cable and 2 axis tilt sensors. The package 
includes a deckbox with computer interface and power 
supply, data logging and display software, and a Hi- 
tachi pentium laptop computer with an extra PCMCIA 
serial card and cables. 

The SPMR was used on the Trichonesia cruise in 
March-May, on NOAA cruises in July and November, 
and a CALCOFI cruise in September. Initially there 
were problems with gain switching and some data 
from the Trichonesia cruise and JUL9SNAN cruise 
were compromised because of this. This has been 
subsequently fixed by Satlantic. The instrument has 
been recalibrated twice at Satlantic and twice at UC 
Santa Barbara since March 1998. It has been shipped 
to the Indian Ocean for use in the Indian Ocean Ex- 
periment (INDOEX). 

7.3 WORK PLAN 

We will continue to work on refining the Tricho- 
desmium algorithm. Specifically we will analyze the 
field data collected during the Trichonesia (1998) and 
INDOEX (1999) cruises to define the Trichodesmium 
specific chlorophyll concentration that can be detected 
by the algorithm and test the sensitivity of the algo- 
rithm. We will work on algorithms specific to the 
OCTS, the MOS, and the MODIS sensors. We will 
participate in the INDOEX cruise in March 1999 and 
a cruise off northern Australia in October 1999. 
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Figure 1. Optical model of Trichodesium remote sensing reflectance for a) 0.1 mg Chl/m 3 , b) 1.0 mg 

Chi/ m 3 , and c) 10.0 mg Chi/ m 3 . 
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Chapter 8 

Stray Light and Atmospheric Adjacency Effects for 
Large-FOV, Ocean-Viewing Space Sensors 

Kendall L. Carder 
University of South Florida, 

St. Petersburg, Florida 


8.1 INTRODUCTION 

The radiance of dark targets (e g. oceans and 
lakes) when measured by large-field imagers such as 
MODIS and SeaWiFS can be enhanced by stray light 
within the sensor and by radiance atmospherically 
scattered from bright, adjacent targets. One of the 
goals of this project is evaluation of this stray light 
and estimation of its effect on derived products of 
remote color imagery. 

8.2 RESEARCH ACTIVITIES 

The first year of this project focused on collecting 
field data to evaluate the performance of SeaWiFS. 
To assure field data quality, a significant effort was 
expended in characterizing and calibrating field in- 
struments as well as standardizing field methodology. 
Field sites were selected where adjacency artifacts 
would be likely. One site is a deep ocean canyon sur- 
rounded by bright carbonate banks, the Tongue of the 
Ocean (TOTO); another is Lake Okeechobee, Florida 
where dark water is surrounded by bright vegetation. 
With funding from other sources, additional field ex- 
periments were applied to this study and help with 
evaluation of SeaWiFS performance. These included 
sampling on Tampa Bay, the West Florida Shelf, 
Exuma Sound, the San Juan Islands in Puget Sound, 
and Lake Tahoe. 

8.3 RESEARCH RESULTS 

The degradation of SeaWiFS infrared bands, 
which are used to evaluate the type and amount of 
aerosol path radiance, caused a perturbation in the 
type of aerosol correction to be made by the SeaWiFS 
atmospheric correction scheme. This resulted in over 
estimation of the blue-richness of the aerosol scatter- 
ing. Furthermore, high cirrus clouds have a tendency 
to exaggerate this effect, due to enhanced aerosol ab- 
sorption in the oxygen band near 762 nm. To obviate 
the worst effects of both problems, we developed an 


aerosol-typing algorithm using bands 6 and 8 over 
clear waters. This approach forced the calibration of 
band 6 relative to band 8 to provide normalized water- 
leaving radiance values at band 6 consistent with the 
Gordon and Clark (1981) for scenes with known aero- 
sol type. For scenes of unknown aerosols, the new 6-8 
band combination delivered appropriate aerosol types 
with SEADAS. This prevented blue exaggeration of 
the aerosol type and prevented excessive removal of 
blue radiance from the scene, which otherwise would 
have resulted in blue-poor water-leaving radiance val- 
ues. This approach provided atmospherically corrected 
scenes that did not became negative at blue wave- 
lengths for gelbstoff-rich (blue-absorbing) areas This 
approach also provided clear water, normalized water- 
leaving radiance values at 555 nm consistent with 
Gordon and Clark (1981) using the calibration factors 
distributed by the SeaWiFS Project on 6 January 
1998. Furthermore, it provided a data distribution for 
clear waters using ratios of bands 1 and 2 versus bands 
2 and 5 that are consistent with those found by Carder 
et al. (1999) using the SeaBASS in situ data set. 

Using the above calibration scheme, SeaWiFS 
images of the TOTO, the West Florida Shelf, and re- 
gions between were compared to field data using the 
default NASA chlorophyll algorithm and the Carder et 
al. (1999), MODIS-Iike, semianalytic chlorophyll al- 
gorithm. For these spring scenes, the pigments were 
somewhat packaged with little in the way of photo- 
protective pigments. Both algorithms performed 
similarly (within 5%) for Case 1 waters. When 
gelbstoff absorption overwhelmed the pigment ab- 
sorption, however, the semianalytic algorithm deliv- 
ered chlorophyll values within 25% of measured val- 
ues; while die default NASA algorithm overestimated 
chlorophyll by as much as a factor of 4. These scenes 
will be re-evaluated^ with the new calibration- 
degradation curves of SEADAS 3.2 during the next 
contract year, but the new results are not expected to 
deviate significandy from these results using the 6-8 
band combination except for turbid waters. 

In terms of evaluating atmospheric adjacency and 
stray light effects due to infrared radiance emanating 
from the bright foliage surrounding Lake Okeechobee 
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or blue-green light reflected from the bright shallow 
banks surrounding TOTO, our preliminary findings 
suggest the following: 

• Lake Okeechobee has so much horizontal vari- 
ability in the water-leaving radiance field at infra- 
red wavelengths that it is at times an impractical 
site at which to evaluate stray infrared light from 
its surrounds. It does offer, however, so much 
gelbstoff absorption at 412 nm that we are con- 
sidering performing atmospheric corrections 
based upon the blue end of the spectrum. 

• Data we collected at the Friday Harbor Labora- 
tory in the San Juan Islands suggest that this site 
provides a viable alternative to Lake Okeechobee 
for evaluation of stray infrared light and/or at- 
mospheric adjacency effects. Hyperspectral 
Imager for Low Light Spectroscopy (PHILLS) 
aircraft data Curtis Davis (NRL) suggests that one 
can find a red-edge effect (exaggerated infrared 
light) in the water near sunlit tree lines in data 
collected at 10,000 feet and with visibility in ex- 
cess of 50 km. These data will be evaluated next 
year along with SeaWiFS data collected at high 
resolution over Griffin Bay (San Juan Island) 
which is surrounded by forest land. 

• Chlorophyll values farther than 2 pixels offshore 
and east from bright, shallow banks of the TOTO 
were not significantly different than those found 
farther offshore. Extreme caution must be taken, 
however, to avoid misinterpreting conditions 
where gelbstoff- and particle-laden waters from 
the shallows are transported offshore. Next year 
an evaluation of the electronic over-shoot, and/or 
hysteresis on the west sides of bright targets will 
also be evaluated on the eastern edge of the 
TOTO. More accurate calibration of the 
SeaWiFS sensor will help in this effort. 


• A recently published paper by Reinersman et al. 
(1998) on calibration of aircraft and small foot- 
print, spacecraft sensors using cloud shadows was 
dependent upon the development of a Monte 
Carlo model of radiance from single spherical 
clouds. This model is being modified to permit 
evaluation of atmospheric adjacency effects of 
clouds on radiance values received by a space 
sensor when viewing the ocean through an open- 
ing (elliptical cylinder) in a cloud bank Algo- 
rithms that are more immune to the adjacency ef- 
fects of clouds will be sought. Much of the light 
scattered into the sensor from clouds will be from 
molecules and stratospheric aerosols and thin cir- 
rus clouds. A strategy using the effect of the oxy- 
gen absorption line on band 7 versus results from 
bands 6 and 8 may provide an indication of the 
effective height of any aerosol scattering for 
SeaWiFS, whereas MODIS has a cirrus discrimi- 
nator built into its data processing scheme 

8.4 WORK PLAN 

The main goal of the project will remain the 
same, to evaluate stray light effects on derived prod- 
ucts from ocean color imagery. Since MODIS and 
OCTS will not be available, the effort will be concen- 
trated on SeaWiFS, Advanced Visible and Infrared 
Imaging Spectrometer (AVIRIS), and PHILLS. 

The TOTO site will be visited in mid- April It is 
planned to use PHILLS overflights of Griffin Bay, 
Tampa Bay, and the West Florida Shelf. The aircraft 
data will be used to refine techniques to detect atmos- 
pheric adjacency effects and bright infrared reflec- 
tance due to vegetation. 
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Chapter 9 

Bio-Optical Measurements at Ocean Boundaries 
in Support of SIMBIOS 

Francisco P. Chavez 

Monterey Bay Aquarium Research Institute 
Moss Landing, California 


9.1 INTRODUCTION 

This project consists of a mooring program and 
supporting cruise-based measurements aimed at quan- 
tifying the spectrum of biological and chemical vari- 
ability in the equatorial Pacific. The mooring program 
was designed to obtain continuous time series of bio- 
logical and chemical properties on a time scale that is 
equivalent to measurements of currents, local winds 
and temperature structure.The project has the follow- 
ing general objectives: 

• to understand the relationships between physical 
forcing, primary production, nutrient supply and 
the exchange of carbon dioxide between ocean 
and atmosphere in the equatorial Pacific; 

• to describe the biological and chemical responses 
to climate and ocean variability; 

• to describe the spatial, seasonal and interannual 
variability in near surface plant pigments, primary 
production, carbon dioxide and nutrient distribu- 
tions, and 

• to obtain near real-time bio-optical measurements 
to ground-truth satellite measurements of ocean 
color. 

9.2 RESEARCH ACTIVITIES 

Moorings 

Bio-optical and chemical sensors were deployed 
on two moorings of the TAO array in the equatorial 
Pacific at 0°, 155°W and 2°S, 170°W. The following 
sensors formed the core set: 

• A Biospherical PRR-620 located 3m above the 
water surface, to measure downwelling irradiance 
at 412, 443, 490, 5 10, 555, 656nm and PAR. 

• A ApCCh system to measure the difference in 
partial pressure of CO2 between ocean and at- 
mosphere (supported by NOAA). 

• Two Satlantic OCR- 100s located approximately 
1.5m below the surface, to measure upwelling ra- 


diance at 412, 443, 490, 510, 555, 670 and 
683nm. 

• Two Wetlabs miniature fluorometers at approxi- 
mately 1.5m and 20m below the surface to meas- 
ure stimulated in vivo fluorescence. 

• Two Biospherical MCP-200s located at 10m and 
30m below the surface to measure downwelling 
irradiance at 490nm. 

• A Biospherical PRR-600 located 20m below the 
surface to measure downwelling irradiance at 
412, 443, 490, 510, 555, 670nm and PAR, plus 
upwelling radiance at 412, 443, 490, 510, 555, 
670 and 683nm. 

Daily noon-time bio-optical data were delivered 
via service ARGOS in near real time to MBARI, and 
then via automated FTP to the SeaBASS database, 
after some processing and quality control Higher fre- 
quency, publication-quality data (10 or 15 -minute 
intervals) were recovered at approximately six month 
intervals, and sent to the SeaBASS database after 
processing and quality control. Some derived prod- 
ucts, such as OC2V2 chlorophyll from surface up- 
welling radiances, and mean chlorophyll over the up- 
per 20m of the water column (Morel, 1988) were in- 
cluded in these data files. Using simple algorithms 
derived from optical profile data the upwelling radi- 
ances measured by the moored radiometers have been 
converted to normalized water-leaving radiances 
(Lwn), for the SeaWiFS calibration/vahdation effort. 

In addition to the two mooring installations de- 
scribed above, MBARI has also maintained smaller 
bio-optical packages at four sites in the equatorial Pa- 
cific: 2°N, 180°; 2°S, 140°W; 2°N, 140°W and 2°N, 
110°W. TTiese bio-optical packages produced meas- 
urements of downwelling incident irradiance at 
490nm, upwelling radiance from ~1.5m depth at 412, 
443, 490, 510, 555, 670 and 683nm, plus stimulated in 
vivo fluorescence. Data, including derived (see 
above) were delivered in near real time to MBARI and 
the SeaBASS database. A bug in the original software 
precluded collection of significant data prior to March 
of 1998. 
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In situ measurements 

All cruises were undertaken aboard the NOAA 
ship Ka 'imimoana, with the exception of GP6-98-RB 
aboard the NOAA ship Ronald II Brown (Table 1). 
Essentially, the cruise-based measurements consisted 
of chlorophyll (using the fluorometric method de- 
scribed by Chavez et al., 1995) and nutrient profiles (8 
depths, 0-200m) obtained at CTD stations between 
8°N and 8°S across the Pacific from 95°W to 165°E. 
On selected cruises (primarily the 155°W and 170°W 
meridional transects), primary productivity measure- 
ments and optical profiles, using the Satlantic Profil- 
ing Multispectral Radiometer, were also performed. 
These data were archived at MBARI and the chloro- 
phyll and optical profile measurements provided to 
NASA’s SeaBASS database. 

9.3 RESEARCH RESULTS 

Chavez et al. (1998) describe the biological- 
physical coupling observed in the central equatorial 
Pacific during the onset of the 1997-98 El Nifio. Fig- 
ure 1 shows the time series of chlorophyll, SST and 
ApC0 2 from the mooring located at 0°, 155°W for the 
time period covering the peak of the 1997-98 El Nifio 
and the subsequent recovery of the phytoplankton 
community. Our data from the four smaller bio-optical 
packages describe the physical and biological effects 
of the passage of Tropical Instability Waves. While 
the corresponding SeaWiFS data in large part agree 
closely with the mooring-derived chlorophyll concen- 
trations, there are time periods where this is not the 


case. This may in part be due to the spatio-temporal 
averaging that is inherent in the SeaWiFS 9km/8-Day 
data. We have performed SeaWiFS chlorophyll com- 
parisons with extracted chlorophyll concentrations for 
all of the 1998 equatorial Pacific cruises, with good 
agreement between the two data types. 

9.4 WORK PLAN 

During 1999 data collection will continue as for 
1998. The two major mooring installations at 0°, 
155°W and 2°S, 170°W are now operating well after 
some initial problems with damage and loss of instru- 
ments. The smaller bio-optical packages at four addi- 
tional locations across the Pacific are also performing 
well and data collection at those sites, as described 
above, will also continue. The program of cruise- 
based measurements of chlorophyll, nutrients, primary 
productivity and bio-optical profiles will continue on 
up to eight equatorial Pacific cruises during 1999, 
with scheduled SeaWiFS LAC where possible, to en- 
hance the probability of obtaining valid matchups. 

The major data processing objective for 1999 is to 
automate the near real-time data provision to the Sea- 
BASS database as much as possible to aid in the cali- 
bration/validation effort. In addition we aim to im- 
prove our data processing routines for the high resolu- 
tion data that are downloaded from the 0°, 155°W and 
2°S, 170°W moorings at approximately 6 month in- 
tervals. This should decrease the lag time between 
data download and provision of products to the Sea- 
BASS database. 
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Table 1. Summary of cruises during which in situ data have been obtained by MBARI in support of SIM- 
BIOS. All cruises were undertaken aboard the NOAA ship Ka'imimoana, with the exception of GP6-98-RB 
aboard the NOAA ship Ronald H. Brown. Meridional transects indicate the lines occupied by the ship. Along 
each line, CTD stations were performed approximately every degree of latitude from 8°N to 8°S. Measurements 
consisted of extracted chlorophyll (Chi) plus nitrate, phosphate and silicate (Nutrients) at 8 depths between 0 and 
200m. On selected cruises, primary productivity (PP) measurements were also made using 'X incubation tech- 
niques, and daily optical profiles with the Satlaintic Profiling Multispectral Radiometer (SPMR) were obtained. 
During GP7-98-KA, a SIMBAD radiometer was also used, on loan from the SIMBIOS instrument pool. 


Cruise ID 

Dates 

Meridional transects 

Measurements 

GP6-97-KA 

27-Sep-97 to 30-Oct-97 

125°W and 140°W 

Chi, PP, Nutrients, SPMR 

GP7-97-KA 

06-Nov-97 to 17-Dec-97 

155°W, 170°W and 180° 

Chl, PP, Nutrients, SPMR 

GP1-98-KA 

05-Feb-98 to 13-Mar-98 

95 °W and 110°W 

Chi, Nutrients 

GP2-98-KA 

18- Apr-98 to 20-May-98 

125°W and 140°W 

Chl, Nutrients 

GP3-98-KA 

02-Jun-98 to 03-Jul-98 

155°W and 170°W 

Chl, PP, Nutrients, SPMR 

GP4-98-KA 

07-M-98 to 03-Aug-98 

165°E and 180° 

Chl, Nutrients 

GP5-98-KA 

05-Sep-98 to 09-Oct-98 

125°W and 140°W 

Chl, Nutrients 

GP6-98-RB 

11 -Oct-98 to 13 -Nov-98 

95°W and 110°W 

Chl, Nutrients 

GP7-98-KA 

19 -Oct-98 to 13 -Nov-98 

155°W and 170°W 

Chl, PP, Nutrients, SPMR 

GP8-98-KA 

18-Nov-98 to 12-Dec-98 

180° and 170°W 

Chl, Nutrients 

GP1-99-KA 

22-Jan-99 to 24-Feb-99 

125°W and 140°W 

Chl, Nutrients 
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Figure 1. Time series of Sea Surface Temperature (SST - heavy shaded line), chlorophyll a averaged over 
the upper 20m of the water column from moored spectroradiometer data (heavy solid line) SeaWiFS OC2V2 
chlorophyll (closed triangles) and ApC0 2 (thin solid line). 
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Chapter 10 

Remote Sensing of Ocean Color in the Arctic: 
Algorithm Development and Comparative Validation 


Glenn F. Cota 
Old Dominion University 
Norfolk , Virginia 


10.1 INTRODUCTION 

The Optical Research Consortium of the Arctic 
(ORCA) is a collaborative project between G.F. Cota 
of Old Dominion University (ODU), T. Platt and 
W.G. Harrison of the Bedford Institute of Oceanogra- 
phy (BIO), S. Sathyendranath of Dalhousie University 
and S. Saitoh of Hokkaido University. ORCA initi- 
ated bio-optical data collections in August, 1994 and 
has conducted 10 high latitude cruises in the Bering 
Sea (7/95 = Ber95, 7/96 = Ber96), the Labrador Sea 
(10-1 1/96 = Lab96, 5-6/97 - Lab97, 6-7/98 = Lab 98), 
the Gulf of Alaska (10/97 = Goa97), and in the 
Northwest Passage near Resolute Bay, NWT, Canada 
(8/94 = Res94, 8/95 = Res95, 8/96 = Res96, 8/98 = 
Res98). ORCA has now collected almost 500 in- 
water optical profiles and most of them have been 
submitted to NASA’s SeaBASS database. 

The ORCA program depends upon coordination 
and cooperation between international co-investigators 
with support from NASA, NASDA and Canadian 
sources. ODU has had primary responsibility for bio- 
optical and biogeochemical data collection with as- 
sistance from all collaborators. ODU has been re- 
sponsible for data processing, distribution and submis- 
sion to SeaBASS. Data analyses are being conducted 
by all co-investigators. Modeling is being done pri- 
marily by BIO and Dalhousie. Our algorithm devel- 
opment efforts encompass most of the potential data 
products listed for SeaWiFS, MODIS and GLI, and 
addresses product validation for a specific area of 
concern, i.e. high latitudes and coastal waters. 

Observational suites depend primarily on field 
program duration, the sampling time, berths and per- 
sonnel available. Longer, dedicated field programs 
(i.e. Resolute) include more comprehensive observa- 
tions. Core observations always include in-water 
spectrorad.ometer profiles for apparent optical prop- 
erties (AOPs) and discrete profiles of fluorometric 
chlorophyll. Absorption spectra for particulate, ex- 
tracted particulates, and soluble materials are also 
available from most cruises. The observational suites 
for Resolute and the Labrador Sea are typically more 
comprehensive with in situ inherent optical property 


(IOPs) profiles (AC -9, Hydroscat -6: Resolute only), 
above-water spectral remote sensing reflectance 
R^X), spectrophotometric and/or HPLC pigments, 
primary production and/or photosynthesis versus irra- 
diance, inorganic nutrients, 15N-nutrient uptake, par- 
ticulate organic carbon and nitrogen. 

10.2 RESEARCH ACTIVITIES 

From May, 1997 through August, 1998 four 
SIMBIOS validation cruises for OCTS or SeaWiFS 
were conducted with one on an early-start supple- 
mental contract (Labrador Sea 5-6/97) and three under 
the main contract (Gulf of Alaska 10/97, Labrador Sea 
6-7/98, and Resolute Bay 8/98). Our scheduled 
August, 1997 Resolute field program was cancelled 
because OCTS had failed and SeaWiFS was not yet 
operational. 

Bio-optical, biogeochemical and experimental ob- 
servations were made as frequently as logistical and 
weather considerations permitted. Data collections 
and processing for biogeochemical and optical obser- 
vations conform closely to Joint Global Ocean Flux 
Studies (JGOFS, 1991) and SeaWiFS protocols 
(Mueller and Austin 1992, 1995). 

Water column profiling consisted of a bottle or 
pump cast for discrete samples preceded or followed 
by optical casts. Discrete water samples were col- 
lected from standard depths and depths corresponding 
to features of particular interest such as fluorescence, 
beam attenuation and/or density maxima. Biogeo- 
chemical observations included fluorometric chloro- 
phyll and phaeopigments, spectrophotometric (Reso- 
lute) or HPLC 0-abrador Sea) pigments, absorption 
spectra of particulates, methanol -extracted particu- 
lates, and dissolved colored organic material 
(CDOM), inorganic nutrients, particulate organic car- 
bon and nitrogen, and selected other observations such 
as total suspended matter and cell count samples. 
Experimental observations at Resolute and the Labra- 
dor Sea included 14 C and 13 C primary production and 
15 N-nutrient uptake of nitrate and ammonium 

Optical stations include hand-held sun photome- 
try and photographs of sea and sky conditions. In- 
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water optical observations were made with a Satlantic 
SeaWiFS Profiling Multichannel Radiometer (SPMR) 
and SeaWiFS Multichannel Surface Reference 
(SMSR). Downwelling irradiance E d (0‘,X) and up- 
welling radiance Lu(0\X.) were measured just below 
the sea surface at 13 wavebands with the SMSR sen- 
sors mounted on a floating, tethered buoy with an um- 
bilical cable for power and data transmission. The 
SPMR also has 13 channels for vertical profiles of 
downwelling irradiance E d (z, X) and upweiling radi- 
ance L u (z, X), and was deployed in a free-fall mode 
with a kevlar cable for power and data transmission. 
The SPMR includes sensors for tilt and pitch, pres- 
sure, conductivity and temperature (Ocean Sensors), 
and fluorescence (WETLabs WetStar). The SMSR 
and SPMR measurements were always made concur- 
rently and the sensors were deployed 20-100 m from 
the vessel to avoid ship-shadow effects. Above-water 
reflectance measurements of the sea surface were 
made with an Analytical Spectral Devices (ASD) Per- 
sonal Spectrometer II (350-1050 nm, 1.4 nm band- 
pass,). Incident spectral irradiance E d (0 + ,X) was esti- 
mated with a gray Spectralon card (10% reflectance) 
and water-leaving radiance L u (0 + ,X) was corrected for 
reflected sky light (K. Carder and C. Davis, pers. 
Comm.). Spectrometers have been calibrated 1-2 
times per year depending upon the number of field 
campaigns. Profiles of IOPs were measured with 1-2 
WET Labs AC-9 absorption and beam attenuation 
meters with 9 channels and 25 cm path lengths for 
total and/or soluble constituents <0.2um. A HOBI 
Labs Hydroscat-6 was used to measure backscatter in 
6 spectral bands. The latter three types of measure- 
ments were only made at Resolute. 

AOPs, fluorometric chlorophyll and absorption 
spectra data from all cruises have been submitted to 
SeaBASS, except the Lab98 data set. The downwel- 
ling irradiance E d sensor malfunctioned during Lab98, 
and the data processing will require modifications of 
Satlantic software or estimation of normalized water- 
leaving radiance L W1 from models. Additional ancil- 
lary data are also being submitted as they become 
available from our lab and collaborators. 

10.3 RESEARCH RESULTS 

High northern latitude waters have unique bio- 
optical properties compared with temperate data and 
models (Sathyendranath et al. 1999b). These waters 


have a large dynamic range of biomass, which must be 
considered in algorithm development. 

Results to date reveal that diatoms usually domi- 
nate phytoplankton assemblages, but prymenesiophyte 
blooms may alter bio-optical signatures (Sathyen- 
dranath et al. 1999a). The relatively large diatom cells 
are highly packed with high chlorophyll-specific at- 
tenuation in the blue. Regional differences have been 
found in soluble attenuation with the highest values in 
the Canadian Arctic Archipelago. Blue-green band 
ratio algorithms provide reasonable chlorophyll and 
net daily primary production retrievals when properly 
timed for high latitudes. However, the current 
SeaWiFS OC2 chlorophyll algorithm underestimates 
biomass at moderate to high concentrations (Figure 1). 
Solar-induced fluorescence should provide reasonable 
instantaneous estimates of primary production with 
robust backscatter corrections. The universality of 
any algorithms requires extensive evaluation and vali- 
dation. 

10.4 WORK PLANNED 

Cruise plans for 1999 are not yet finalized, but 
should include a spring (May-June) Labrador Sea 
cruise and one more to be determined. With 
NOAA/NDBC, a time-series of bio-optical observa- 
tions from the Chesapeake Light Tower is being initi- 
ated. As proposed ORCA is conducting 1-3 cruises 
per year, and contributing critical core observational 
data as well as a variety of ancillary bio-optical and 
biogeochemical data. 

Core observations have been delivered within 90 
days and include profiles of AOPs made with our 
Satlantic profiling spectroradiometer, fluorometric 
chlorophyll, and normally particulate and dissolved 
absorption spectra. Additional observations for which 
ODU has primary responsibility for are being deliv- 
ered within 12 months, while contributions from in- 
ternational collaborators are submitted as they become 
available. 

Our goal is to make raw data and derived prod- 
ucts generally available via the web but our site is still 
under construction. Efforts in 1999 will concentrate on 
data analyses, modeling, and publication of results. 
Two manuscripts have been submitted and several 
more are in preparation. 
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Figure 1. The current best-fit line of our high latitude data (minus Lab98) versus the current SeaWiFS OC2 
chlorophyll algorithm. 


49 



SIMBIOS Project Annual Report 


Chapter 11 

High Frequency, Long-Time Series Measurements from 
the Bermuda Testbed Mooring in Support of SIMBIOS 

Tommy Dickey, S. Zedler and D. Sigurdson 
University of California at Santa Barbara , 

Santa Barbara , California 


11.1 INTRODUCTION 

The primary goal of the present activity is to pro- 
vide high frequency, long-term time series bio-optical 
data from the Bermuda Testbed Mooring (BTM) in 
support of SIMBIOS (Dickey et al., 1998). The site is 
located about 80km southeast of Bermuda at -31° 43’ 
N, 64° 10* W in waters of ~4530m depth. 

The BTM study, which is funded by NASA, NSF, 
ONR, and NOPP, supports a variety of scientific re- 
search including SIMBIOS calibration and validation 
activities. Our study contributes to the development 
of methodologies and capabilities for synthesizing 
data derived from several satellite ocean color mis- 
sions (Esaias et al., 1995), particularly SeaWiFS. Key 
NASA supported BTM optical measurements (every 
15min during daylight) presently include: surface 
downwelling spectral irradiance (7 X's, SeaWiFS 
matched) and subsurface (nominal depths of 14 and 
32m) downwelling spectral irradiance and upwelling 
spectral radiance (7 Vs for each, SeaWiFS matched). 

Products include: radiance and irradiance ratios, 
spectral attenuation coefficients, and upwelling spec- 
tral radiances just beneath the surface (0-). These data 
provide necessary links to and interpolation of radio- 
metric data for remotely sensed observations of ocean 
color (suite of satellite color imagers) which can be 
used to estimate biomass and primary productivity 
globally. 

Our newly developed telemetry system is now 
providing near real-time radiometric data to our lab. 
We are in the process of developing system software 
which will make key data available to the SIMBIOS 
project on a daily basis. 

Highly complementary, value-added optical, me- 
teorological, and physical data products are provided 
through our complementary NSF BTM study (e.g., 
Dickey et al., 1998). Our comprehensive measure- 
ments can be used to relate bio-optical signals to in- 
herent optical properties and to biogeochemical meas- 
urements and to test interdisciplinary models of phy- 
toplankton biomass and productivity and carbon flux. 


11.2 RESEARCH ACTIVITIES 

A major objective of the SIMBIOS program is to 
estimate normalized spectral water-leaving radiance, 
Lwn, by using in situ measurements. Our approach is 
to utilize moored instruments to obtain high temporal 
resolution measurements of fundamental radiometric 
variables to estimate upwelled spectral radiance just 
beneath the ocean surface, Lu(O-), which can be used 
to determine Lwn. 

The mooring methodology, which allows for col- 
lection of very high temporal resolution, long-term 
data, has several advantages including the ability to 
obtain high numbers of M match-up" data for compari- 
sons with SeaWiFS (e g., Mueller and Austin, 1992). 
However, the use of moorings for obtaining optical 
data is a relatively new approach and has been used by 
very few groups (e.g., Dickey. 1991; Smith et al., 
1991; Clark et al., 1997; Dickey et al., 1998). 

Moored optical observations are more challenging 
than their ship-based counterparts for several reasons 
(e.g., less frequent calibration, biofouling, etc ). Tlius, 
it is desirable to intercompare fundamental radiomet- 
ric data and derived products based on nearly concur- 
rent, co-Iocated data from moored and ship-based in- 
strumentation. We have been collaborating with Dr. 
David Siegel (UCSB), who has been collecting com- 
parable vertical profile data during approximately 
monthly (every other week in springtime) shipboard 
cruises near the BTM site. 

Our BTM data products are being intercompared 
with ship-based data collected as part of Dr Siegel's 
Bermuda Bio-optical Program (BBOP) program (Sie- 
gel et al., 1995a, b; Siegel and Michaels, 1996). It 
should be noted that the BTM and BBOP measure- 
ments are displaced in space and to a lesser extent in 
time. In addition, there are several other significant 
differences. Thus, correlations between the BTM and 
BBOP data sets should not necessarily be expected to 
be great. 

The BTM Deployment 7, 8, and 9 (May 3, 1997 - 
March 31, 1998) data sets encompass periods when 1) 
only the Japanese ocean color imager (OCTS) was 
collecting data, 2) neither OCTS nor SeaWiFS was 
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collecting data, and 3) when only SeaWiFS was col- 
lecting data. The time series of meteorological, physi- 
cal, and optical data are described in Dickey et al. 
(1998). Time series measurements of surface and near 
surface optical properties were made using moored 
optical radiometer systems. Measurements include 
surface spectral downwelling irradiance (Ed (0+, X) at 
X = 412, 443, 490, 510, 555, 665, and 683nm (lOnm 
bandwidth at half power) at 6Hz for 20sec every 
15min during daylight and at midnight. Subsurface 
(nominal depths of 14 and 32m) downwelling irradi- 
ance (Ed (z, X) and upwelling radiance, Lu (z, X) data 
at X = 412, 443, 490, 510, 555, 665, and 683nm are 
collected following the same regimen. 

Our work is demonstrating that moored radiomet- 
ric instrumentation can be used to provide a large 
number of in situ match-up data for comparisons with 
SeaWiFS observations in near real-time. The present 
approach uses a multi-purpose mooring, the BTM, and 
is relatively inexpensive. The use of such moorings 
for optical measurements is relatively new, and thus 
several complicating issues need to be considered and 
evaluated. We are working with data from the BTM 
and are doing statistical comparisons of data obtained 
from the BTM and those obtained nearly concurrently 
by the BBOP program. Despite several complicating 
factors, generally favorable correlations have been 
obtained for fundamental radiometric variables. Re- 


sults including estimates of Lu(O-) have been trans- 
ferred to the SeaBASS system. More detailed analy- 
ses and comparisons will continue to be done in the 
future. Particular areas of investigation will include 
calibration, sampling issues (e g., BBOP spatial vs. 
BTM temporal), and biofouling. In addition, it is 
noteworthy that our near real-time data telemetry sys- 
tem has been functioning during the present Deploy- 
ment 10. Lu(O-) estimates for this deployment have 
also been transferred to SeaBASS. 

11.3 WORK PLAN 

It is our intention to place an additional radiome- 
ter at 5-7m depth for a future deployment to determine 
if our estimates of Lu (0-) can be improved signifi- 
cantly. Deployment 10 full bandwidth data will be 
recovered in March 1999 and Deployment 11 will be 
initiated. 
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Chapter 12 

The High-latitude Intercomparison and 
Validation Experiment (HIVE) 

Dave L. Eslinger 

Institute of Marine Science , University of Alaska Fairbanks 
(Present affiliation : NOAA , Charleston , South Carolina) 


12.1 INTRODUCTION 

The HIVE project is a three-year series of cruises 
in both the Bering Sea and Prince William Sound/Gulf 
of Alaska. Both locations are at similar latitudes 
(~60°N) and have similar planktonic communities; 
however, the physical environments are dramatically 
different. These locations provide excellent test sites 
to determine the performance of bio-optical algo- 
rithms at high latitudes. We are making measurements 
of pigments (fluorometric and HPLC), radiometric 
profiles, total suspended matter, and ocean aerosols. 
The HIVE project builds on additional funding from 
the Japanese ADEOS program and the EXXON Val- 
dez Oil Spill Trustee Council, which are supporting 
research projects in the Bering Sea and Prince William 
Sound, Alaska, respectively. 

We are contracted to perform two or more cruises 
per year in the Gulf of Alaska, Prince William Sound, 
and/or Bering Sea depending on availability of sam- 
pling opportunities on US and Japanese vessels. In 
1997, we received additional funding from the NASA 
SIMBIOS program to participate in a collaborative 
Japanese cruise in the Bering Sea prior to the start of 
the actual SIMBIOS project. In the first year of the 
SIMBIOS contract, we participated in three cruises, 
two in Prince William Sound and one in the Bering 
Sea. This technical memorandum reports on all four 
HIVE cruises to date. 

The goal of the HIVE project is to validate high- 
latitude, satellite data products derived from imagery 
from different ocean color sensors. This goal is to be 
accomplished through meeting the following objec- 
tives: 

• Make accurate measurements of in situ data fields 
of potential ocean color data products in northern- 
hemisphere, high-latitude waters. 

• Compare measured data fields with data fields 
determined from ocean color satellites. 

• Intercompare data fields retrieved from different 
ocean color satellites. 

• Intercompare data fields retrieved from the same 
satellite at different times of day. 


The failure of the ADEOS spacecraft precluded 
us from getting any match-up OCTS data during our 

1997 cruise, which was prior to the SeaWiFS launch 
In 1998, clouds, fog and poor weather have prevented 
any match-ups with the SeaWiFS data sets. Therefore 
poor weather and the ADEOS failure have prevented 
us from making any progress towards objectives 2 and 
4. However, we have been successful in making field 
measurements, even under less-than-ideal conditions. 

12.2 RESEARCH ACTIVITIES 

We are making field measurements of phyto- 
plankton pigments using both fluorometric and HPLC 
techniques, profiles of upwelling radiance and down- 
welling irradiance, total suspended matter, and ocean 
aerosols. These activities are described in detail be- 
low. We describe our activities in four cruises The 
first occurred in July, 1997, in the Bering Sea, aboard 
the Japanese training ship Oshoro Maru and has a 
cruise identifier of OS97. The second was in June, 

1998 was in Prince William Sound, AK, aboard a 58- 
foot charter vessel, the M/V Auklet and has a cruise 
identifier of AUK9806. The third was in July and was 
back to the Bering Sea aboard the T/S Oshoro Maru 
and has a cruise identifier of OS98. The fourth HIVE 
cruise was in August, again in Prince William Sound 
aboard the MIWAuklet, and has identifier AUK9808. 

On all cruises, water samples were collected from 
6 depths in the upper water column using Niskin bot- 
tles. During AUK9806, samples were taken from the 
surface and 5, 10, 15, 25 and 50 meter depths. During 
the OS97, OS98, and AUK9808 cruises, sampling 
depths were selected based on standard light levels 
and were generally taken at the depths corresponding 
tol, 10, 20, 36, 50, and 95% of surface light. Depths 
were altered when the water column was too shallow' 
to permit sampling to the one-percent light level. Sur- 
face samples generally came from one-meter depth, 
but were occasionally deeper under rough sea condi- 
tions, when necessary to avoid firing bottles above the 
surface. 
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Phytoplankton Pigments 

Water samples were filtered under low vacuum 
onto GF/F filters, followed by freezing in liquid nitro- 
gen for HPLC, and by grinding and extraction in 90% 
acetone for immediate fluorometric analyses. Chloro- 
phyll and phaeophytin concentrations were measured 
in the acetone extracts aboard the ship immediately 
after centrifugation using the technique of Strickland 
and Parsons (1972). HPLC samples were shipped in 
liquid nitrogen back to Fairbanks, where they were 
stored. Our group does not have the necessary facili- 
ties to perform HPLC analyses, therefore, after stor- 
age, samples from the first three cruises were shipped 
via FedEx to the Grice Marine Lab, Charleston, SC, 
for analysis. A combination of inadequate packing 
and shipment delays in both pick-up and transit led to 
the samples arriving in Charleston at room tempera- 
ture. HPLC analysis of a subset of these filters indi- 
cated significant degradation of pigments had oc- 
curred and the samples were discarded. The filters 
from the fourth cruise remain in Fairbanks, in liquid 
nitrogen, while we determine exactly where the delays 
in handling and timing occurred and how to prevent 
them in the future. 

Total Suspended Mater Concentrations 

Water samples were taken for determination of 
total suspended matter (TSM) according to the JGOFS 
procedures (JGOFS, 1991) as outlined in the SeaWiFS 
protocols (Mueller and Austin, 1995). A modification 
of this was to use pre-combusted glass fiber filters. 
Filtered samples were dried and returned to Fairbanks 
for weighing on a precision analytical balance. 
Weighed dry filters were then recombusted and 
reweighed to give both total suspended matter and 
organic and inorganic fractions. 

Radiance and Irradiance Profiles 

For OS97, we were able to borrow a Biospherical 
Instruments MER-2048 system from the NOAA 
Coastal Services Center. Using this MER, we made 
measurements of E d at 340, 380, 395, 412, 443, 455, 
490, 510, 532, 555, 565, and 671 nm. Eu and U pro- 
files measurements were made at 380, 412, 443, 490, 
510, 555, and 683 nm. For surface reference, E„ was 
measured at 340, 380, 395, 412, 443, 490, 510, 555, 
565, 671, 780, and 875 nm. Temperature was also 
measured during these profiles. 

For AUK9806 cruise Satlantic loaned us an 
SPMR system (serial # 022). The profiler measured 
E d at 405.4, 411.9, 435.6, 442.8, 455.4, 489.6, 510.4, 
531.9, 554.8, 590.5, 665.6, 669.6, 700.4 nm and U at 

405.1, 411.9, 435.1, 443.5, 456.3, 489.4, 510.3, 531.5, 
554.3, 590.5, 665.3, 670.2, and 700.3 nm. The surface 


reference measured L s at 406.5, 412.2, 435.3, 443 4, 
455.9, 489.9, 510.4, 531.6, 554.6, 590.3, 665.1, 670.0, 
and 700.6 nm and E s 405.1, 412.4, 435.6, 442.9, 

456.1, 489.3, 510.4, 531.5, 554.5, 590.4, 664.8, 670.0, 
and 700.6 nm. This system was only used for the 
AUK9806 cruise. Due to failure of a PCMCIA serial 
port card, both the profiler and surface reference could 
not be used simultaneously. Instead, we made a 
measurement of E s and L s> did three profiles, and then 
repeated the surface reference measurements. Given 
the heavy overcast conditions during this cruise, we 
do not expect the surface light values to have varied 
significantly during the time we were profiling. 

For the OS98 and AUK9808 cruises, we used our 
newly purchased SPMR system, which simultaneously 
measured E d , L u , E s temperature, conductivity, and in 
vivo fluorescence profiles. E d was measured at 412.3, 

443.1, 490.7, 510.0, 555.1, 665.6 and 682.4 nm. L u 
was measured 412.3, 442.1, 490.9, 510.1, 554,9, 
665.4, and 682.9 nm. E, was measured at 412.3, 
442.6, 490.3, 510.2, 554.4, 665.0, and 683.8 nm 

Instrument Calibration 

All instruments have been calibrated as required 
by the SeaWiFS protocol and SIMBIOS contract. 
Calibration files have been submitted with the bio- 
optical data to SeaBASS. 

12.3 RESEARCH RESULTS 

Table 1 summarizes the cruise time, location, sta- 
tion, and measurement data. All cruises had complete 
overcast at almost every station. Although these high 
latitudes are generally cloudy, this was an exception- 
ally bad year for our cruises. In 1994 through 1996 
cruises, we had much better weather at these times and 
expected similar intermittent clouds these last two 
years. 

12.4 WORK PLAN 

This last year we participated in three cruises in- 
stead of only the two in which we were contracted to 
participate. Next year, we plan to consolidate our two 
5-day Prince William Sound cruises into one 10-day 
cruise. By increasing our at-sea time to ten days, we 
hope to get the satellite match-ups that we have 
missed up to the present time. At the present, we plan 
to schedule the cruise for middle to late April, to cap- 
ture the high, spring-bloom phytoplankton concentra- 
tions. However, we will review the previous years 
AVHRR data to try and determine if expecting cloud- 
free days during this time period is realistic. We will 
still plan to participate in the T/S Oshoro Maru cruise 
in July. 
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Table 1 . Summary of cruises and in situ data collected. 


Cruise ID 

Location 

Cruise Dates 

Number 
of Stations 

Non-optical 

measurements* 

OS97 

Bering Sea 

18 July -2 Aug. 1997 

20 

Pigments (fluor), 

TSM, T 

AUK9806 

Prince William 
Sound 

5-9 June 1998 

21 

Pigments (fluor), 

TSM, T 

OS98 

Bering Sea 

18 July -1 Aug. 1998 

17 

Pigments (fluor), 

TSM, T, S, Fluor 

AUK9808 

Prince William 
Sound 

8 Aug. - 2 Sept., 1998 

15 

Pigments (fluor & 
HPLC), TSM, T, S, 
Fluor 


* TSM = Total Suspended Matter, T = temperature profile, S = Conductivity profile, Fluor = in vivo fluo- 
rescence profile 
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Chapter 13 

Satellite Ocean Color Validation Using Merchant Ships 

Robert Frouin and David L. Cutchin 
Scripps Institution of Oceanography 
University of California San Diego , La Jolla, California 

Pierre-Yves Deschamps 
Laboratoire d f Optique Atmospherique 
Universite des Sciences et Technologies de Lille , Villeneuve d'Ascq , France 


13.1 INTRODUCTION 

The usual approach to satellite ocean color vali- 
dation is to organize dedicated experiments at sea, e g. 
in regions that are likely to cause atmospheric correc- 
tion errors, and to measure concomitantly water- 
leaving radiance and aerosol optical properties at the 
time of satellite overpass (Clark et al., 1997). The data 
are necessary not only to compare satellite estimates 
with in situ measurements, but also to interpret the 
differences and, eventually, adjust atmospheric cor- 
rection algorithms. These experiments, unfortunately, 
are expensive; they cannot be carried out over the full 
range of expected atmospheric and oceanic conditions. 
Furthermore, because of cloud cover only a few 
match-ups are generally obtained during a campaign, 
making the approach even less cost-effective. 

A complementary approach -and the focus of the 
present project— is to use ships of opportunity, in par- 
ticular, merchant ships traveling routes in the world's 
oceans. For the approach to be cost-effective, the ships 
must participate on a volunteer basis or at a very small 
cost, and the measurement program must not interfere 
with the ship’s activities. Importantly, the ship should 
not be required to stop, i.e. one should be able to col- 
lect the data en route, while the ship is steaming. Be- 
cause of these limitations, one cannot expect to meas- 
ure the full set of ocean-color variables using ships of 
opportunity. Still, to be useful the measurement pro- 
gram should meet the basic requirements for satellite 
ocean color validation. 

It is in this context, and also according to a study 
of the measurements required for satellite ocean-color 
validation (Schwindling et al., 1998), that the SIM- 
BAD concept was developed. SIMBAD is a portable 
radiometer designed for use onboard ships of opportu- 
nity. It measures normalized water-leaving radiance 
(or diffuse marine reflectance) and aerosol optical 
thickness in typical spectral bands of satellite ocean 
color sensors. Knowledge of both variables is needed 
and generally sufficient to evaluate atmospheric cor- 


rection algorithms (Schwindling et al., 1998) The 
instrument is easy to operate, so that any ordinary 
crew can learn quickly how to make measurements. 
To achieve adequate sampling, i.e. acquire data in a 
wide range of oceanic and atmospheric conditions, a 
series of ten radiometers was built for use in two com- 
plementary networks of ships/routes, one operated by 
the Scripps Institution of Oceanography, the other by 
the University of Lille. 

13.2 RESEARCH ACTIVITIES 

The SIMBAD radiometer was designed and built by 
the Laboratoire d’Optique Atmospherique (LOA) of 
the University of Lille. It measures water-leaving ra- 
diance from above the surface as opposed to below the 
surface the usual way. This eliminates constraints as- 
sociated with deploying underwater instruments at 
sea, e.g. stopping the ship. Viewing of the ocean sur- 
face is made from the side of the ship lit by the sun, 
outside the glitter region, at a nadir angle of 45 de- 
grees. A vertical polarizer reduces the skylight re- 
flected by the surface in the instrument’s field-of- 
view. The best location to operate the instrument is 
generally the bow, making it easy to avoid ship- 
generated foam. The radiometer also measures aerosol 
optical thickness, by viewing the sun like a standard 
sun photometer. The same optics (2.5 degree field-of- 
view) and detectors are used in both ocean- and sun- 
viewing modes, but at different electronic gain. 

The radiometric measurements are made simulta- 
neously in five spectral bands centered at 443, 490, 
560, 670, and 865 nm. These bands are not exactly 
those of the ocean color sensors that are flying or 
scheduled to fly, but result from a compromise. They 
constitute, however, a basic set of bands for atmos- 
pheric correction and bio-optical algorithms. Viewing 
of the ocean and sun is made sequentially, but the 
non-simultaneity of the measurements is not limiting. 
SIMBAD does not measure incoming solar irradiance, 
a variable that must be known to normalize water- 
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leaving radiance. It is calculated accurately in clear 
sky conditions or when the sky is partly cloudy 
(<30%) with the sun not obscured by clouds. This 
defines the sky conditions under which SIMBAD 
should be operated. In addition to radiometry, SIM- 
BAD automatically acquires viewing angles (using a 
magnetometer and two inclinometers), time and geo- 
graphic location (using a GPS), pressure, and tem- 
perature. It requires about 15 minutes to collect a 
complete data set (ocean, sun, optical zero), including 
deploying the instrument and logging ancillary data 
(wind speed, sea state, cloud cover, etc.). Internal 
memory and batteries allow for 3 months of opera- 
tions in normal mode (typically 5 sets of measure- 
ments per day). At the end of the campaign or voyage, 
after the ship returns to port, the data are downloaded 
onto computer disk and diskette and the radiometer is 
calibrated. In the case of long campaigns, the data are 
mailed from selected ports along the ship’s route, al- 
lowing for a more rapid data turnover. 

Justification and Verification 

The principle of the SIMBAD radiometer, in par- 
ticular its ability to reduce reflected skylight with the 
help of a vertical polarizer, has been justified theoreti- 
cally and verified experimentally (Fougnie et al., 
1999b). Theoretical calculations show that, for a nadir 
angle of 45 degrees and a relative azimuth angle be- 
tween solar and viewing directions of 135 degrees, the 
recommended geometry, reflected skylight is reduced 
to typically 10' 3 in reflectance units at 443 nm. This 
represents 2 to 10% of the diffuse marine reflectance, 
the signal of interest. Furthermore, the effects of sur- 
face roughness on skylight reflection and, hence, un- 
certainties in sea state (wind speed) are minimized. 
Taking into account typical uncertainties in wind 
speed and geometry, the residual reflected skylight are 
correctable to a few lO' 4 in reflectance units. For most 
oceanic waters, the resulting error on the diffuse ma- 
rine reflectance in the blue and green is less than 1%. 
The theoretical results have been verified experimen- 
tally at the Scripps Institution of Oceanography pier, 
by viewing the ocean surface with a scanning polari- 
zation radiometer (Fougnie et al., 1999b). The various 
angular and spectral effects predicted by theory have 
been evidenced in the measurements. 

Since the water body polarizes incident sunlight 
(e g., Ivanoff, 1975, Frouin et al., 1994), the measured 
polarized radiance differs from the total radiance. For 
the recommended geometry, the underwater scattering 
angle varies only between 148 and 158 degrees, and 
the maximum polarization factor varies between 0.83 
and 0.92 (Fougnie et al., 1999b). One should be able 
to correct the polarization effects to within 5% relative 
accuracy, which compares with other errors, such as 
those due to bi-directional effects and radiometric 
calibration. 


The SIMBAD radiometer was evaluated at sea, 
during California Cooperative Fisheries Investigation 
(CalCOFI) cruises in the southern California Bight 
The diffuse marine reflectance measured by the radi- 
ometer was compared with that measured by an un- 
derwater MER optical system. The agreement be- 
tween the two types of measurements was good, to 
within 1.6 10' 3 in reflectance units. Aerosol optical 
thickness measured by the SIMBAD radiometer also 
compared well with that measured by sun photometers 
during the October-November 1996 CalCOFI cruise 
(Nakajima et al., 1999) and during the Second Aerosol 
Characterization Experiment (ACE -II). 

The SIMBAD radiometer proved useful to check 
vicariously the radiometric sensitivity of the Polariza- 
tion and Directionality of the Earth’s Reflectance 
(POLDER) instrument onboard ADEOS (Fougnie et 
al., 1999a). The satellite radiance was compared with 
that computed for the same geometry using a radia- 
tion-transfer model. SIMBAD measurements of aero- 
sol optical thickness and marine reflectance at the time 
of satellite overpass were used as input to the model 
The accuracy of the vicarious calibration coefficients 
was estimated to be better than 3%. A large decrease 
of the POLDER instrument response was found in the 
blue, confirming the results previously obtained using 
alternative calibration techniques. 

Data Format 

The SIMBAD data sets are fairly complete 
for satellite ocean color validation. They include near- 
concurrent values of spectral aerosol optical thickness 
and diffuse marine reflectance, necessary variables for 
checking radiometric calibration and evaluating at- 
mospheric correction algorithms. They also include 
fractional cloud coverage and wind speed (useful to 
estimate whitecap reflectance). The data sets can be 
accessed at the following web addresses: 
http://genius.ucsd.edu/~frouin/ and http://www-loa- 
univ-lille 1 ,fr/~poteau/. 

13.3 RESEARCH RESULTS 

Since the launch of the SeaWiFS in September 
1997 and until December 1998, SIMBAb measure- 
ments have been made during 12 research cruises and 
9 merchant ship voyages (Table 1). A total of 156 
SIMBAD data sets have been collected under clear or 
mostly clear sky (fractional cloud coverage less than 
0.2) within a few hours of satellite overpass A variety 
of oceanic and atmospheric conditions have been 
sampled, but no data have been acquired in the South 
Atlantic and Indian Oceans. The research cruises and 
merchant ship voyages are listed below in chronologi- 
cal order with location or route, date, name of the 
SIMBAD operator, and network (SIO or LOA). 
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Table 1, SIMBAD cruise measurements from 9/1997 to 12/1999 


CalCOFI 9709, Southern California Bight, 09/20/98-10/08/98, Wieland, SIO. 

M TV Asian Challenger, Valparaiso-Auckland, 10/12/98-10/18/98, Cutchin, SIO. 

R/V Melville, San Diego- Acapulco-Callao-Valparaiso-Easter-Papeete-Honolulu, 1 l/02/97-06/07/98,Comer, SIO. 
M/VNacre, Honolulu-Coronel, 12/20/98-01/02/98, Cutchin, SIO. 

CalCOFI 9802, Southern California Bight, 01/23/98-02/10/98, Subranamian, SIO. 

CalCOFI 9804, Southern California Bight, 04/02/98-04/24/98, Mitchell, SIO. 

R/V Revelle, Eastern Tropical Pacific Ocean, 04/09/98-04/27/98, Dupouy, LOA. 

M/V Toucan, Le Havre-Cayenne-le Havre, 05/04/98-05/23/98, Poteau, LOA. 

M/V Chevron Mississippi, Honolulu- Valdez, 05/18/98-05/23/98, Cutchin, SIO. 

M/V Toucan, Le Havre-Cayenne-le Havre, 06/08/98-06/20/98, Poteau, LOA. 

CalCOFI 9807, Southern California Bight, 07/09/98-07/27/98, Wieland, SIO. 

ATMARA 9807, English Channel, 07/21/98-08/06/98, Wartel, LOA. 

NAN 9807, Cape Cod, 07/06/98-07/10/98, Shieber, SIO. 

M/V Mokihana, Oakland-Guam-Taiwan-Japan-Oakland, 07/20/98-08/1 5/98,Bousseloub, SIO. 

CalCOFI 9809, Southern California Bight, 09/13/98-10/01/98, Wieland, SIO. 

BREST 9809, Northeastern Atlantic Ocean, 09/25/98-10/01/98, Gohin, LOA. 

M/V Mokihana, Oakland-Guam-T ai wan- Japan-Oakland, 08/25/98-09/08/98, Bousseloub, SIO. 

M/V Micronesian Navigator, Califomia-Micronesia-Philippines-Califomia, 08/25/98- 10/28/98, Sumogat, SIO. 
R/V Cape Hatteras, Beaufort- Jacksonville, 10/21/98-12/01/98, Subranamian, SIO. 

M/V Micronesian Navigator, Cahfornia-Micronesia-Philippines-California, 1 0/30/98- 1 2/31/98, Sumogat, SIO. 
R/V Ka’imimoana, Honolulu-Fiji-Honolulu, 10/19/98-12/11/98, Strutton, SIO. 


13.4 WORK PLAN 

A collaborative LOA-SIO measurement program 
for satellite ocean color validation based on ships of 
opportunity has been developed. This has been ac- 
complished in a series of successive, logical steps. 
First, a study of validation requirements was per- 
formed, indicating that measurements of spectral wa- 
ter-leaving radiance and aerosol optical thickness at 
the time of satellite overpass are generally sufficient to 
evaluate satellite-derived ocean color (Schwindling et 
al., 1999). Second, in view of the study and limitations 
of dedicated validation experiments at sea, the SIM- 
BAD radiometer was designed, for use of ships on 
opportunity such as merchant ships traveling regular 
routes in the world’s oceans. Third, the instrument’s 
innovative principle (use of a vertical polarizer to re- 
duce skylight reflection) was verified experimentally 
at the Scripps Institution of Oceanography pier (Foug- 
nie et al., 1999b). Fourth, the instrument was tested at 
sea, on station and en route, during various oceano- 
graphic campaigns (Fougnie, 1998, Nakajima et al., 
1999). Fifth, suitable merchant ships and routes in the 
Pacific and Atlantic Oceans were identified and se- 
lected officers were trained. Some of the ships, for 
example M/V Micronesian Navigator (SIO network) 
and M/V Toucan (LOA network), are now providing 
SIMBAD data routinely. Sixth, SIMBAD data were 
used to check vicariously the radiometric calibration 
of the POLDER instrument onboard ADEOS (Fougnie 
et al., 1999a). 

Regarding measurement opportunities, the pro- 
gram’s focus has been on research campaigns during 


which water-leaving radiance and/or aerosol optical 
thickness were measured by other instrumentation 
(e.g., CalCOFI campaigns), so that comparisons could 
be made. Although more evaluation of the SIMBAD 
data is necessary, and should be done on a continuous 
basis, in order to increase the number of sur- 
face/satellite match-ups the focus should be shifted to 
ships that do not cany out the same type of measure- 
ments as the SIMBAD radiometer (e g., more mer- 
chant ships). The program also needs to be extended 
to ships traveling in other oceanic regions, in particu- 
lar the southern Oceans and the Indian Ocean, to com- 
plement the present coverage. There are many possi- 
bilities, and the demand for SIMBAD radiometers has 
been increasing. The realization of a new series of 
radiometers is envisioned, with the instruments avail- 
able for the next generation of satellite ocean color 
sensors. 
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Chapter 14 

High Attitude Measurements of Radiance at High Spectral 
and Spatial Resolution for SIMBIOS Sensor 
Calibration, Validation, and Intercomparisons 

Robert Green and Thomas G. Chrien 
NASA, Jet Propulsion Laboratory 
Pasadena, California 


14.1 INTRODUCTION 

The successful combination of data from different 
ocean color sensors depends on the correct interpreta- 
tion of signal from each of these sensors. Ideally, the 
sensor measured signals are calibrated to geophysical 
units of spectral radiance, and sensor artifacts are re- 
moved and corrected. The calibration process re- 
samples the signal into a common radiometric data 
space so that subsequent ocean color algorithms that 
are applied to the data are based in physical processes 
and inherently sensor independent. The objective of 
this SIMBIOS effort: “High Altitude Measurements of 
Radiance at High Spectral and Spatial Resolution For 
SIMBIOS Sensor Calibration, Validation, and Inter- 
comparisons is to conduct a series of direct calibration 
comparisons using the NASA Airborne Visi- 
ble/Infrared Imaging Spectrometer (AVIRIS),” is to 
assess the calibration of the satellite sensors which 
participate in SIMBIOS. 

_ The AVIRIS se nsor fl ies at 20 km altitude on a 
NASA ER-2 near the top of the absorbing and scat- 
tering portions of the atmosphere AVIRIS measures 
the solar reflected spectrum from 370 nm to 2500 nm 
through 224 contiguous 10 nm wide spectral channels. 
Data are acquired as images of 11 by up to 800 km 
extent with 20 m spatial resolution. The stability and 
repeatability of AVIRIS calibration has been and con- 
tinues to be validated through a series of inflight cali- 
bration experiments. With pre and post flight calibra- 
tions of AVIRIS, coupled with the on-board calibra- 
tor, calibration accuracy of better than 2% spectral, 
3% radiometric and 3% spatial have been achieved. 

The activities and results presented here are pri- 
marily for SeaWiFS. However, related work has been 
carried out for ADEOS OCTS. In the future, it is 
planned that AVIRIS will underfly EOS. 

14.2 RESEARCH ACTIVITIES 

On the 2nd of October 1997 AVIRIS under flew 
the SeaWiFS sensor of the coast of Southern Califor- 
nia near latitude 35° and longitude 127°. The 


SeaWiFS sensor was tilted to the Northeast to avoid 
sun glint. In order for the AVIRIS spectral images to 
overlap the area, the azimuth, and the zenith angles of 
the SeaWiFS measurements, AVIRIS was flown in a 
circle beneath SeaWiFS. The AVIRIS scan angle (± 
15 degrees) plus the aircraft roll angle (20 degrees) 
assured overlap of both area and observation geome- 
try. The AVIRIS data from the 2nd of October 1997 
flight were spectrally, radiometrically and spatially 
calibrated through the AVIRIS data system algo- 
rithms. The SeaWiFS data and the SEADAS proc- 
essing software were requested and received. Both 
the AVIRIS and the SEAWIFS data sets include the 
position and pointing information at the time of acqui- 
sition. The simultaneous acquisition of the AVIRIS 
and SeaWiFS with overlapping observation geometry 
provides the essential data set to complete the calibra- 
tion objective. 

14.3 RESEARCH RESULTS 

Determination of the areas of the SeaWiFS and 
AVIRIS images with the same observation geometries 
requires projection of the data based on the position 
and pointing knowledge of the sensor. A set of algo- 
rithms were developed to accept the SeaWiFS position 
and point information and calculate the azimuth and 
zenith angles for the 2nd of October SeaWiFS data. 
All angles were calculated from the perspective of the 
surface. 

The observation geometry at the surface is the 
correct geometry for calibration comparison. AVIRIS 
receives and records both Global Positioning System 
(GPS) and Inertial Navigation System (INS) informa- 
tion from the ER-2 aircraft. The GPS provides lati- 
tude, longitude and altitude information, while the 
INS provides roll, pitch, and yaw information. These 
position and pointing data were used with an analo- 
gous set of algorithms to calculate the AVIRIS spec- 
tral image azimuth and zenith angles were calculated 
With both the AVIRIS and SeaWiFS image azimuth 
and zenith angles determined, the target area of over- 
lapping observation geometry was identified. This 
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area contains the AVIRIS spectra and SeaWiFS multi 
spectral data of the same area on the surface with the 
same observation azimuth and zenith (Table 1). 

The average AVIRIS spectrum from the area of 
observational and temporal coincidence was extracted. 
A correction factor was applied to the AVIRIS spec- 
tra for the transmittance from the 20 km AVIRIS alti- 
tude to the top of the atmosphere. This correction 
factor was calculated with the MODTRAN radiative 
transfer code using Ozone values from the TOMS 
total ozone archive for the area. The SeaWiFS multi 
spectral band data were extracted directly for this tar- 
get. A graphical comparison of the AVIRIS spectrum 
and SeaWiFS bands showed good agreement. In order 
to quantitatively compare the AVIRIS and SeaWiFS 
data, a convolution algorithm was developed to con- 
volved the AVIRIS spectra to the SeaWiFS bands. 
The algorithm was applied to the AVIRIS spectra for 
the target area and the AVIRIS predicted radiances 
compared to the SeaWiFS reported radiances (Table 
2). These results are preliminary. A review of the 
version of SEADAS used to process the SeaWiFS data 
is underway. A high priority activity has begun to 
review and compare the absolute radiometric standard 
used for SeaWiFS calibration and AVIRIS calibration. 

The primary objective of this SIMBIOS effort has 
been achieved. On the 2nd of October 1997, AVIRIS 
underflew the SeaWiFS sensor. The data set was ac- 
quired such that the observation geometries of 
SeaWiFS and AVIRIS could be match for direct com- 
parison of measured radiance. A set of algorithms 
were developed to calculate the surface azimuth and 
zenith angles for both the SeaWiFS and AVIRIS data 
sets. These algorithms required as input position and 
pointing information for the two sensors. The initial 
graphical comparison of the AVIRIS spectrum and 
SeaWiFS bands showed good agreement for the ocean 
target of overlapping observation geometry. An addi- 
tional set of algorithms were developed to convolved 
the AVIRIS spectra to the SeaWiFS spectral bands. 
Direct comparison of the AVIRIS measured and 
SeaWiFS measured radiance shows the agreement is 


best for SeaWiFS bands 3,4,5, and 6 and least good 
for SeaWiFS bands 1, 2, 7, and 8. These are initial, 
preliminary results and additional review and uncer- 
tainty analysis is required to strengthen the basis of 
comparison. Never-the-less, these initial results vali- 
date the objective, approach, and activities of this 
SIMBIOS Project. 

14.4 WORK PLAN 

In 1999 three additional underflights of 
SEAWTFS are planned. The first underflight is 
planned for the Spring of 1999 from the west coast. A 
second underflight will be attempted from the east 
coast in the late Spring or early Summer. The third 
underflight is planned for the west coast in the late 
summer early fall time frame. Each underflight will 
include two flight circle acquisitions to assure a match 
to the SeaWiFS viewing geometry. The first flight 
circle will be timed to coincide with the satellite over- 
pass. The second flight circle acquisition enables in- 
vestigation of the time sensitivity of the measure- 
ments. Efforts will be made to coordinate all under- 
flights with other participating SIMBIOS investiga- 
tors. Following each successful underflight, the 
AVIRIS runway calibration standard will be used to 
place the maximum accuracy on the AVIRIS radio- 
metric calibration for the SeaWiFS underflight 

In April of 1999, the new SIMBIOS SXR-II 
transfer radiometer will be used at JPL to further 
compare the SIMBIOS and AVIRIS radiometric stan- 
dards. Results of all the 1999 activities of this project 
will be provided to SIMBIOS. Work will proceed on 
writing up the 1998 and 1999 results for publication. 
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Table 1 . Preliminary results of AVIRIS and SeaWiFS radiance comparison. 


BAND 

SeaWiFS 

AVIRIS 

AVIRIS 

Difference 


Radiance 

Radiance 

Uncertainty 


1 

7.800 

6.791 

5.1% 

12.9% 

2 

6.852 

6.311 

3.5% 

7.9% 

3 

4.899 

4.861 

2.1% 

0.8% 

4 

3.995 

3.972 

2.5% 

0.6% 

5 

2.753 

2.711 

2.9% 

1.5% 

6 

1.248 

1.193 

2.3% 

4.5% 

7 

0.641 

0.601 

3.6% 

6.2% 

8 

0.449 

0.412 

3.5% 

8.3% 


The units of radiance are microWatts/cm 2 /nm/sr. 


Table 2. Website locations for AVIRIS SeaWiFS match up example 


1 First AVIRIS flight circle, ftp://makalu.jpl.nasa.gov/pub/rog/seawifs/pic01.jpg 

2 Corresponding SeaWiFS image ftp://makalu.jpl.nasa.gov/pub/rog/seawifs/pic04.jpg 
SeaWiFS azimuth ftp://makalu.jpl.nasa.gov/pub/rog/seawifs/pic05.jpg 
SeaWiFS zenith ftp://makalu.jpl.nasa.gov/pub/rog/seawifs/pic06.jpg 
AVIRIS azimuth ftp://makalu jpl.nasa.gov/pub/rog/seawifs/pic07.jpg 
AVIRIS zenith ftp://makalu jpl.nasa gov/pub/rog/seawifs/pic08 jpg 
AVIRIS and SeaWiFS data measured within 5 seconds. 

ftp : //makalu . jpl . nasa . gov/pub/rog/seawifs/pic09 . jpg 
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Chapter 15 

Merging Ocean Color Data from Multiple Missions 

Watson W. Gregg 

NASA/Goddard Space Flight Center 
Greenbelt , Maryland 


15.1 INTRODUCTION 

In response to the potential importance of phyto- 
plankton in the global carbon cycle and the lack of 
comprehensive data, the National Aeronautics and 
Space Administration (NASA) and the international 
community have established high priority satellite 
missions designed to acquire and produce high quality 
ocean color data. Seven of the missions are routine 
global observational missions: the Ocean Color and 
Temperature Sensor (OCTS), the Polarization and 
Directionality of the Earth's Reflectances sensor 
(POLDER), SeaWiFS, Moderate Resolution Imaging 
Spectrometer- AM (MODI S- AM), Medium Resolution 
Imaging Spectrometer (MERIS), Global Imager 
(GLI), and MODIS-PM. In addition, there are several 
other missions capable of providing ocean color data 
on smaller scales. 

Any individual ocean color mission is limited in 
ocean coverage due to sun glint and clouds. For ex- 
ample, one of the first proposed missions, the 
SeaWiFS, can provide about 45% coverage of the 
global ocean in four days and only about 15% in one 
day (Gregg and Patt, 1994). 

The proliferation of missions suggests that better 
ocean coverage can be obtained in less time if the data 
are combined. In addition to improved, and faster 
coverage, data can be taken from different local times 
of day if the missions are placed in different orbits, 
which they are. This can potentially lead to informa- 
tion on diel variability of phytoplankton abundances. 

We propose to investigate, develop, and test algo- 
rithms for merging ocean color data from multiple 
missions. We seek general algorithms that are appli- 
cable to any retrieved ocean color data products, and 
that maximize the amount of information available in 
the combination of data from multiple missions. Most 
importantly, we will investigate merging methods that 
produce the most complete coverage in the smallest 
amount of time, nominally, global daily coverage. We 
will also assess the ability to produce fuller coverage 
in larger time increments, including 4-day, 8-day 
(weekly), monthly, seasonal and annual. Secondarily, 
we will investigate the ability of the missions to pro- 
duce coverage at different times of day, for diel vari- 
ability and dynamical evaluations, and develop algo- 
rithms to produce this information. 


15.2 RESEARCH ACTIVITIES 

Work in the first two years has focused on: a) 
defining the problems and opportunities provided by 
the existence of multiple sensors; b) analysis of past 
sensors to provide insights into the characteristics, 
drawbacks, and advantages of individual sensor re- 
sponses in the context of merging; and c) analysis, 
development, and testing of candidate merger algo- 
rithms. 

a) Defining the Problems 

First we assessed the global coverage improve- 
ments possible by the six global missions These re- 
sults show that significant scientific advantages can 
accrue from assembling and merging data from the 
multiple satellite platforms proposed for ocean color 
in the next decade. The principal advantages are in- 
creased ocean coverage in less time, and new obser- 
vations of the daily dynamics of phytoplankton abun- 
dances, resulting from different observation times of 
co-located ocean areas. Data from three satellites as 
opposed to one, can increase ocean coverage by 58% 
for one day, and 45% for four days. Additional satel- 
lites produce diminishing returns, however. This latter 
point is not necessarily an adverse finding, since each 
mission has a limited life expectancy and in-flight 
problems are, unfortunately, still not rare in the satel- 
lite business. Using observations from pairs of mis- 
sions, as much as 14 hour time differences at co- 
located points can be realistically achieved at high 
latitudes, even considering the distribution of land 
masses and ice cover. Smaller time differences are 
observed at lower latitudes, but 5 to 7 hour differences 
are still available. Furthermore, massive numbers of 
these co-located pairs are available, suggesting that 
routine scientific studies of diel phytoplankton vari- 
ability can be supported (Gregg et al., 1998). More 
detailed analyses emphasized seasonal and regional 
coverage improvements and was limited to the 
SeaWiFS/MODIS combination, since these are the 
two missions planned for next launch. The results 
showed that the launch of EOS AM-1 provides an 
opportunity to improve ocean color observations by 
combining data from the SeaWiFS and MODIS mis- 
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sions. The sensors have different scanning character- 
istics and are flying on different platforms in different 
orbits. The results suggested that very large im- 
provements in coverage frequency (daily to four-day) 
can be obtained by combining data from both sensors: 
40-47% increases in global coverage over SeaWiFS 
alone in one day, and > 100% in areas near the solar 
declination. Four-day increases are slightly smaller 
for global coverage, 29-35%, but meridional percent 
increases are similar to the one-day case. These dif- 
ferences are due to reduced sun glint contamination 
obtained by tilting (SeaWiFS), and scanning away 
from the maximum glint region (MODIS),due to its 
10:30 AM equator crossing time, and due to the large 
scan width of MODIS. The results show that 
SeaWiFS and MODIS are very complementary ocean 
color missions that can provide more complete obser- 
vations of ocean processes at high frequency if data 
are combined (Gregg and Woodward, 1998). A more 
detailed paper with more fully defined results was 
published as a NASA TM (Woodward and Gregg, 
1998). Third, assessed potential capability for im- 
proving ocean color observations by selecting com- 
plementary mission orbits. Considering that observa- 
tions are severely hampered by cloud cover and sun 
glint, a possibility exists for using multiple missions to 
improve the coverage of the oceans in shorter time 
scales. In fact, only about 10-18% of the oceans are 
observed in a single day, even by so-called global 
observational missions, due to these two ocean color 
contaminants. Analysis of a 7-day period of cloud 
cover from the International Satellite Cloud Climatol- 
ogy Project (ISCCP), show that 12% new surface area 
is available for viewing each day. This translates to 
an increase of about 0.5% per hour of separation in 
viewing times. Thus if the ocean could be viewed by 
2 different satellites 4 hours apart, 2% more ocean 
area could be observed. This represents a coverage 
increase of about 13%. If 2 satellites were placed in 
Earth orbit with equator crossing times 4 hours apart, 
the improvement in coverage by these 2 satellites over 
a single one is about 60%. Furthermore, by managing 
the orbits of the 2 satellites, nearly complete sun glint 
avoidance can be achieved, further improving the 
ocean coverage. The best improvements can be made 
by 2 satellites in the same node, whose orbital posi- 
tions are adjusted to view the sun glint contaminated 
areas of the other. If scan edges are useful for quanti- 
tative ocean color observations (the validity of which 
is unknown at this time), then only 2 polar-orbiting 
satellites are necessary. If not, a third satellite is nec- 
essary, preferably in a low inclination orbit where 
losses in coverage in the tropics by the polar orbiters 
can be best compensated (polar orbiters overlap cov- 
erage at high latitudes, and the scan edges are neces- 
sary only near the equator). However, the best con- 
figuration is 3 geostationary satellites, which provide 
complete global coverage routinely with a viewing 


time separation that can maximize cloud and sun glint 
avoidance, with a single polar orbiter to provide high 
latitude coverage (Gregg, 1999a). 

b) Analysis of Past Sensors 

We continued investigations into OCTS data. 
These investigations began before the failure of 
ADEOS, in an effort to characterize OCTS data for 
use with a merging activity with SeaWiFS. The loss 
of data from the sensor precludes its use in merging, 
but the similarity of some aspects of sensor design 
with MODIS suggests that significant understanding 
of the capabilities and deficiencies of MODIS data can 
be facilitated through a thorough analysis of OCTS 
data. In fact, we have found this to be true. Initial 
observations just after failure of OCTS suggested 6 
problem areas: band registration, image striping, cloud 
noise, navigation, calibration, and bright target re- 
sponse. Analysis of the first three led to solutions that 
improved imagery substantially (Gregg, 1998; Gregg, 
1999b). Particularly the methods for reducing image 
striping is a useful method for a problem expected to 
occur in MODIS imagery. Analysis of geolocation, 
radiometric stability and accuracy, and bright target 
responses were characterized and submitted for publi- 
cation (Gregg et al.,1999). The results here also have 
implications for MODIS. 

c) Merging Algorithms 

Efforts have began on developing merging algo- 
rithms. We are investigating 4 possible approaches: a) 
simple splicing, where 1 data set is deemed superior 
and the other is patched into gaps, b) sensor weight- 
ing, where specific dependences and deficiencies are 
identified and co-located observations are merged 
using different weighting functions for the sensors, c) 
the Conditional Relaxation Analysis Method, where 
the best data are selected as interior boundary condi- 
tions into a merged set using Poisson's equation (Rey- 
nolds, 1988), and d) optimum interpolation, where 
merging occurs by weighting individual sensor data to 
minimize spatial covariance function. The Condi- 
tional Relaxation Analysis Method (CRAM) was im- 
plemented and tested using CZCS data and in-situ 
data. Although these results were funded by separate 
grant, the results enabled us to begin to understand the 
strengths and weaknesses of the method and its poten- 
tial use for satellite data (Gregg and Conkright, 1 999). 
Several modifications are required before effective use 
with in-situ data, but it is not clear at this time whether 
these modifications will be necessaiy for merging of 
multiple satellite data, since data paucity may not be 
an issue. It may be an issue, however, if significant 
deficiencies are identified. 
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Chapter 16 


Validation of Bio-optical Properties in Coastal 
Waters: a Joint NASA - Navy Project 

Richard L. Miller 

NASA, Stennis Space Center , Mississippi 


16.1 INTRODUCTION 

This project is a joint effort between NASA’s 
Earth System Science Office and the Naval Oceano- 
graphic Office (NAVOCEANO) at the John C. Sten- 
nis Space Center, MS. The collaboration reflects a 
partnership developed by NASA and NAVOCEANO 
over six years to share resources and expertise in the 
collection of ocean optics, hydrographic, biological, 
and remotely sensed data to support each agency’s 
goals in coastal research. NAVOCEANO provides 
ship time, instrumentation, and operational support 
during oceanographic surveys while NASA provides 
expertise in data collection, computer programming 
and modeling, and data processing and analysis. This 
SIMBIOS project was established to provide data 
collected in the South China Sea, Gulf of Thailand and 
adjacent coastal waters. This project is designed to 
provide NASA with cost-effective access to in-water 
optics and ancillary data in coastal waters representing 
diverse turbid environments, and if successful, will 
help provide an effective mechanism for interagency 
cooperation and data exchange. 

16.2 RESEARCH ACTIVITIES 

NAVOCEANO hydrographic surveys 260797 (27 
April - 19 May) and 260897 (23 May - 16 June) were 
completed in the South China Sea, and the Gulf of 
Thailand / South China Sea, respectively. The SIM- 
BIOS project investigators comprise the optics team 
of these surveys, one of three primary operations: op- 
tics, bioluminescence, and basic hydrography. Vertical 
profiles were acquired by each operations team at a 
dense grid of selected stations. In addition, XBTs were 
dropped at regular intervals along the cruise track 
between the profiling stations. Three optics packages 
were deployed. These were designated as the day 
(IOPs and radiometer), night (IOPs), and the Satlantic 
system. The day and night packages (Table 1) used the 
Mermaid II data acquisition system (Miller et al., 
1995). System power and communications network 
for the day and night packages were provided by 
MOD APS (Wet Labs, Inc). 


Vertical Profiles and IOP 

Several power and communication problems were 
encountered with the MOD APS deck and subunits of 
each package. Power delivery from the subunit to the 
MER 1032 failed after 1 1 profiles and hence the com- 
bined package was abandoned. The subunit for the 
IOP package could not establish communications with 
the Seabird CTD. The CTD was disconnected from 
the MODAPS system and was deployed in standard 
profiling mode (internal data collection). CTD Data 
were downloaded independently following each cast. 
Unfortunately, this approach decoupled the CTD and 
AC-9 data streams. A program was developed and 
used to interactively merge AC-9 and CTD profiles 
during preprocessing and quality control. Two AC-9s 
were on the IOP package. The inflow tubes of one 
AC-9 were fitted with 0.2 pm inline filters at several 
stations to partition the total spectral absorption and 
beam attenuation into dissolved and particulate frac- 
tions. At these stations, total (unfiltered) and dissolved 
(filtered) a and c profiles were obtained. A compari- 
sion / calibration of the AC-9s were performed when 
both AC-9s were unfiltered. The IOP package was 
deployed at 54 stations (1 1 1 profiles) and 120 stations 
(248 profiles), during survey 260797 and survey 
260897, respectively. Two PC-based programs were 
written in IDL 5.0 (castDisplay, castMapper) to apply 
quality control checks and routine data correction and 
display. Temperature, salinity and scattering correc- 
tions were applied to all AC-9 data Scattering correc- 
tions used either Method 1 of the Zaneveld Method 
(AC-9 User Guide, Wet Labs, Inc). A branching algo- 
rithm was used based on the a715 nm / b7 1 5 nm ratio 
to select which scattering method was applied. 
Method 1 was used for ratios greater than 0.35. The 
Zaneveld Method was used for ratio values <= 0.35. 
The castDisplay and castMapper (a geographic con- 
touring package for vertical profiles) will be used on 
subsequate surveys for near-realtime processing and 
analysis. All data have been translated to SeaBASS 
format. 
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Satlcmtic Profiles 

A SeaWiFS Profiling Multichannel Radiometer 
(SPMR) and a SeaWiFS Multichannel Surface Refer- 
ence (SMSR) instrument were deployed at 42 and 36 
stations during survey 26079 and 260897, respec- 
tively. In general, three profiles were acquired at each 
station. The central wavelengths for the SMPR and 
SMSR channels were: 405.3, 412.9, 435.1, 443.8, 
456.2, 490.9, 509.1, 532.7, 555.2, 590.5, 665.8, 683.8, 
and 700.4 nm. Satlantic’s Prosofit 4.0 (revision f) was 
used to produce derived products of K, Reflectance, 
Rjs, L wn , and modeled pigments. The Satlantic system 
performed without problems. 

Discrete Water Samples 

At most stations, discrete water samples were 
collected at select depths (based on profiling transmis- 
someter and fiuorometer data) using Niskin Bottles 
mounted on the hydrographic CTD rosette. These 
water samples were filtered for analysis of fluoromet- 
ric chlorophyll a and phaeopigments, HPLC pigments, 
CDOM absorption, particle absorption (QFT), and 
total suspended mater. All samples were collected and 
processed following the procedures of Mueller and 
Austin (1995). HPLC pigments determined were chlo- 
rophyll a, chlorophyll b, peridinin, fucoxanthin, 19 
hex, prasinoxanthin, alloxantin, lutein, zeaxanthin, 
diadinoxanthin, canthaxantin, and P-carotene. Ab- 
sorption spectra were determined for total, detrital, 
and phytoplankton components following the hot- 
methanol extraction technique (Kishino et al , 1985). 
Corrections for path-elongation were performed fol- 
lowing the method of Mitchell (1990). CDOM ab- 
sorption spectra were collected using a Perkin-Elmer 
Lambda 3 spectrophotometer. The number of samples 
collected for each analysis were (Survey 260797; 
260897): HPLC (71; 100), fluorometic chi a (86; 149), 
total suspended mater (50; 50), CDOM absorption (5; 
20), and particle absorption (0; 61). 

16.3 RESEARCH RESULTS 

The project was able to follow the proposed work 
plan with regards to data collection and processing. In 
general, instrument calibration and data collection 
followed the guidelines provided by the proposal with 
the exception that the MER 1032 needed to be phased 
out of the instrument suite. The problems with the 
MOD APS created unexpected delays in preprocessing 
and analysis. In addition to collection problems and 
limitations, considerable effort was invested to de- 
velop a reliable and robust method to merge CTD and 


AC-9 data. Raw and intermediate data products were 
created for most profiler data. The volume of data 
collected at sea exceeded our expectations. This re- 
sulted in a significant increase in the time required to 
process data. The NAVOCEANO provide operational 
support to the US fleet by conducting military surveys 
worldwide. These surveys are conducted under the 
provisions of military ordinances, and as such, do not 
require diplomatic permission from local govern- 
ments. 

The data collected during this project were mostly 
obtained in the Exclusive Economic Zone (EEZ) of 
several countries. Therefore, as a consequence of 
“The Law of the Sea”, these data are currently re- 
stricted to DOD and DOD contractors use only. 
Hence, specific results and findings can not be pre- 
sented at this time. 

Members of the project team at NASA and NA- 
VOCEANO are actively engaged in determining how 
these data may be made available to the SIMBIOS 
project office. It is the project’s intent to establish a 
precedent for the release and use of this data for civil- 
ian projects without compromising the Navy’s inter- 
national agreements and mandates. If successful, this 
project may establish the means by which an impor- 
tant resource of oceanographic data may be made 
available to the research community. 

16.4 WORK PLAN 

The project will follow the proposed project plan 
with the exceptions that above water estimates of Rrs 
(A.) will be obtained using a GER (Geophysical Envi- 
ronmental Research, Inc.) 1500 hand-held spectrora- 
diometer (300 -100 nm, effective 1 nm bandwidth) 
and the MER 1032 will not be used. 

A survey to the South China Sea and Gulf of 
Thailand is planned for the summer 1999. A duplicate 
Satlantic system has been purchased to provide system 
redundancy. A Wet Labs WetStar profiling fiuorome- 
ter will be configured on a profiling package to direct 
the selection of depths for discrete water samples. 
SeaWiFS LAC data will be requested during filed 
activites and upon processing using SeaDAS will be 
compared with in situ observations. 

The project has acquired a AC-9 to ensure cali- 
bration history and performance for project deploy- 
ments. A MicroTops II sun photometer will be used 
during future cruises. Increased effort will be devoted 
to secure the release of the data for SIMBIOS project 
use. 
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Table 1. Profiling Instrument Packages. 


INSTRUMENTS 

MEASUREMENTS 

Day Package 


MER 1032 In-water Radiometer 1 

Spectral Upwelled and Downwelled Irradiance, and Spectral Up- 
welled Radiance 

MER 21 Deck Radiometer 1 b 

Surface Incident Irradiance 

AC-9 2 

Spectral absorption and beam attenuation coefficients 

LS-6000 Light Scattering Sensor 3 

Light scattering at 880 nm 

Sea-Bird SBE 19 CTD 4 

Conductivity, temperature, and depth 

PSA-900 Sonar Altimeter" 

Height (m) of package above bottom 

Night Package 


AC-9 2 (x2) 

Spectral absorption, beam attenuation, backscattering coefficients 

LS-6000 Light Scattering Sensor 1 

Light scattering at 880 nm 

Sea-Bird SBE 19 CTD 4 

Conductivity, temperature, and depth 

PSA-900 Sonar Altimeter 5 

Height (m) of package above bottom 

Remote - Tethered Package 


SeaWiFS Profiling Multichannel Radi- 
ometer 6 

Spectral Upwelled Radiance, Spectral Downwelled Irradiance 

SeaWiFS Multichannel Surface Refer- 
ence 6 

Downwelled Spectral Surface Irradiance 


'Biospherical Instruments, Inc.; lb Biospherical Instruments, Inc., MER 21 is a deck-mounted sensor, 
2 WET Labs, Inc. (Western Environmental Technology Laboratories, Inc); 3 SeaTech, Inc.; 

4 Sea-Bird Electronics, Inc.; 5 Datasonics, Inc.; 6 Satlantic, Inc. 
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Chapter 1 7 

Bio-Optical Measurement and Modeling 
of the California Current and Polar Oceans 

B. Greg Mitchell and Piotr Flatau 
Scripps Institution of Oceanography 
University of California San Diego , La Jolla , California 


17.1 INTRODUCTION 

This SIMBIOS project contract supports in situ 
oceanic optical observations in the California Current 
and Southern Ocean The principal objectives of this 
research are to validate standard or experimental 
products through detailed bio-optical and biogeo - 
chemical measurements, and to combine ocean optical 
observations with advanced radiative transfer model- 
ing to contribute to satellite vicarious radiometric 
calibration and algorithm development. 

Our sampling efforts have been directed towards 
obtaining measurements in both the California Current 
and Antarctic polar waters, with the purpose of gener- 
ating a high-quality, methodologically consistent data 
set encompassing a wide-range of oceanic conditions. 
The combined data base includes stations which cover 
the clearest oligotrophic waters to highly eutrophic 
blooms and red-tides, and provides a coherent set of 
observations to validate bio-optical algorithms for 
pigments and primary production. This unique and 
comprehensive data is utilized for development of 
experimental algorithms (e g. high-low latitude pig- 
ment transition; phytoplankton absorption, photosyn- 
thesis, cDOM. 

17.2 RESEARCH ACTIVITIES 

The Southern California Bight region, from San 
Diego to just north of Point Conception, has one of the 
longest, most comprehensive time-series of marine 
observations; the California Cooperative Oceanic 
Fisheries Investigation (CalCOFI). This region expe- 
riences a large dynamic range of coastal and open 
ocean trophic structure, and has been extensively 
studied with respect to its regional optical properties 
in an effort to develop regional ocean color algorithms 
(e.g. Smith and Baker 1978, Mitchell and Kiefer 1988, 
Sosik and Mitchell 1995). During the first year of our 
contract, we participated in 4 quarterly cruises in the 
California Current region as part of the CalCOFI pro- 
gram. 

The Southern Ocean is a large, remote region 
which plays a major role in global biogeochemical 


cycling. Despite evidence that bio-optical relation- 
ships in these waters can diverge significantly from 
lower-latitude waters (e.g. Mitchell and Holm-Hansen 
1991), Antarctic waters have not been represented in 
the databases (e.g. SeaBAM) used to formulate and 
test modem ocean color algorithms. During the past 
year, we participated in 3 cruises to the Southern 
Ocean as part of the US JGOFS program. One cruise 
was located within the Ross Sea Polyna during the 
annual spring phytoplankton bloom, with two subse- 
quent cruises covering the region of Antarctic Polar 
Front Zone along 170° W. 

On all cruises, an integrated underwater profiling 
system was used to collect optical data and to charac- 
terize the water column. The system included an un- 
derwater radiometer (Biospherical Instruments MER- 
2040 or MER-2048) measuring depth, downwelling 
spectral irradiance (E d ) and upwelling radiance (L u ) in 
13 spectral bands. A MER-2041 deck-mounted refer- 
ence radiometer (Biospherical Instruments Inc) pro- 
vided simultaneous measurements of above-surface 
downwelling irradiance. Details of the profiling pro- 
cedure, characterization and calibration of the radi- 
ometers, data processing and quality control are de- 
scribed in Mitchell and Kahru (1998). The underwa- 
ter radiometer was also interfaced with 25 cm trans- 
missometers (SeaTech or WetLabs), a fluorometer, 
and SeaBird conductivity and temperature probes. 
When available, additional instrumentation integrated 
onto the profiling package included AC9 absorption 
and attenuation meters (Wetlabs Inc ), and a Hydro- 
scat-6 backscattering meter (HobiLabs). 

In conjuction with in situ optical measurements, 
discrete water samples were collected from a CTD- 
Rosette immediately before or after each profile for 
additional optical and biogeochemical analyses. Pig- 
ment concentrations were determined flurometrically 
and with HPLC. Spectral absorption coefficients 
(300-800 nm) of particulate material were estimated 
by scanning particles concentrated onto Whatman 
GF/F (Mitchell 1990) in a dual -beam spectropho- 
tometer (Varian Cary 1). Absorption of soluble mate- 
rial was measured in 10cm cuvettes after filtering 
seawater samples through 0.2 um pore size polycar- 
bonate filters. 
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17.3 RESEARCH RESULTS 

Validation of satellite-retrieved normalized water- 
leaving radiances (Lwn ) was done by comparing 
SeaWiFS images with in situ data collected concur- 
rently (± 4 hours). Satellite data received at the Mon- 
terey Bay Research Institute, the University of Cali- 
fornia Santa Barbara (CalCOFI region) and McMurdo 
Station, Antarctica (Southern Ocean region) were 
processed to Lwn using SeaDAS 3.2 software (Fu et 
al. 1998). A total of 16 matching sets of Lwn were 
found between 2-Oct-1997 and 21-Apr-1998 for the 
CalCOFI region. Because of persistent cloud cover in 
the Southern Ocean, only 3 matchups were possible in 
the Ross Sea region (all on 1 -Dec- 1997). Satellite 
values were derived as averages over 3x3 pixel areas 
centered at the in situ measurement. 

In both regions, these comparisons reveal signifi- 
cant under-estimation of the SeaWiFS-retrieved Lwn 
compared to in situ measurements. The differences 
were generally smallest in the 555 nm band, and larg- 
est at shorter wavelengths. The magnitude of under- 
estimation in the shorter wavlength bands increases at 
high Chi concentration. 

O’Reilly et al. (1998) describe the Ocean Color 2 
(OC2) chlorophyll algorithm that is used by NASA in 
the operational processing of SeaWiFS data (Fu et al. 
1998). This algorithm uses the ratio of remote sensing 
reflectances (R„ ) at 490 and 555 nm to estimate chlo- 
rophyll a concentration, with the coefficients derived 
by a statistical fit to a data set of 919 bio-optical 
measurements comprising the SeaBAM data set. 
More recently (August 1998), NASA announced a 
revised version of the OC2 (OC2-v2) which was in- 
tended to reduce the drastic over-estimation of Chi in 
high biomass waters produced by the original OC2 
algorithm. 


Figure 1 compares the performance of the OC2- 
v2 algorithm with our present data base of measure- 
ments from CalCOFI and the Southern Ocean When 
applied to the CalCOFI data, this algorithm overesti- 
mates chi a at very high chi a and underestimates 
elsewhere (with the exception of the extreme low chi 
a). A similar pattern is seen with the Southern Ocean 
data, although in general the degree of underestima- 
tion is greater and the transition to overestimation 
occurs at lower Chi. These results underscore the 
eventual need for specific regional empirical algo- 
rithms to obtain more accuate estimates of Chi and 
primary production from ocean color remote sensing. 
We have recently developed an improved empirical 
chlorophyll algorithm for the California Current 
(CAL-P6), which utilizes a sixth order polynomial of 
the ratio of Lwn at 490 and 555nm (Kahru and 
Mitchell, 1999). 

17.4 WORK PLAN 

Our participation in the quarterly CalCOFI 
cruises in the California Current will continue 
throughout the next period. We are also initiating 
field programs in the Indian Ocean (INDOEX) and the 
Sea of Japan (JES) to increase the regional scope of 
our data base. Modeling efforts include the pursuit of 
regional bio-optical algorithms for in water optical 
properties and their relationship to biogeochemical 
parameters, as well as the development of semi- 
analytical models for the retrieval of inherent optical 
properties from satellite data. We anticipate that these 
efforts will lead to an improved understanding of the 
variability observed in empirical satellite algorithms. 
We will also initiate analyses to determine the ele- 
ments of the SeaWiFS processing that lead to the un- 
derestimates of which is of particular concern 
for high chlorophyll waters. 
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Figure 1. A comparison of the OC2-v2 algorithm (solid line in left panels) with in situ measurements of 1 

chlorophyll a from CalCOFI and the Southern Ocean. The right panels illustrate quantile-quantile plots of the ! 

differences between modeled and measured chloroDhvIl a. 
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Chapter 18 

SIMBIOS Normalized Water Leaving Radiance 
Calibration and Validation: Sensor Responses, 
Atmospheric Corrections, Stray Light and Sun Glint 

James L. Mueller 

Center for Hydro-Optics and Remote Sensing 
San Diego State University, San Diego California 


18.1 INTRODUCTION 

The objectives of research under this contract are 
to improve the uncertainty budgets and protocols for 
normalized water leaving radiances Lwn(^) and re- 
mote sensing reflectance R^nC^-) determined from in 
situ measurements, and to apply in situ reflectance 
measurements to validate Lwn(^) and RonM deter- 
mined from radiances measured by the SIMBIOS en- 
semble of orbiting ocean color sensors (i.e. SeaWiFS, 
MODIS, GLI, POLDER and others). Our approach to 
these objectives includes protocol experiments to 
evaluate uncertainties in reflectances derived from 
above-water and in-water radiometry, field expedi- 
tions to acquire in situ radiometric, optical and bio- 
logical ground truth data for validating SIMBIOS sat- 
ellite ocean color sensors and match-up analyses, 
combining the in situ data with satellite radiometric 
data and derived parameters for vicarious radiometric 
calibration, and for evaluating sun glint effects and 
stray light artifacts in the image data. 

18.2 RESEARCH ACTIVITIES 

There is a general consensus that poorly under- 
stood optical variability associated with surface waves 
contributes strongly to the persistently significant dis- 
crepancies between remote sensing reflectances de- 
termined from above-water and underwater radiomet- 
ric measurements. In above-water radiance meas- 
urements, wave roughness modifies the surface re- 
flectance relative to a flat surface. Wave roughening 
also widens the effective field-of-view of reflected 
skylight components of the measured signal, and thus 
may incorporate angular variability in skylight which 
is not quantified by measurements using current pro- 
tocols (Mueller and Austin, 1995). Finally, wave 
slope variations repolarize reflected skylight in ways 
which may not be properly accounted for using these 
and proposed alternative protocols. 

In underwater radiometric profiles near the sea 
surface, surface waves radomly focus and defocus 


incident radiance and introduce “noise” which in- 
creases exponentially with proximity to the surface. 
Additional uncertainties result from strong wave- 
induced accelerations on an underwater radiometric 
package, which accelerations also increase exponen- 
tially with proximity to the surface and introduce 
quasi-random, difficult to measure fluctuations in 
sensor orientation. It is very difficult to separate these 
wave-related sources of uncertainty in comparing re- 
flectances calculated from in-water and above-water 
measurements. It is correspondingly difficult, there- 
fore, to partition the overall observed uncertainty 
between the two types of measurement. We are 
working towards comparative outdoor tank and field 
measurements which attempt to separately quantify 
near-surface in-water and above-water wave related 
uncertainties under semi-controlled conditions. 

To support the difficult measurements associated 
with these protocol experiments, we are developing a 
unique hyperspectral above-water and profiling spec- 
troradiometer and inherent optical property (IOP) 
system. The system is designed to allow repeated 
radiometric and IOP profiling at controlled ascent and 
descent rates, concurrent with measurements of 
above-water irradiance and radiance. The profiler 
package is built around a tethered variable buoyancy 
Autonomous Profiling Vehicle (APV) manufactured 
by Ocean Sensors Inc. of San Diego. In addition to 
the APV’s built-in Conductivity, Temperature and 
Depth (CTD) sensors, sensors are attached to the pro- 
filer package to measure irradiance E d (\,z), radiance 
L u (?.,z), absorption a(X.,z), beam attenuation c(/.,z) 
and chlorophyll a fluorescence. E d (L,z) and L u (a,z) 
are measured at 3.3 nm intervals from approximately 
380 to 800 nm using radiometers based on Zeiss 
miniature fiber-optic spectrometers. Chlorophyll a 
fluorescence is measured using a WETSTAR fluo- 
rometer, and a(X,z) and c(X,z) are measured at 3.3 nm 
intervals using a HISTAR underwater spectropho- 
tometer - both instruments are manufactured by 
WETLABS of Philomath, Oregon. (The WETSTAR 
and HISTAR are obtained under a ONR/ECOHAB 
research grant at CHORS.) Additional radiometers 
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are being integrated into the system for above-water 
measurements of irradiance E s (X,z), total water- 
leaving radiance L sfc (X,z) and sky radiance L sky (X,z). 

Validation Field Experiments 

We are acquiring in situ data for validating 
SeaWiFS (and in the future MODIS and other ocean 
color sensors) data during field expeditions aboard the 
Mexican Research Vessel Francisco de Ulloa (FdU) 
in the Gulf of California (GoCAL). The GoCAL 
cruises are in collaboration with Drs. Ron Zaneveld 
and Scott Pegau of Oregon State University (SIM- 
BIOS contract NAS5-97129) and Mexican scientists 
(Drs. Helmut Maske, Ruben Lara-Lara, Saul Alvarez- 
Borrego, and their students) at CICESE (an oceanog- 
raphy PhD granting postgraduate university in En- 
senada, BC, Mexico). 

The FdU is owned and operated by CICESE, and 
they have thus far contributed ship time for cruises in 
October 1997 (GoCal-97), March 1998 (GoCal-98A) 
and Nov-Dec 1998 (GoCal-98B). Measurements on 
each cruise include profiles of E d (X,z) and z) at 
SeaWiFS wavelengths, fluorometric chlorophyll a, 
and HPLC pigments. IOP measurements are made by 
the OSU investigators, and additional bio-optical and 
ancillary measurements are contributed by CICESE. 
The next GoCal cruise aboard the FdU (GoCal-99A) 
is tentatively scheduled for 21 April - 2 May 1999. 

All radiometric, optical and phytoplankton pig- 
ment measurement and data analysis procedures fol- 
low protocols in Mueller and Austin (1995). Radio- 
metric profiles are analysed to obtain E^O*), 
Ea(^),0 + ), and L u (k ? 0*), using the integral method of 
Mueller (1995), as implemented using the CHORS 
interactive Kfit routine. 

We are also planning to acquire validation data 
during day-cruises near San Diego and other Califor- 
nia coastal sites using a small research boat. We have 
recently ordered and will outfit a trailerable research 
boat (a Parker-25 sport-fishing boat) to support this 
and other research projects at CHORS. Electronic 
equipment will be installed in the boat’s enclosed pi- 
lothouse and forward cabin, while water sam- 
pling/filtering stations will be installed aft of the pi- 
lothouse under a permanent awning. Costs (acquisi- 
tion, maintenance, and operating) associated with the 
CHORS research boat will be shared between this 
contract, other CHORS research projects using the 
boat, and SDSUF cost-sharing funds. 

Calibration and Characterization 

To support these experiments and analyses, 
CHORS carries out calibrations and characterizations 
of radiometers used in this project, including those of 
our collaborators at OSU and CICESE. NIST trace- 
ability of CHORS calibrations is maintained through 


tertiary standards (Optronics, Inc) and assured 
through participation in SIRREX and independent 
comparisons with other laboratories calibrating SIM- 
BIOS instruments. 

18.3 RESEARCH RESULTS 

Radiometric data from GoCal-97 and GoCal-98A 
(including the results of the CHORS KFTT analysis of 
OSU’s SPMR profiles) have been processed, analysed 
and archived in SeaBASS. Pigment data analyses and 
quality control have recently been completed for these 
two cruises, and the data will be archived in SeaBASS 
after appropriate reformatting and insertion of meta- 
data. Processing and analysis of data from GoCal-98B 
(24 Nov through 8 Dec 1998) are in progress. 

18.4 WORK PLAN 

During the current contract year we will complete 
work in progress including: 

• data processing, analysis and SeaBASS archi- 
val/documentation of data from GoCal-98B, 

• Analysis and publication, in collaboration with 
investigators from CICESE and OSU, of data 
from completed GoCal cruises in 1995, 1996, 
1997, 1998 A and 1998B, (including SeaWiFS 
data for the 1997 and 1998 periods); 

• Integration, laboratory characterization and cali- 
bration, and field-testing of the CHORS hyper- 
spectral radiometric and IOP profiling and re- 
flectance system; 

• Continued routine calibrations of radiometers 
used by CHORS and our collaborators (OSU and 
CICESE) in this project. 

During the current contract year, planned field cam- 
paigns include: 

• GoCal-99A aboard the FdU in the Gulf of Cali- 
fornia between 20 April and 2 May 1999, and 

• Day cruise, clear-sky validation sorties (within 20 
nautical miles of the California coast) aboard the 
CHORS research boat ( circa June 1999). 

Reflectance protocol experiments will include 
tank and field in-water versus above-water reflectance 
measurement comparisons, as well as experimental 
evaluations of plaque reflectance versus direct meas- 
urements of incident irradiance and instrument self- 
shading tests. 
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Chapter 19 


Validation of Carbon Flux and Related Products for 
SIMBIOS: the CARIACO Continental Margin 

Time Series 

Frank Muller-Karger 
Department of Marine Science 
University of South Florida 

Ramon Varela 

Fundacion La Salle de Ciencias Naturales 
Punta de Piedras , Estado Nueva Esparta, Venezuela 


19.1 INTRODUCTION 

This SIMBIOS investigation focuses on validating 
ocean color satellite products from waters near conti- 
nental margins. Specifically, our program focuses on 
collecting monthly observations within a coastal up- 
welling site, and seasonal extreme measurements 
within the plume of the Orinoco River, which is classi- 
fied as having the third largest discharge in the world. 

Initial focus of this project has been on setting up 
the logistics for field data collection and initiating 
routine observations. The program incorporates an 
effort to validate SeaWiFS products and products from 
the future MODIS (AM and PM), MISR, MERIS, 
GLI, and ROCSAT missions. Our goals are to evalu- 
ate algorithms used to compute biological and geo- 
physical data products, develop new products aimed at 
assessing carbon flux, estimate the spatial and tempo- 
ral errors in these products, and identify strategies that 
will help minimize such errors and eliminate possible 
biases between data sets from different missions. 

19.2 RESEARCH ACTIVITIES 

The core of this SIMBIOS project is an estab- 
lished oceanographic time series stations in the 
CARIOCO basin. In November 1995, the U S. Na- 
tional Science Foundation, in collaboration with the 
Consejo Nacional de Investigaciones Cientificas y 
Tecnologicas (CONICIT) of Venezuela, implemented 
a comprehensive program of continuous observations 
in the Cariaco Basin, which is located in the southeast- 
ern Caribbean Sea. This multidisciplinary program is 
referred to as the CArbon Retention In A Colored 
Ocean (CARIACO) program. Specifically, the pro- 
gram conducts monthly observations at a station lo- 
cated at 10.5 N, 64.67 W. 


The importance of the Cariaco Basin as a natural 
laboratory to study seasonal upwelling, sediments and 
anoxic effects has been appreciated for over 40 years. 
The CARIACO Program seeks to define a budget 
which assesses the total C02 upwelled with deep, 
nutrient-rich water with respect to the annual export of 
organic carbon from surface waters along this conti- 
nental margin. This basin is ideal for a carbon flux 
study because it is a large depression on a continental 
shelf which experiences vigorous wind-driven up- 
welling, and oligotrophic conditions at different times 
of the year. As the circulation below 200 m in this 
basin is restricted, it forms a natural sediment trap 
The Cariaco Basin is the only permanently anoxic 
oceanic basin, and as such serves as an oceanic analog 
of the Black Sea.The series station is located at 10.50 
N, 64.66 W and consists of: 

• 1 to 3 monthly cruises with a fully-equipped, 
modem oceanographic vessel; 

• Mooring with 4 sediment traps (200, 400, 800, 
1200 m; capturing bi-weekly sample integra- 
tions); 

• Varved (laminated) sediment core analyses. 

The monthly oceanographic cruises provide in- 
formation on: 

• Composition, levels, and light absorption proper- 
ties of organic particulate and dissolved matter; 

• Pigments, HPLC, taxonomy of phytoplankton and 
general classification of bacteria; 

• Biological productivity (phytopl. and bacteria); 

• The physical/chemical characteristics of the water, 
including nutrient and oxygen levels; 

• Carbonate system. 
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Shore-based, continuous meteorological and tide 
gauge observations are available, and we are trying to 
obtain these in collaboration with the local Air Force, 
Navy, and Environment Ministries. This project incor- 
porates subsurface and above-surface spectral radi- 
ance, irradiance, and aerosol optical thickness meas- 
urements to the CARIACO series. The in situ bio- 
optical observations are all done following the 
SeaWiFS protocols. We have also participated in 
other SIMBIOS activities, including the first SIM- 
BIOS PI meeting in Salomon's Island and a cal/val 
cruise to the Bahamas organized under Ken Carder's 
SIMBIOS program. 

19.3 RESEARCH RESULTS 

Each month since December 1995, bio-optical 
measurements have been part of the suite of measure- 
ments made at the CARIACO site. These observations 
were initially limited to above-water hyperspectral 
reflectance measurements. MER casts were performed 
only a few months out of this period because there was 
no continuous access to such an instrument. 

Since October 1997, our suite of observations has 
included MER subsurface radiance/irradiance profiles. 
The optical measurements are used to validate 
SeaWiFS. The measurements are scheduled around 
the time of the SeaWiFS overflight. The monthly 
cruises last about 2 days each. 

During each cruise, we collect bio-optical meas- 
urements using the underwater profiling Biospherical 
Instruments MER2048 in conjunction with a 
MER2041 Deck Cell. Above water measurements are 
made with a Photo Research Hyperspectral Color- 
imeter Model PR650, which is a calibrated, hand-held 
radiometer (8 nm resolution). Derived products in- 
clude Lw (Water-Leaving Radiance), Rrs (Remote- 
sensed reflectance) and K (attenuation coefficient). In 
addition, a full suite of measurements is made which 
includes: particulate material and pigment absorption, 
HPLC, fluorometric determinations of Chi concentra- 
tion, pH, Alkalinity, Primary Productivity, POC, DOC 
absorption and concentration, nutrients, sun photome- 
try, oxygen, and salinity. 

Cruise reports for each cruise starting in Septem- 
ber 1997, have been submitted to the SIMBIOS Proj- 
ect. Approximately two years of bio-optical data have 
been submitted to the SeaBASS archive. 

Preliminary results suggest marked seasonal 
changes in optical properties occur in this basin as 
phytoplankton concentration increases and decreases 
in response to coastal upwelling. Figure 1 shows the 
seasonal variation in surface hyperspectral remote 
sensing reflectance at the Cariaco station. A clear pro- 
gression is seen from clear water conditions during 
boreal summer to blue-poor reflectances during boreal 
winter-spring, the time of upwelling. The decrease in 
blue reflectance is accompanied by a marked increase 


at 683 nm, clearly defining the signal of solar- 
stimulated fluorescence. We plan to examine this fea- 
ture in relation to the primary productivity data col- 
lected concurrently during our CARIACO cruises. 

We also find that there is a close correlation be- 
tween phytoplankton specific absorption coefficients, 
primary productivity, and chlorophyll a. These data 
clearly are preliminary and need to be examined closely 
to understand variations from season to season. 

The Orinoco River Plume 

Above-water spectral remote-sensing reflectance 
(Rrs 380-700 nm) was measured including corrections 
for sky-radiance. We also conducted underwater 
measurements of remote sensing reflectance at 7 
wavelengths (412, 443, 455, 475, 490, 510, 532, 560, 
589, 625, 671, 683, and 700 nm) using submersible 
MER2048 instrumentation (Biospherical Instr ). Some 
observations were conducted away from the coast, 
where the color of the plume is influenced by colored 
dissolved organic matter as well as phytoplankton, but 
where suspended terrestrial material is expected to be 
minimal. We also collected data very close to the 
coast where the plume is extremely turbid. As of Feb- 
ruary 1, 1999, we are currently submitting bio-optical 
data collected during the first and second 1998 Ori- 
noco plume cruises to SeaBASS 

The Orinoco Plume data span a very broad range of 
remote sensing reflectance values. The river plume 
data are also being used by Joe Salisbury and Charles 
Vorosmarty at the University of New Hampshire for 
automated modeling of river water impact off the conti- 
nent of South America, in a model linked to terrestrial 
hydrology. 

19.4 WORK PLAN 

We will continue to occupy the CARIACO sta- 
tion on a monthly basis. In addition, we plan to con- 
duct two cruises to the Orinoco River plume, specifi- 
cally in February and October 191)9. Analysis of data 
and samples will continue. We are currently finalizing 
a manuscript that describes the seasonal cycle in hy- 
drographic properties, phytoplankton biomass and 
productivity, and vertical carbon flux at the 
CARIACO station. We will now initiate work on a 
manuscript that summarizes the bio-optical data from 
this location. 
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Figure 1 . Remote sensing reflectance at the CARIACO station over the course of one year. 
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Chapter 20 


Bermuda Bio-Optics Program (BBOP) 

David Siegel 

University of California at Santa Barbara , 

Santa Barbara, California 


20.1 INTRODUCTION 

The Bermuda BioOptics Project (BBOP) is a col- 
laborative effort between the Institute for Computa- 
tional Earth System Science (ICESS) at the University 
of California at Santa Barbara (UCSB) and the Ber- 
muda Biological Station for Research (BBSR, D. A. 
Siegel [UCSB] and N. B. Nelson [BBSR] P i 's). This 
research program is designed to characterize light 
availability and utilization in the Sargasso Sea and in 
doing so, provide an optical link by which biogeo- 
chemical observations may be used to evaluate bio- 
optical models for pigment concentration, primary 
production, and sinking particle fluxes from satellite- 
based ocean color sensors. 

20.2 RESEARCH ACTIVITIES 

The Bermuda Bio-Optics Program (BBOP) col- 
lects detailed profiles of AOPs and IOPs in conjunc- 
tion with the US JGOFS Bermuda Atlantic Time- 

series Study (BATS) at 31 50* N, 64 10’ W in the 
mesotrophic waters of the Sargasso Sea. There are up 
to 16 cruises per year, conducted monthly with addi- 
tional cruises during the spring bloom period, January 
through May. 

Continuous profiles of apparent optical properties 
(AOPs) are collected in the upper 200m using a Mul- 
tispectral Environmental Radiometer (MER- 
2040/2041, Biospherical Instruments Inc., San Diego 
CA; Smith et al., 1984). The primary optical meas- 
urements are downwelling vector irradiance and up- 
welling radiance, Ed(z,t,X) and L u (z,t, X), respec- 
tively. The instrument samples these properties at 12 
wavebands (including SeaWiFS wavelengths) plus 
broadband natural fluorescence. The BBOP sampling 
package also includes sensors for temperature and 
conductivity (SeaBird, Bellevue WA) chlorophyll 
fluorescence and red beam transmission (SeaTech, 
Corvallis OR) and a second mast-mounted spectrora- 
diometer (MER-2041) with wavebands similar to 
those on the underwater instrument for measuring 
incident downwelling vector irradiance, Ed(0 + ,t, X). 
Instrument deployments are planned to optimize 


match-ups with the BATS primary production incuba- 
tions and, when possible, with SeaWiFS overpasses. 

The BBOP group also collects several measure- 
ments of inherent optical properties (IOPs). Continu- 
ous profiles of spectral beam absorption and attenua- 
tion, a(z, X) and c(z, X), are measured using the WET 
Labs AC-9 (Brody, 1998), and discrete samples for 
determining the absorption spectra of particulates, 
aph(z, X) and ad(z, X), and CDOM (ag(z, X)) are col- 
lected according to Nelson et al (1998). Spectral 
beam absorption measurements have been hampered 
by the fact that the baseline absorption of Sargasso 
Sea water is nearly the same as that of the deionized 
water used for the calibration. Particulate absorption 
spectra are determined using the quantitative filter 
technique (Mitchell, 1990) and CDOM absorption 
according to Nelson et al (1998). In 1997, we initi- 
ated the collection of above water R^(l) spectra using 
the Analytical Spectral Devices FieldSpec spec- 
trometer (ASD, Boulder Colorado). We are currently 

comparing this data set with values of Rrs(l) calcu- 
lated from MER data to evaluate its usefulness for 
SIMBIOS. 

The profiling radiometers are calibrated three 
times per year at UCSB and have demonstrated re- 
markable calibration stability for several years, al- 
lowing the use of long-term averages for calibration 
coefficients (O’Brien et al., 1999). This stability, 
combined with streamlined procedures for data analy- 
sis, has made it possible for us to report many of our 
derived products with -98% reliability in near-real 
time, including profiles of remote sensing reflectance 
(R rs (z, k) = L u (z, X)/Ed(z, X)), and down- and up- 
welled attenuation coefficients (Kd(z, X), K u (z, X)). 
Our analysis procedures are fully documented in Sie- 
gel et al. (1995b). In addition to radiometer data, 
samples for fluorometric chlorophyll-a are collected 
before dawn and near noon (local time) and results are 
delivered to the SIMBIOS project with the AOP data, 
usually within one week after a cruise. Table 1 con- 
tains a list of data products collected by BBOP and/or 
BATS, which are relevant to SIMBIOS. A total of 
472 AOP, 16 IOP and 36 pigment profiles were deliv- 
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ered from the start of this contract (July 1997) through 
1998. 

20.3 RESEARCH RESULTS 

Since its inception in 1992, BBOP has made a 
number of important contributions toward the scien- 
tific understanding of the relationship between light 
and ocean biogeochemistry. The Sargasso Sea near 
Bermuda is mesotrophic, characterized by both eutro- 
phic and oligotrophic conditions during different times 
of the year. Although low chlorophyll a stocks and 
primary production rates prevail for the most of the 
year, there is a short spring bloom characterized by 
somewhat higher concentrations and rates which is 
controlled by winter mixing and the associated nutri- 
ent influx. 

Our research has shown that the BATS/BBOP site 
has an interesting and important ocean color signal 
associated with its colored dissolved organic material 
(CDOM) content. We have found that CDOM may 


comprise up to 50% of the light attenuation budget 
(Kd(A.)) at 410 nm (Siegel et al. 1995a, Siegel and 
Michaels 19%) and that CDOM concentrations can be 
determined from in situ ocean color spectra (Garver 
and Siegel, 1997). The ratio of Kd(410) to K<i(488) 
can be used as an index of CDOM concentration to 
illustrate its time-depth dynamics (Siegel et al. 1995a, 
Siegel and Michaels 1996). During the summer, 
CDOM is simultaneously produced within the sea- 
sonal pycnocline and photo-oxidized in the surface 
mixed layer (Figure lc). Fall and winter mixing ho- 
mogenizes the profile. The seasonal distribution of 
CDOM is unrelated to the those of chlorophyll or dis- 
solved organic carbon, however it may be that the 
summer CDOM source is microbial (Siegel and 
Michaels 1996, Nelson et al. 1998). Since CDOM is a 
significant source of light attenuation in the Sargasso 
Sea that varies independently of chlorophyll, it ap- 
pears that the optical definition of Case I waters which 
is based solely on chlorophyll should be reconsidered 


Table 1. A partial list of measurements made by BBOP and BATS 


Direct Measurements: 

Ed(z, A) Downwelling vector irradiance (410,441,465,488,510,520,555,565,589,625,665 & 683 nm) 

Ed(0 + >?l) Incident irradiance (340,390,410,441,465,488,520,545,565,589,625,665 & 683 & 350-1050 nm) 

L u (z, X ) Upwelling radiance (410,441,465,488,510,520,555,565,589,625,665 & 683 nm) 

a(z, X) In situ absorption spectrum using WetLabs AC-9 (410,440,490,520, 565, 650, 676&715 nm) 

c(z, A.) In situ beam attenuation spectrum (same A's as above) 

atp(A) Particulate absorption spectrum by QFT 

a<j(A) Detrital particle absorption spectrum by MeOH extraction 

a ys (A) Colored dissolved absorption spectrum 

E 0 (z, A.) Scalar irradiance at 441 and 488 nm 

Ff(z) Natural chlorophyll fluorescence using a broadband upwelled radiance sensor 

chl-fl(z) Chlorophyll fluorescence with a SeaTech fluorometer 
c(z,660) Beam attenuation coefficient at 660 nm with SeaTech 25 cm transmissometer 
T(z) & S(z) Temperature and conductivity with SeaBird probes 

chl-a(z) Discrete chlorophyll a determinations via Turner fluorometry (for next day delivery) 


Primary Derived Products: 

L wN (A) Normalized water leaving radiance (410,441,465,488,510,520,555,565,589,625,665 & 683 nm) 
Rr S (0-,X) In-water remote sensing reflectance (410,441,465,488,510,520,555,565,589,625,665 & 683nm) 
Rr8(0 + ,/„) Above- water remote sensing reflectance (350 to 1050 nm) 

a p h(A) Phytoplankton absorption spectrum (= a p (A) - adet(^)) 

Kd(z, A.) Attenuation coefficient for Ej(z, X) (410,441,465,488,510,520,555,565,589,625,665 & 683 nm) 

Kl(z, A) Attenuation coefficient for L u (z, A) (410,441,465,488,510,520,555,565,589,625,665 & 683 nm) 

<PAR(z)> Daily mean photosynthetically available radiation at depths of the in situ C 14 incubations 
b(z. A) Spectral scattering coefficient (= c(z, A) - a(z, A)) 
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Chapter 21 

SIMBIOS Data Product and Algorithm Validation 
with Emphasis on the Biogeochemical 
and Inherent Optical Properties 

J. Ronald V. Zaneveld and Andrew H. Barnard 
Oregon State University 
Corvallis \ Oregon 


21.1 INTRODUCTION 

The purpose of the SIMBIOS program is to pro- 
vide routine, long-term ocean color data from the 
various international ocean color missions. Our 
funded component of this program is to address and 
quantify the relative accuracies of the products from 
the international missions by means of product and 
algorithm validation. Our long-term objectives are to 
a) collect optical and biochemical data in oceanic and 
coastal regions including a cooperative international 
bio-optical monitoring program in the Gulf of Califor- 
nia; b) provide collected data to the NASA database, 
SIMBIOS project office, and other interested users; 
and c) determine spatial and temporal error fields for 
the biological and geophysical data products from the 
various ocean color missions. 

21.1 RESEARCH ACTIVITIES 

Over the past 4 years we have been involved in a 
multi-national research effort focused on determining 
the physical and biogeochemical variability in the 
Gulf of California. This research is done in coordina- 
tion with Drs. J. Mueller and C. Trees from San Diego 
State University (SDSU), Drs. S Alvarez Borrego, R. 
Lara-Lara, G. Gaxiola and H. Maske from the Centro 
de Investigacion Cientifica y de Educacion Superor de 
Ensenada (CICESE), Ensenada, Mexico, and Dr. E. 
Valdez from the University of Sonora, Hermosillo, 
Mexico. Our component has been to measure the in- 
herent and apparent optical properties as well as the 
traditional CTD parameters on vertical scales less than 
0.5 m. We have completed 5 cruises in the Gulf of 
California, with 3 of these cruises occurring in the last 
year and half during our SIMBIOS contract period. 
Optical data collected during these cruises are cur- 
rently being used to validate SeaWiFS algorithms. 
We have also completed a cruise of opportunity off- 
shore of Oregon in September 1998. Data from this 
cruise is currently being processed and analyzed. 


Two sampling platforms are typically used during 
these cruises; a SLOW Descent Rate Optical Profiler 
(SLOWDROP) to measure the inherent optical prop- 
erties as well as the physical parameters, and a profil- 
ing radiometer system to measure the downwelling 
irradiance and upwelling radiance. Typically, 
SLOWDROP and radiometer profiles are made within 
an hour of each other at the same location 

The SLOWDROP platform is free-falling and 
slightly negatively buoyant and provides 10-20 cm 
scale vertical resolution. A typical instrument con- 
figuration includes a SeaBird SBE-25 CTD for meas- 
uring physical parameters, a fluorometer, and two 
WETLabs ac-9 (9 wavelength absorption and attenua- 
tion meters) for measuring particulate and dissolved 
component of the inherent optical properties (IOP). 
CTD and IOP data are collected simultaneously dur- 
ing a single profile using a Wet Labs, Inc MOD APS 
system for data collection and integration. Calibra- 
tions of the ac-9 instruments were performed several 
times during the contract period. Field calibrations of 
the ac-9 meters were typically done once per sampling 
day using a clean water standard produced using a 
Bamstead nanopure water system according to the 
SIMBIOS protocols for this instrument. 

A Satlantic Inc. SPMR system is used to measure 
profiles of the apparent optical properties (AOP) in- 
cluding downwelling irradiance and upwelling radi- 
ance at 412, 444, 490, 533, 555, 590, and 684 run 
The Satlantic radiometer is calibrated at SDSU by 
Drs. Mueller and Pegau following NIST protocols and 
calibrated twice a year by Satlantic Inc. The optical 
filters on this radiometer were upgraded by the manu- 
facturer (Satlantic) over the past year. 

One meter binned profiles of all of the data, are 
provided to the SeaBASS database within 90 days of 
completion of each cruise. Also provided are the cali- 
bration histories of the ac-9's and the radiometer, and a 
detailed logbook. In the future we plan to make these 
data sets available via the web. 

The dates of each cruise and the number of pro- 
files made using each system are listed in Table 1. 
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Cruises in 1997 and 1998 were SIMBIOS supported. 
Figure 1 shows the station locations sampled during 
each cruise. Data collected during the fall cruise in 
1998 is currently being processed and quality con- 
trolled. As is such, the locations of the profiles are not 
shown on Figure 1. Currently, we are receiving all of 
the SeaWiFS HRPT Level 1 A data for the West Coast 
of North America including the Gulf of California 
from the DAAC. This data is being archived and 
processed to various levels for use in the validation of 
ocean color products. We have been contracted to 
maintain an optical profiling system for use in the 
SWIMBIOS instrument pool. The system contains 
two hyperspectral absorption and attenuation meters 
(WET Labs Histars), a CTD (SeaBird Electronics 25), 
a data acquisition system (WET Labs MODAPS+), 
cage, and necessary cabling. We have assembled the 
above system and have been testing the system in 
various local environments. Currently there are still 
some issues to be resolved as to the Hi Star measure- 
ments in comparison to the ac-9 measurements, which 
at this point seems to be calibration dependent. We 
are currently working closely with WET Labs to re- 
solve these issues. A full time technician has been 
hired to accompany this system, to assist users in de- 
ployment and operation and to process the data. We 
also maintain 3 deionizing water filtration systems 
(Bamstead Nanopure) which are for use by the SIM- 
BIOS authorized users. These systems are used to 
produce a clean water standard for ac-9 and HiStar 
field calibrations. These systems are currently in hand 
and each has been tested and is ready for general use. 
We expect that the calibration issues with the profiling 
package will be fully addressed and the system ready 
for general usage by the SIMBIOS authorized users by 
late February 1999. 

21.3 RESEARCH RESULTS 

Over the past contract period, we have begun 
analyzing our field data collected over the past 4 
years. Data collected in the Gulf of California as well 
as other locations were used to examine the spectral 
relationships between the inherent optical properties in 
Barnard et al. (1998). The results of this research 


were used to compare and contrast the inherent optical 
properties observed in the Gulf of California. In the 
shallow northwest region of the Gulf, Case II waters 
are consistently observed due to the resuspension of 
sediments. In the central portion of the Gulf, the larg- 
est variability in optical properties is associated with 
coastal upwelling. Away from the coasts the spectral 
optical properties are fairly constant. Validation of the 
normalized water leaving radiance measured by 
SeaWiFS in the Gulf of California during the 1997 
and spring 1998 cruises was examined. The results of 
this research indicated that atmospheric correction was 
problematic, often causing negative water leaving 
radiances at412nm. 

In coordination with our SeaWiFS funding, we 
have undertaken the development of a new remote 
sensing algorithm involving the triple ratio of the re- 
mote sensing reflectance. We have also compared the 
in situ radiance with the upwelling radiance resolved 
by the OCTS sensor off the northeastern U S. coast. 

21.4 WORKPLAN 

We will be participating in two research cruises in 
the Gulf of California during 1999, tentatively planned 
for the spring and fall. We anticipate participating in 
other cruises in the next year, as the opportunity 
arises. We will also be working on the installation of 
an in-line fluorometer on the trans-Gulf ferry to pro- 
vide a monthly transect of chlorophyll fluorescence 
and beam transmission data in the Gulf of California. 

In the next year, we anticipate that the hyper- 
spectral profiling package will receive a high amount 
of usage by the authorized users. We plan to maintain 
these systems as well as provide technical support to 
the users. Again we anticipate that this system will be 
available for general use by the end of February 1999 
We plan to continue to investigate the inherent and 
apparent optical properties of the Gulf of California, 
as well as continuing our work with ocean color algo- 
rithm validation and development. We will also be 
addressing issues of atmospheric correction over the 
Gulf of California, which have been found to be 
problematic. 


Table 1 Dates and number of profiles made during each of the Gulf of California cruises. 


Cruise Name 

Dates 

# IOP profiles 

# AOP profiles 

GOC95fall 

11/25/95 to 12/6/95 

28 

69 

GOC96fall 

10/30/96 to 11/7/96 

25 

25 

GOC97fall 

10/16/97 to 10/29/97 

39 

24 

GOC98spring 

3/6/98 to 3/16/98 

27 

7 

GOC98fall 

1 1/27/98 to 12/6/98 

45 

11 


1 


■ 

E 
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Figure 1. Map of the Gulf of California with sampling locations. The circles denote sampling locations 
during the 1995 cruise, the stars the 1996 cruise, the triangles the 1997 cruise, and the squares are for the spring 
1998 
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Chapter 22 

Validation of the SeaWiFS Atmospheric Correction 
Scheme Using Measurements of 
Aerosol Optical Properties 

Mark A. Miller 

Brookhctven National Laboratory 
Upton , New York 


22.1 INTRODUCTION 

The focus of the research described below is 
comparisons between surface and satellite measure- 
ments of aerosol optical properties as a medium for 
validating the assumptions made in the aerosol models 
used in the retrieval of water-leaving radiance. Two 
research projects are being conducted as part of the 
SIMBIOS contract: The first project is to develop a 
Marine Fast Rotating Shadow-band Spectral Radi- 
ometer (MFRSSR) that can be used for continuous, 
ship-board measurements of the spectral aerosol opti- 
cal depth and the spectral direct-normal and spectral 
diffuse irradiance. The second project is to measure 
and interpret aerosol properties in regions under sur- 
veillance by an ocean color satellite, thereby providing 
a data set from which to evaluate the atmospheric cor- 
rection algorithm. This second project entails the use 
of AERONET and ship-board photometers and radi- 
ometers. Progress toward the completion of each of 
these projects and work planned for the next period is 
described in the subsections below. 

22.2 RESEARCH ACTIVITIES 

MFRSSR 

There are two primary techniques used to meas- 
ure the aerosol optical depth: sun photometry and ro- 
tating shadow-band radiometry. Just as sun photome- 
ters must be accurately oriented toward the sun, an 
obvious pitfall at sea, conventional shadow-band radi- 
ometers have required exact orientation. The devel- 
opment of fast-rotating, shadow-band radiometers, a 
hybrid form of the original shadow-band technique, 
has removed orientation requirement, whereupon 
shipboard use of shadow-band radiometers has been 
encouraged. As compared to land-based units, ship- 
board fast-rotating shadow-band radiometers require 
much faster rotation of the occulting arm, fast- 
response silicon detectors, higher data sampling rates, 
and more sophisticated data analysis. As part of a 


multi-agency effort between the Department of En- 
ergy (DOE) and NASA SIMBIOS, the Brookhaven 
National Laboratory has developed an MFRSSR. It 
measures spectral global and spectral diffuse irradi- 
ance, which are used to calculate the spectral di- 
rect-normal irradiance as a function of solar zenith 
angle. From these data, the aerosol optical thickness 
can be computed during clear periods using the Lan- 
gley regression technique, or continuously if high 
quality extraterrestrial calibration coefficients are 
known. 

Initial data comparisons of optical depths from 
hand held sunphotometers have been encouraging and 
an example of a time series of optical depth from the 
MFRSSR and MicroTops is shown in Figure 1. In 
this case, a well-defined decrease in the optical depth 
during a 9-hour measurement period was observed by 
both instruments. Periods of noisy optical thickness 
are the result of clouds, which have not been filtered 
in this plot. The MFRSSR is designed to operate for 
long periods with minimal attention and it is sched- 
uled, along with two new units, to be at sea for a total 
of 7-8 months during 1999. During this period meas- 
urements will be taken over three oceans aboard the 
R/V Ron Brown and other ships, and from an island in 
the Tropical Western Pacific. At the conclusion of 
these cruises, all three units will spend several weeks 
at Mauna Loa for calibration. An important element 
of the MFRSSR design philosophy is automated op- 
eration and timely data analysis. Toward this end, 
efforts have been made to develop a cloud-filtering 
algorithm that can be used to identify cloud-free con- 
ditions. The initial version of this algorithm was 
based on the work of Long et al. (1999), but it has 
since been modified for more efficient use in a marine 
environment. 

AERONET SeaWiFS Comparisons 

Another approach to validating the atmospheric 
correction algorithms is to use land-based sunpho- 
tometers deployed in coastal locations, such as 
AERONET. These systems make a full range of 
measurements that can be used to compute the atmos- 
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pheric correction parameters defined by Gordon et al. 
(1980). It is not possible, with current instrumenta- 
tion, to measure the full range of variables necessary 
to compute the atmospheric correction parameters at 
sea. One drawback of the land-based approach, how- 
ever, is that near-shore measurements may not faith- 
fully represent conditions further from the coast, 
whereupon an approach that combines ship-board and 
AERONET measurements is probably the best alter- 
native. Both approaches are part of this study. 

The AERONET and SeaWiFS optical properties 
for some selected locations are being compared using 
the SeaWiFS/AERONET match-up code developed at 
NASA GSFC. The initial form of this code selects 
dates and times of coincident SeaWiFS and 
AERONET measurement and extracts the aerosol op- 
tical depth from the respective data files, providing the 
following are true: (1) AERONET data that have 
passed a cloud filter algorithm are available during the 
two-hour period centered on the satellite overpass 
time, and (2) the SeaWiFS data in each individual 
pixel in a scene pass all quality control checks. 

The match-up code identifies coincident SeaWiFs 
and AERONET measurement times at several loca- 
tions and detailed examinations of these data are in 
progress. There are several points to be made. At pre- 
sent, we have analyzed the matchup data for the fol- 
lowing locations: Dry Tortugas, Bermuda, San 
Nicholas Island, Lanai, Kaashidoo, and Bahrain. Pre- 
liminary results suggest that the pixel-to-pixel vari- 
ability in the aerosol optical thickness at most of the 
sites is unrealistically large, perhaps due to an inade- 
quate cloud filtering algorithm applied to the 
SeaWiFs data. The best match-up data sets analyzed 
to date are from Bermuda, Bahrain, and Kaashidoo. 
Comparing Bermuda and Kaashidoo, both island sites, 
the aerosol optical thickness as measured by both 
AERONET and satellite suggests that the aerosol op- 
tical thickness is generally higher at Kaashidoo, which 
lies in the Indian Ocean, and it tends to be more vari- 
able than at Bermuda. In general, more data tend to 
pass the SeaWiFs quality control algorithms at 
Kaashidoo, which may reflect less cloudiness more 
than any other single factor. The data suggests that 
the comparison between the SeaWiFS and AERONET 
is somewhat better at Bermuda than at Kaashidoo. 
The reason for the less favorable comparison data at 
Kaashidoo is under investigation, but some evidence 
suggests that the AERONET sunphotometer calibra- 
tion may have impacted this result. Before this analy- 
sis is complete, however, the wind vector must be 
considered in match-up pixel selection. If the wind 
speed is greater than a few meters per second, the cur- 
rent two-hour sunphotometer data averaging period 
suggests that pixels outside the scene area may need to 
be considered for an appropriate comparison. Con- 
versely, the sunphotometer averaging window could 


be reduced, although it is likely that this problem is 
not a factor in many of the analyses. 

The largest amount of match-up data available is 
from the Dry Tortugas. After extensive analysis of the 
Dry Tortugas data, there appear to be some problems 
with data quality, most likely the result of AERONET 
sunphotometer calibration or filter decay This prob- 
lem may be corrected in the latest match-up data set. 
In general, data from Bermuda suggest that the aerosol 
optical thicknesses measured by AERONET and 
SeaWiFs are highly compatible. In contrast, match-up 
data from Kaashidoo suggest potential problems with 
either the AERONET or SeaWiFS aerosol optical 
thickness measurements, or perhaps both. It is rea- 
sonable to conclude on the basis of these results that 
the atmospheric correction algorithm is highly suspect 
in the vicinity of Kaashidoo. Data from other match- 
up locations are being analyzed and preliminary re- 
sults suggest that SeaWiFS tends to measure a higher 
aerosol optical thickness at almost all locations than 
the AERONET. Additional ship and satellite aerosol 
optical thickness comparisons await a more extensive 
match up data base, which should be in place by the 
end of this year. 

23.3 WORK PLANNED 

There are still some issues to be addressed in the 
design of the MFRSSR. Controlled experiments need 
to be conducted on a motion simulator. The required 
equipment for these tests has been assembled and the 
experiments will be conducted when a second 
MFRSSR under construction for DOE is complete. In 
addition, pitch and roll information need to be incor- 
porated into the data processing algorithms used to 
process MFRSSR, especially when the instrument is 
deployed on small ships. 

Data from the R/V Ron Brown’s 7-month cruise 
for INDOEX and DOE Nauru-99 campaigns will be 
used to facilitate match-ups with SeaWiFS and sub- 
mitted to SeaBASS. In addition, two new MFRSSR 
units will be deployed in support of Nauru-99: one on 
the Japanese R/V Mirau and the other on the island of 
Nauru in the Tropical Western Pacific. These data 
will include over 200 measurement days aboard the 
R/V Ron Brown and additional 90 measurement days 
from the Mirai and Nauru. Several quality match-ups 
are anticipated from a variety of geographic locations 
and aerosol regimes. At the conclusion of INDOEX 
and Nauru-99, in August, all three MFRSSR units (the 
original and the two new units) will be offloaded at 
Hawaii and taken to Mauna Loa for a three week cali- 
bration, along with the microtops used in these ex- 
periments. Before shipment back to Brookhaven Na- 
tional Laboratory, one MFRSSR will be used in a 
collaborative cruise with Dr. John Porter at the Uni- 
versity of Hawaii. Hence, over 300 ship-board meas- 
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urement days are planned this fiscal year and, once tions, such as Kaashidoo, where the analysis suggests 
these data have been processed, satellite match-ups problems with the atmospheric correction algorithm 
will be made. The AERONET/SeaWiFS matchup In addition, attempts will be made to interpret the 

analysis is continuing. Data from all of the match-up AERONET- SeaWiFS match-up data in the context of 

sites will be analyzed. Particular emphasis will be the full set of atmospheric correction parameters as 
placed on interpreting the match-up data from loca- described by Gordon et al. (1980). 
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Figure 1 . Spectral optical thickness measurements made with the MFRSSR at 415 nm, 500 nm, 610 nm, 

660 nm, and 862 nm on January 27, 1999 over the Atlantic Ocean near the equator. Indicated in the boxes are 
instantaneous hand-held sun photometer measurements at 440 nm, 500 nm, 670 nm, and 870 nm. 
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Chapter 23 

Measurements of Aerosol, Ocean and Sky Properties 
at the HOT Site in the Central Pacific 


John N. Porter and Ricardo Letelier 
Hawaii Institute of Geophysics and Planetology, 
University of Hawaii, Honolulu, Hawaii 


23.1 INTRODUCTION 

Monthly cruises have been made to the Hawaii 
Ocean Time Series (HOT) site (~100 km north of 
Oahu) since October 1988. The goal of these cruises is 
to make hydrography, chemistry, and biology obser- 
vations (PI: Dave Karl and Roger Lucas). Measure- 
ments are typically made at a near coastal station 
(Kahe) on the first day to test the equipment and to 
obtain coastal shallow water (~ 1500m) observations. 
The second and third days are spent at the HOT site. 
On the morning of the fourth day, measurements are 
made at the HOT site and noon time measurements 
are made at the Hale-Aloha station near the mooring 
before returning to port by early the next morning. 
The locations of the three stations are: 

• Kahe (21.34° N, 158.27° W) 

• HOT (22.75° N, 158.0°W) 

• Hale-ALOHA buoy (22.43°N, 158.0°W). 

The routine HOT measurements are available 
during the summer following the year of the observa- 
tions. The 1998 measurements will therefore become 
available during the summer of 1999. The data sets 
can be obtained at the National Oceanographic Data 
Center (NODC) or from the Hawaii HOT web site 

(www.hahana.soesthawaii.edu/hot/hotJgofs.html). 

Selected parts of this data set will be submitted to the 
SeaBass archive as the data comes on line. 

23.2 RESEARCH ACTIVITIES 

Microtops Sun Photometers 

Beginning with the HOT89 cruise, aerosol optical 
depth measurements were made with two Microtops 
sun photometers. Examples of the data sets collected 
are shown in Figure 1. Aerosol optical depth meas- 
urements were made at 380, 440, 500, 675, 870 and 
1020 nm. Column integrated water vapor and ozone 
concentrations were also derived. These measure- 
ments have been submitted to the SeaBass archive. In 
general we find the aerosol optical depths to be 


roughly the same (-0.05) at all wavelengths for typical 
trade wind conditions where sea salt dominates. Asian 
dust periods also have a flat spectral aerosol optical 
depth but can have optical depths of 0. 1 and higher. 
Under VOG (volcanic smog) conditions, the shorter 
wavelengths have larger optical depths such as the 
case shown in Figure 1 (left panel). Calibration has 
been a major concern and numerous calibration efforts 
have been cariied out at Mauna Loa and Haleakala. 
The results of these calibrations are good and they are 
reported in a technical report submitted to SIMBIOS. 
The same report goes into detail on calibration prob- 
lems, measurement problems, and suggested protocols 
for using the Microtops hand held sun photometers. 
Comparisons between the ship sun photometer meas- 
urements and the SeaWiFS aerosol optical depths at 
870 are shown in Figure 2. In general the agreement is 
about V0.04 with the exception of two corrisponding 
cases which had flat spectral optical depths higher 
than 0.06, suggesting Asian dust may have been pres- 
ent. For VOG we found a systematic under-prediction 
of the aerosol optical depth by about 20%. 

Marine Shadowband Radiometer 

During the past year we have been working on a 
marine shadowband radiometer which will measure 
aerosol optical depths and downwelling irradiance at 
-20 wavelengths. This system uses a gimbaled cosine 
response detector and a rotating shadowband arm 
which shadows the detector. A CVI Laser spectrome- 
ter is used to collect the light and a PC104 computer is 
used to collect the data. The system has been tested on 
several cruises and the gimbal works well. An algo- 
rithm has been written to collect the total and diffuse 
light and the direct solar component is obtained from 
the difference. A cloud screening algorithm has also 
been developed which detects clouds by the variation 
in the downwelling irradiance just before and after the 
detector is shadowed. Although we believe we have 
made significant progress in many areas, we still have 
not resolved the cosine response problem as we would 
like near perfect cosine response with the cosine re- 
sponse known to 0.1%. We have tried diffusers made 
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of teflon, flashed opal, ground glass, spectralon, and 
integrating spheres. We now believe that the integrat- 
ing sphere or the spectralon offer the best perform- 
ance. Several new designs are being machined. 

Radiometric Calibration Facility 

In the past we have carried out radiance calibra- 
tion of our systems at the Moby calibration facility 
(thanks to Dennis Clarke). Now we have purchased a 
NIST traceable irradiance lamp and integrating sphere 
from Optronics. These systems were set up in a dark 
room at the University of Hawaii for optical calibra- 
tion. Our plans are to maintain this facility with rou- 
tine checks. Unfortunately, during a SIRREX com- 
parison, we found our calibration was approximately 
5.5% too low. We could not understand this differ- 
ence. Therefore we have purchased a shunt resistor 
which was recently calibrated by a NIST traceable 
calibration facility. We will soon use this calibrated 
resistor to check the current output of our power sup- 
ply. The Optronics power supply is the first suspect as 
it had a transistor failure which was replaced prior to 
our SIRREX comparison. We expect this problem will 
soon be solved. 

In Water Optics Measurements 

As part of this SIMBIOS effort, in water PRR 
(Profiling Reflectance Radiometer) and TSRB (Teth- 
ered Spectral Radiometer Buoy) measurements were 
begun on HOT90 cruise (Dec. 1997). PRR measure- 
ments are made by Karl et al. and processed by Ri- 
cardo Letalie. Calibration of this instrument was car- 
ried out and the coefficients have been submitted to 
the SeaBASS archive. A description of the data and 
calibration is available at 

http://picasso.oce.orst.edu/users/jasmine/ORSOO/. 

Aerosol Phase Function Measurements 

As part of our initial SIMBIOS proposal, we pro- 
posed to develop a polar nephelometer which could 
measure the aerosol phase function. Throughout the 


past year we have been developing this system and 
have recently collected our first measurements from 
this prototype instrument. Our design follows the sys- 
tem described by Winchester Jr. (in Optical Engi- 
neering, 1983, pg. 40). Our first measurements were 
made at Bellows beach directly downwind of small 
breaking waves. We have recently purchased a pulsed 
laser which will allow for phase function measure- 
ments at 3 wavelengths (1064, 532, 355 nm) The 
system is being designed to install in the door of a 
light aircraft and we expect to hope to make column 
aerosol phase function measurements in the future. 

Sky and Surface Radiometer 

We are developing a hand radiometer to sky radi- 
ance and above surface upwelling radiance The sky 
radiance will be used to derive aerosol size distribu- 
tion information and the surface radiance will be used 
to validate satellite measurements and to study the 
diffuse reflectance problem. This system uses a wide 
angle camera and a tilt meter to determine the viewing 
geometry and a spectrometer to measure the sky or the 
aureole. The scattering angle is derived from the posi- 
tion of the sun in the camera image and a tilt meter. 
The system will be in a weather proof enclosure using 
a PC/ 104 computer to control the two axis motor, the 
camera, and the spectrometer. We had expected to 
have a prototype system working by now but delays 
have pushed it back to Spring 1999. 

23.3 WORKPLAN 

We have three goals: 1) continue collection of 
existing measurements on the HOT cruises, 2) finish- 
ing the construction of new instruments (shadowband 
radiometer, polar nephelometer, and sky radiometer), 
and 3) cany out modeling efforts. The modeling effort 
will focus on developing new algorithms to derive 
column aerosol properties from surface measurements 
and estimating the impact of aerosol on upwelling 
radiance at the top of the atmosphere. 



AEROSOL OPTICAL DEPTH 


SMBIOS Project Annual Report 





SIMBIOS Project Annual Report 


Chapter 24 

Assessment of the Contribution of the Atmosphere to 
Uncertainties in Normalized Water-Leaving Radiance: 
a Combined Modelling and Data Analysis Approach 


Knut Stamnes and Bingquan Chen 
Geophysical Institute 

University of Alaska Fairbanks , Fairbanks , Alaska 


24.1 INTRODUCTION 

Atmospheric correction, i.e. the process of deriv- 
ing the normalized water-leaving radiances from sat- 
ellite-measured imagery of the oceans, has been es- 
tablished as one sub-goal of the SIMBIOS program 
for which the main purpose is to provide routine, long 
term ocean color data from a variety of international 
ocean color missions. Our research for the SIMBIOS 
program is focused on the study of atmospheric cor- 
rection. The overall objective is to use a sophisticated, 
state-of-the-art radiative transfer model for the cou- 
pled atmosphere-ocean system, with existing and 
planned algorithms, and in conjunction with data 
taken by other investigators in the SIMBIOS program, 
to help quantify the uncertainties in the water-leaving 
radiance due to inadequacies in atmospheric correc- 
tions applicable to OCTS, POLDER, SeaWiFS, 
MODIS, and other ocean color sensors. 

Efficient radiative transfer codes for the atmos- 
phere-ocean system are not generally available, and 
the popular and intuitively appealing Monte Carlo 
method (Morel and Gentili, 1991) is too time- 
consuming to be used in routine retrieval algorithms. 
To solve the radiative transfer equation for the cou- 
pled atmosphere-ocean system, we employ the dis- 
crete-ordinate method (Jin and Stamnes, 1994). This 
model for computing the transfer of radiation in the 
coupled atmosphere-ocean system is reliable and effi- 
cient. It has the accuracy of the Monte Carlo method, 
but overcomes its computational inefficiency. 

24.2 RESEARCH ACTIVITIES 

Water-leaving Radiance at N1R Wavelengths 


provide realistic simulations of the atmospheric con- 
tribution to the TOA radiance. We also employ a bio- 
optical model of ocean water with ocean particle con- 
centrations of 1.5, 5.0, and 10.0 mg/m 3 , respectively, 
and we employ phase functions for phytoplankton 
particles (i.e., cyanobacteria and algae, Volten et al., 
1998) in the simulations conducted in this study. 

Forward modeling is used to simulate or predict 
radiances that ocean color sensors would measure 
From the deviation between the TOA radiances pre- 
dicted by a model of the coupled atmosphere-ocean 
system (Jin and Stamnes, 1994) and the same model in 
which scattering by ocean particles (and hence water- 
leaving radiances) at NIR wavelengths has been ig- 
nored, one can determine under what conditions the 
water-leaving radiance at MR wavelengths can be 
ignored in atmospheric corrections. For each sensor 
channel we use forward modeling simulations to ac- 
curately predict the impact of both atmospheric aero- 
sols and biogenic constituents in the ocean on TOA 
radiances for any desired sun-satellite geometry, solar 
zenith angle 0 0 > sensor viewing polar angle 0, and 
relative azimuthal angle <)> with respect to the sun. 

To describe scattering by particles in the ocean 
we employ the Henyey-Greenstein phase function p HG 
(Henyey and Greenstein, 1941): 


PHG^> cos0 ) 




\3/2 


L + g - 2g cos © J 


This phase function has the desirable feature that 
it yields complete forward scattering when the asym- 
metry factor g=l, isotropic scattering when g=0, and 
complete backward scattering when g=-l. The linear 
combination 


The main goal of this study is to assess the un- 
certainty or error incurred in atmospheric correction 
by ignoring the water-leaving radiance caused by 
scattering from particles in the near-surface water at 
NIR wavelengths between 660-870 nm. Aerosol mul- 
tiple scattering is included in the forward modeling to 


p ( cos ©) = bp HG (g, cos ©) + (1 - b)p HG (g \ cos 0) 

is sometimes used to simulate a phase function with 
both a forward and a backward scattering component 
(g>0 and g'<0). Here 0 is the scattering angle, b is the 
fraction of forward-scattered light (0<b<l), and g and 
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g' are usually different. Based on laboratory measure- 
ments by Volten et al. (1998), we have used equation 
(2) with different values of g, g\ and b to simulate 
phase functions corresponding to ocean waters with 
different types of particles. 

The measurements by Volten et al. (1998) were 
carried out at 633 nm. Lacking better information we 
assume for our present purpose that they can be used 
to approximate the phase functions at wavelengths 
between 630 nm and 870 nm. Previous studies show 
that far from sources of pollution and/or sources of 
desert aerosols the aerosol optical depth at 865 nm 
over the Pacific Ocean lies in the range between 0.08 
and 0.11 (Villevalde et al., 1994). In this study we 
employ a Tropospherical aerosol model with relative 
humidity of 80% (Shettle and Fenn, 1979) and use a 
Mie code to calculate the aerosol optical properties. 
From forward-modeling results based on these syn- 
thetic phase functions we conclude that: 

• The TOA radiance deviation depends on the 
scattering phase function of particles in the near- 
surface water. 

• As the aerosol optical depth decreases, the TOA 
radiance deviation increases. When the aerosol 
optical depth is large (e g. x a =0.2), TOA radi- 
ances, predicted by a model in which water- 
leaving radiances at NIR wavelengths are ig- 
nored, do not deviate significantly from those 
obtained when including scattering by ocean par- 
ticles. Thus, at large aerosol optical depths, we 
may ignore water-leaving radiances at some NIR 
bands (e g. 865 nm) without introducing signifi- 
cant errors. 

• The TOA radiance deviation increases with in- 
creasing ocean particle concentration. When the 
ocean particle concentration is large and the aero- 
sol optical depth is low, water-leaving radiances 
at NIR wavelengths contribute substantially to 
TOA radiances and should not be ignored. 

• The TOA radiance deviation that results from 
neglecting water-leaving radiances depends on 
the sun-satellite geometry. 

• The TOA radiance deviation at A,=765 nm is 
larger than that at X =865 nm. 

• The TOA radiance deviation at X =665nm is 
much larger than that at X =865nm. Since the 665 
nm wavelength has been used in atmospheric cor- 
rection for CZCS, this deviation can have led to 
significant uncertainties in ocean color retrievals 
from the CZCS measurements. 

To avoid the uncertainties incurred by ignoring 
water-leaving radiances at NIR wavelengths retrievals 
and atmospheric corrections in ocean color imagery 
should be based on a radiative transfer model for the 
coupled atmosphere-ocean system (Jin and Stamnes, 
1994). 


Atmospheric Correction Algorithm 

The existing atmospheric correction algorithm is 
mainly based on the approximate exponential relation- 
ship of the single scattering aerosol reflectance over 
the range 412-865 nm. Two types of lookup tables for 
N selected aerosol models then constitute the main 
part of the algorithm. The first type of lookup table is 
for the coefficient ( X , 865) which relates the single 

scattering aerosol reflectance at wavelength X to that 
at 865 nm. The second type of lookup table is for the 
coefficient k[X y p^,)] which makes a correction for 
multiple scattering. For each type of lookup tables, 
about twelve candidate aerosol models are employed; 
and table entries are computed for each of these mod- 
els for solar zenith angles between 0°and 80° in in- 
crements of 2.5°, for 33 values viewing zenith angles, 
and for eight values of the aerosol optical depth be- 
tween 0.05 and 0.8. The total number of separate so- 
lutions to the radiative transfer equation used in the 
preparation of these lookup tables exceeds 33,000 
(Gordon, 1996). Thus, the computational effort in- 
volved in the construction of these lookup tables is 
very large. 

Based on a state-of-the-art radiative transfer 
model for the coupled atmosphere-ocean system, we 
simulate the process of atmospheric correction. Using 
four types of aerosol models (Tropospheric, Maritime, 
Coastal and Urban with RH = 50, 80 and 98 %, 0= 7°) 
we find that there exists an approximate exponential 
relationship of the aerosol reflectance over the range 
412-865 nm not only in the single scattering case, but 
also in the multiple scattering case. Thus, a simple and 
efficient algorithm for atmospheric correction is pro- 
posed and discussed by comparing it with the existing 
algorithm. In our algorithm, the aerosol contribution 
in the visible in the multiple scattering can be ex- 
trapolated directly from that in the near-infrared by 
use of the approximate exponential relationship 
Therefore it is not necessary to construct lookup tables 
for the coefficient k[X, p^X)]. The new coefficient e 
m5 (X, 865), which represents the ratio of the multiple 
scattering aerosol reflectance at X to that at 865 nm, is 
introduced in the new algorithm instead of e ms ( X , 
865) used in the existing algorithm. Future testing of 
this algorithm will be carried out by use of satellite 
data and field measurements. The most important con- 
clusion of this study is that one should examine 
whether or not water-leaving radiances at NIR wave- 
lengths can be ignored when developing algorithms 
either for the retrieval of atmosphere-ocean parame- 
ters or for atmospheric corrections. Otherwise, errone- 
ous results can be obtained if such algorithms are ap- 
plied to satellite-received data, even for open oceans 
A simple and efficient algorithm for atmospheric cor- 
rection is also proposed and compared with the exist- 
ing algorithm. 
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Chapter 25 

Validation of the Water-Leaving Radiance Data Product 

Oleg V. Kopelevich 
Ocean Optics Laboratory 

P.P. Shirshov Institute of Oceanology , Moscow , Russia 


25.1 INTRODUCTION 

The overall objective of the research is to assess 
an accuracy of retrieval of the water-leaving radiance 
from satellite ocean color data. Specific goals include 
an assessment of accuracy of the atmospheric correc- 
tion algorithms for variety of meteorological and 
oceanological conditions, effects of the rough sea sur- 
face reflectance and of the small-scale spatial inho- 
mogeneities in the subsurface layer (sub-pixel vari- 
ability), an accuracy of the water-leaving radiance 
values obtained from in situ measurements. Combined 
approach is used which includes direct comparison 
between the values of water-leaving radiance derived 
from satellite data and in situ measurements, computer 
simulation, theoretical analysis. 

25.2 RESEARCH ACTIVITIES 

Field Methods and Data 

The instrumental set constructed in the 
P.P. Shirshov Institute of Oceanology Russian Acad- 
emy of Sciences (SIO RAS) for field studies com- 
prises deck and floating spectroradiometers and 
monitor photometer. The deck spectroradiometer is a 
modification of the device described before (Goldin et 
al. 1983); the floating spectroradiometer and the 
monitor photometer are new instruments (Artemyev et 
al., 1999). The floating spectroradiometer measures 
the upwelling radiance and downwelling irradiance in 
the spectral range 390-700 nm with spectral resolution 
2.5 nm; the upwelling radiance is measured just be- 
neath the sea surface and the downwelling irradiance 
just above the sea surface at a distance about 50 m 
from ship; thus the difficulties connected with influ- 
ence of sun glints, correction for light reflected at the 
sea surface, influence of ship shadow are precluded. 
The monitor photometer is designed to check calibra- 
tion of irradiance at 554 nm and measure continuously 
the surface irradiance at this wavelength during spec- 
tral measurements by the spectroradiometers. The 
radiometric accuracy of the deck and floating spectro- 
radiometers and the monitor photometer is respec- 
tively 10%, 5% and 3%. The field testing of the float- 
ing spectroradiometer were carried out on the oceano- 


grafic platform of the Marine Hydrophysical Institute 
NASU in Katsiveli, September 1996. They showed an 
ability of the instrument to work under bad weather 
conditions (strong wind, cloudiness) when measure- 
ments by the deck radiometer are impracticable. 

The field studies under satellite ocean color sensors 
SeaWiFS and MOS-IRS were performed during the 
scientific cruise of RV Akvanavt in the Black and Ae- 
gean Seas, October 6-24, 1997 (Artemyev et al., 
1999). Along with measurements of spectral upwell- 
ing radiance and surface irradiance by the deck and 
floating radiometers, the studies included measure- 
ments of vertical profiles of nadir upwelling radiance 
and downwelling irradiance by means of the sub- 
mersible radiometer MER-2040 (these measurements 
were performed by specialists from the Polish Institute 
of Oceanology, Sopot); determination of the atmos- 
pheric spectral optical thicknesses by means of the sun 
photometer; measurements of vertical profiles of the 
seawater beam attenuation coefficient by the sub- 
mersible transmissometer, sampling and conservation 
of water samples for determination of phytoplankton 
pigment concentration. Five drift stations were earned 
out during the cruise, three in the Black Sea and two 
in the Aegean Sea. The information on the stations is 
given in Table 1. 

Absolute radiometric calibration of the radiance 
and irradiance channels was performed before the 
expedition. During the cruise, intercalibration of the 
floating spectroradiometer and MER-2040 was carried 
out; an example of intercalibration is shown in Fig 1. 
Discrepancy between values of L u and E d measured by 
these two instruments was mainly in limits of 5%. 

In August-September 1998 the field studies in the 
eastern part of the Barents Sea were carried out during 
the 13 th and 14 th scientific cruises of R/V Akademik 
Sergey Vavilov. They included the measurements 
listed above, excluding measurements by MER-2040, 
and additionally measurements of the diffuse attenua- 
tion coefficient at wavelength 530 nm by the IQ - 
meter and also determinations of the particulate matter 
concentration, primary production, dissolved organic 
carbon concentration, CTD profiles. The information 
on the stations with optical measurements is given in 
Table 2. 
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25.3 RESEARCH RESULTS 

The validation cruise of RV Akvanavt in the 
Black and Aegean Seas 

The validation of the SeaWiFS atmospheric cor- 
rection algorithm was performed with measured data 
from Station 2 and Station 4 where the measurements 
were carried out under appropriate weather conditions. 
Before validation, an agreement of the values of the 
top-of-the-atmosphere radiance L/Aj) measured by 
SeaWiFS sensor and calculated by a Monte Carlo 
technique with in situ data on pfAJ and r a (ty was 
checked. The procedure developed for the Monte 
Carlo calculation is described elsewhere (Kopelevich 
et al. 1998; Burenkov et al., 1999). The results of 
comparison between the measured and calculated val- 
ues of L t (Ai) are given in Table 3. As it is seen, the 
agreement is quite well for Station 2; it should be 
borne in mind that the observed discrepancies result 
not only from errors of the measurements but from 
also some assumption taken for the calculation. The 
agreement is worse for Station 4; it can be explained 
by some atmospheric inhomogeneity due to the 
proximity of the cloud edge to the pixel with this sta- 
tion. 

The SeaWiFS atmospheric correction algorithm 
was validated by comparison of the Lm/^i) and r a (Ai) 
values derived from SeaWiFS data and from in situ 
measurements. The results of comparison are given in 
Table 4. As it is seen, the agreement is quite satisfac- 
tory both for the LmjfAJ and z a . 

The SeaWiFS algorithm for retrieving chlorophyll 
concentration and the analytic algorithm developed by 
the specialists of SIO RAS were validated by com- 
parison of retrieved and measured data. The SeaWiFS 
algorithm overestimates significantly the chlorophyll 
concentration for Case-2 waters with predominance of 
absorption by yellow substance; the analytic algorithm 
results to reasonable agreement in all cases (Burenkov 
et al. in press). 

The Akademik Sergey Vavilov Cruises in the 
Barens Sea 

The preliminary analysis revealed some cases of 
incorrect SeaWiFS data products; an example is given 
in Figure 2, The negative values of 1^443) are seen 
over most part of the image (left); the chlorophyll 
concentration derived from SeaWiFS data agree rea- 
sonably with the measured values only in open part of 
the Barent Sea and differ drastically (by a factor 20 
and more) in the Pechora Basin (right image). 


The measured spectral dependences of the radiance 
reflectance exhibit a wide variety of forms (Figure 3): 
the spectra measured in the north-eastern part of the 
Barents Sea have the maximum near 420 nm and the 
ratio Zww(490)/Z,# w (555) about 3 (St. 1209), whereas in 
the southern part of the Pechora Basin they are re- 
spectively near 580 nm and less than 0.5 (St. 1095). 
According to the SeaWiFS data, the chlorophyll con- 
centration varies over a wide range of values: from 
less than 0.5 mg-m 3 in the northern part to more than 
20 mg-m' 3 in the Pechora Basin (Figure 2, right im- 
age). The in situ data show much more narrow range 
of its changes. The further analysis is planned to es- 
tablish the reasons for the revealed errors resulted 
from the SeaWiFS algorithms. 

Effects of the Sub-pixel Variability 

The mathematical simulation of effects of the 
small-scale inhomogeneities in the subsurface layer 
(sub-pixel variability) on the retrieval results from 
satellite data has been performed. The horizontal dis- 
tribution of chlorophyll concentration with the sharp 
change was assumed, and the averaged spectral up- 
welling radiance was calculated. Then the averaged 
spectra were used to retrieve the values of chlorophyll 
concentration. Comparison between the obtained and 
true mean values has shown that discrepancies be- 
tween them are no more 5 % if the amplitude of 
change is less 50%, but they increase up to 15% with 
the twofold change and 45% with the fourfold. Such 
sharp changes are real in the coastal waters, and po- 
tential errors should be taken into account. 

25.4 WORK PLAN 

• The further analysis of data obtained during the 
13 th and 14 th cruises of R/V Akademik Sergey 
Vavilov in the Barent Sea, August-September 
1998. 

• Preparation of the scientific equipment for work 
in the 15-th cruise of R/V Akademik Sergey Vav- 
ilov in the Norwegian and Greenland Seas 

• The field studies in the 15 th cruise of R/V 
Akademik Sergey Vavilov in the Norwegian and 
Greenland Seas, May-June 1999. 

• Analysis of data obtained during the 15 th cruise of 
R/V Akademik Sergey Vavilov in the Norwegian 
and Greenland Seas, May-June 1999. 

• Validation of the SeaWiFS algorithm for retriev- 
ing the chlorophyll concentration in the Black Sea 
by comparison of the satellite and available in situ 
data. 
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Figure 1. Results of intercalibration between the floating spectroradiometer and MER-2040 at St 2; upper 

curves - surface irradiance, lower - upwelling radiance; ♦ - floating spectroradiometer, • - MER-2040. 



Figure 3. Spectra of the water radiance reflectance measured in different parts of the Barents Sea during 

the 13-th and 14-th cruises of R/VAcademik Sergey Vavilov: St. 1209 -the north-eastern part of the Barents Sea; 
St. 1088 and St. 1095 - the northern and southern parts of the Pechora Basin (see Table 2). 
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Table 1. The information on the stations of Akvanavt ’97. 


Date 

St. 

Coordinates (degrees) 
at the measurement 
moments 

Local 

time 

Difference with 
GMT, h 

10/07/97 

1 * 

42.5 N, 39.5 E 

09:20- 14:40 

+3 

10/08/97 

2 V* 

43.0 N, 35.6 E 

10:05-15:40 

+3 

10/09/97 

3 

42.9 N, 31.6 E 

10:48 - 15:11 

+3 

10/11/97 

4 * ** 

39.3 N, 25. IE 

11:30-14:55 

+3 

10/16/97 

5 ** 

39.6 N, 25.8 E 

10:19-14:58 

+3 


* - under SeaWiFS; ** - under MOS-IRS. 


Table 2. The information on the stations of R/V “ Akademik Sergey Vavilov” 1 3 and 14 cruises. 


Date 


Coordinates (degrees) 
at the measurements 
moments 

Local time 

Difference with 
GMT, h 

08/1 1/98 


70.42 N, 47.58 E 


+4 

08/12/98 

■MSS 

70. 18 N, 52.42 E ' 


+4 

08/14/98 

tmm 

68.97 N, 58.47 E 


+4 

08/14/98 


69.07 N. 58.42 E 1 

16:20-17:30 

+4 

08/15/98 

■ 

69. 10N, 58.03 E 


+4 

08/19/98 

III EM 

69.09 N, 58.29 E 


+4 

08/19/98 

mi 

69.09 N, 58.29 E 

13:55-14:25 

+4 

08/23/98 

1123 

69.06 N, 58.16 E 

07:00 - 10:00 

+4 

08/23/98 

1125* 

69.50 N, 57.25 E 

13:10 - 14:20 

+4 


| 

69.67 N, 57.24 E 


+4 

liEfcirkii 

■ 

69.77 N, 56.28 E 


+4 


■ 

70.33 N, 55.30 E 


+4 

08/29/98 

HEM 

69.33 N, 50. 12 E 

08:10-09:30 

+4 

08/30/98 

WSSSM 

69.62 N, 50.57 E 

10:30 - 14:50 

+4 

08/31/98 

■IBM 

70.54 N, 52.79 E 


+4 

09/02/98 

wmm 

68.40 N, 50.70 E 


+4 

\mmum 

inn 

69.25 N, 41.00 E 

08:40 - 10:40 

+4 

09/11/98 

1183 


mmiBniM 

+4 

09/12/98 

1190 

73.25 N, 41.00 E 

10:00-12:00 

+4 

09/13/98 

1196 

74.75 N, 41.00 E 

09:30- 11:30 

+4 

09/15/98 

1209 

78.00 N, 41.00 E 

10:20-13:00 

+4 

09/27/98 

1281 

76.00 N, 42.27 E 

10:20-12:40 

+4 


* - under SeaWiFS 


Table 3. Comparison between the SeaWiFS measured and Monte Carlo calculated values of the 
top-of-the-atmosphere-radiance U(Xj\ A is the relative discrepancy. 



Station 2 

Station 4 


X, 

L, SeaWiFS, 

L, MC, 

A, 

L SeaWiFS, 

L, MC, __ 

A, 

nm 

uW cm' 2 sr' 1 nm' 1 

uW cm' 2 sr' 1 nm' 1 

% 

uW cm" 2 sr’ 1 nm -1 

jiW cm' 2 sr' 1 nm' 1 

% 

412 

7.40 

7.42 

0.3 

8.10 

7.89 

2.6 

443 

6.52 

6.57 

0.7 

7.24 

6.86 

5.0 

490 

4.84 

4.83 

0.2 

5.35 

5.16 

3.5 

510 

4.12 

4.08 

1.0 

4.50 

4.21 

6.4 

555 

2.93 

2.80 

4.4 

3.18 

2.96 

6.9 

670 

1.25 

1.16 

7.2 

1.55 

1.44 

7.1 

765 

0.55 

0.59 

7.2 

0.77 

0.73 

4.8 

865 

0.36 

0.33 

8.3 

0.55 

0.43 

22 
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Table 4. Comparison between values derived from SeaWiFS and in situ data for the normalized 
water-leaving radiance Lwu(ty and the aerosol optical thickness r a (865)\ A is the relative discrepancy. 



Station 2 

Station 4 

A,, nm 

Lwn SeaWiFS, 

Lwn measured, 

A, 

Lwn SeaWiFS, 

Lwn measured, 

A, 


pW cm' 2 sr' 1 nm' 1 

pW cm* 2 sr* 1 nm* 1 

% 

pW cm’ 2 sr* 1 nm* 1 

pW cm* 2 sr' 1 nm" 1 

% 

412 

0.536 

0.492 

9 

0.770 

0.732 

5 

443 

0.752 

0.708 

6 

0.886 

0.891 

0.6 

490 

1.036 

0.998 

4 

0.882 

0.784 

6 

510 

0.991 

0.925 

7 ~ 

0.592 

0.537 

10 

555 

0.770 

0.772 

0.3 

0.271 

0.272 

0.4 

t„(865) 

0.065 

0.095 

31 

0.160 

0.142 

13 
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Chapter 26 

OCI Data Validation in the Waters Adjacent to Taiwan 


Hsien-Wen Li, Chung-Ru Ho and Nan- Jung Kuo 
Department of Oceanography 
National Taiwan Ocean University , Keelung, Taiwan 


26.1 INTRODUCTION 

Taiwan successfully launched her first experi- 
mental satellite, ROCSAT-1, on January 27 1999. It is 
a three-axis stabilized, low-earth orbit satellite with a 

35° inclination. ROCSAT-1 is designed to carry out 
three scientific experiments: ocean color imaging, 
ionospheric plasma and electrodynamics, and Ka-band 
communication. 

The payload instrument for ocean color imaging 
on ROCSAT-1 is called Ocean Color Imager (OCI). 
OCI is a multi-spectral push-broom imager, which is 
designed to map six reflected spectral radiances from 
ocean surfaces. The six spectral bands are identical to 
the six of the eight spectral bands on SeaWiFS. The 
band characteristic comparison of OCI and SeaWiFS 
is listed in Table 1. 

In order to have successful operation of OCI, 
there are two scientific groups under the OCI project, 
the Science Team (ST) and the Science Data Distribu- 
tion Center (SDDC). The purposes of the ST are to 
calibrate the sensor, to correct for the atmospheric 
influence, to develop bio-optical algorithms, and to 
calibrate and validate OCI data. The SDDC at Na- 
tional Taiwan Ocean University (NTOU) is responsi- 
ble for processing OCI data and then distributing it to 
users. 

OCI will perform imaging at any time between 
9:00 and 15:00 local time The normal mission opera- 
tion of OCI is to map ocean surface pigments when- 
ever weather permits. 

Ideally, OCI will acquire ocean surface pigment 
data when cloud coverage is less than 50% of the area 
along the track. The actual OCI operation will ulti- 
mately be decided by the OCI-ST. A preliminary op- 
eration priority has been established based on geo- 
graphic locations, request sequences, weather condi- 
tions, and cross calibration and field validation areas. 

For calibration and validation purposes, images 
are schedule to be taken when the OCI tracks overlap 
with SeaWiFS, or the calibration buoy location or 
location of the experimental cruises. 


26.2 RESEARCH ACTIVITIES 

It is believed that the OCI calibration at each 
wavelength will change in some unpredictable man- 
ners as a function of time. Experience with previous 
sensors, such as Coastal Zone Color Scanner (CZCS), 
has shown that it is very difficult to determine a sen- 
sor’s calibration once it has been launched (Mueller 
and Austin, 1992). Unlike the SeaWiFS instrument, 
OCI has no on-board solar calibrator. Therefore, a 
cross-calibration program with SeaWiFS will be per- 
formed after OCI is available. To validate OCI and 
SeaWiFS data, field measurements of optical proper- 
ties are carried out with a Tethered Spectral Radi- 
ometer Buoy (TSRB-II) and a SeaWiFS Profiling 
Multi-channel Radiometer (SPMR) made by Satlantic 
Company. TSRB-II is an optical buoy, which can 
measure the in-water upwelling radiance just beneath 

the sea surface L u (0 , X) and the incident spectral 

irradiance above the sea surface E s (X) . SPMR is a 
profiling radiometer that can measure the in-water 
upwelling radiance L U (Z,X) and downwelling 
spectral irradiance E d (Z , X) with depth Z. 

To obtain L W (X) that is measured by the OCI it 

is necessary to propagate L u (0 , X) upward through 
the sea surface as 

L w (A) = L u (0~ , A) (1) 
n wW 

where PC, &) and (^) are the Fresnel reflec- 
tance at the solar zenith angle and the refractive index 
for seawater, respectively. In order to remove the in- 
fluence of view angle, sun angle and the solar irradi- 
ance, the water-leaving radiance is generally trans- 
ferred to normalized water-leaving radiance L W n as 

L m W = L.W^P- ( 2 , 

EM) 
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where F 0 (/l) denotes the mean extraterrestrial solar 
irradiance (Neckel and Labs, 1984). For calculating the 
remote sensing reflectance just above the sea surface 

R rs ( 0 + > i^e following equation was used 

*„(0*,-t) = 0.544(0-, X)/E,{A) (3) 

Since the band characteristics of OCI are similar 
to those of SeaWiFS, some of bio-optical algorithms 
from the SeaWiFS Bio-optical Algorithms Mini- 
Workshop (SeaBAM) are evaluated using the data sets 
that are collected in the water adjacent to Taiwan. The 
evaluated algorithms are listed in Table 2. We have 
collected a total of 29 in situ data including optical 
properties and water samples in order to derive the 
relationship between optical properties and chloro- 
phyll a concentration. All samples, except the surface 
seawater, were collected by Go-Flo 2.5/5 L sampler 
which was attached on the rosette multi-sampler. Sur- 
face seawater was collected with a bucket. One liter of 
seawater was immediately filtered by Whatman GF/F 
(25 mm) filter, under pressure (<0.2 bar). To deter- 
mine the chlorophyll a concentration, the filtered 
sample was frozen in liquid nitrogen on board, while 
in the laboratory the N, N-dimethylformamide (DMF) 
extracting solvent was added in darkness at -20°. 
Compared with the acetone extracting solvent, DMF 
has some advantages (Suzuki and Fujita, 1986; Porra 
et al, 1989; Suzuki and Ishimaru, 1990) and was 
therefore used in this study. The measurement of chlo- 
rophyll a was undertaken in the laboratory according 
to the fluorescence method described by Strickland 
and Parsons (1972). 

26.3 RESEARCH RESULTS 

A total of 29 bio-optical data were made in the 
time period of 1998 in the Kuroshio region near Tai- 
wan for testing the SeaBAM algorithms The field 

measurements include E s (Z) ? L u (Z,Z) > 

E d (Z,A) f and water sample for determining the 


chlorophyll a concentration. The Lwn is computed 
using Eqs. (1), and (2). The Rrs is obtained from Eq. 
(3). Most of the water samples were taken in the Case 
I water. The range of measured chlorophyll a concen- 
tration is from 0.08 ug/1 to 1.25 ug/I in this data set 
The statistical results between calculated chlorophyll a 
concentration using some of SeaBAM algorithms and 
measured chlorophyll a concentration from water 

2 

samples are listed in Table 3. The results show that R 
of Aiken-C, CalCOFI 2-band linear, CalCOFI 2 -band 
cubic, Morel 2, and SeaWiFS algorithms are higher 
than 0.8. The RMS is smaller than 0.015 for these 
methods. This suggests that the SeaWiFS algorithm is 
also suitable for the waters adjacent to Taiwan 

26.4 WORK PLAN 

In situ measurements of optical properties and 
chlorophyll a concentration will be continuous in the 
waters adjacent to Taiwan, especially in Taiwan Strait 
and East China Sea, The schedule for normal opera- 
tion of OCI is from April 1999. Therefore the inter- 
comparison between SeaWiFS and OCI will be the 
most important plan for next period. Some of 
SeaWiFS project members will also take bio-optical 
data for validating OCI data. OCI data will be avail- 
able in April 1999, 
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Table 3. Summary of statistical results of algorithm evaluations. 


Method 

1 

2 

3 

4 

5 

6 

7 

8 

9 

R 2 

0.859 

0.833 

0.816 

0.772 

0.825 

0.769 

0.794 

0.792 

0.848 

RMS 

0.012 

0.014 

0,015 

0.019 

0.015 

0.020 

0.017 

0.018 

0.013 
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Table L A Comparison of characteristics between SeaWiFS a nd OCI. 



SeaWiFS 

OCI 

Inclination 

98.25° 

35° 

Altitude (km) 

705 

600 

Period (min) 

98.9 

96.6 

Orbital repeat time (days) 

16 

52 

Spectral bands (nm) 

B1 402-422 

B1 433-453 


B2 433-453 

B2 480-500 


B3 480-500 

B3 500-520 


B4 500-520 

B4 545-565 


B5 545-565 

B5 660-680 


B6 660-680 

B6 845-885 


B7 745-785 

B7 545-565 


B8 845-885 


2 

Nadir pixel (m ) 

1130 x 1130 

800 x 800 

Swath width (km) 

2801 

702 

Redundancy 

No 

555 nm 

Color sensing 

Scanner 

push broom 

Crossing equator time 

12:00 

9:00 ~ 15:00 

(local time) 



Bits 

10 

12 

Tilt 

-20°, 0°, 20° 

No 

Launch date (year/month) 

1997/8 

1999/1 


Table 2, Bio-optical algorithms performed for the data taken in the waters adjacent to Taiwan. 


Method 

Algorithms 

Band ratio(R), Coefficient(s) 

1 . Aiken-C 

C=EXP(a0+a 1 *ln(R)) 
C=(R+a2)/(a3+a4*R) 

R=L wn490/Lwn5 5 5 
a=f0.464, -1.989, -5. 29,0. 7 19,-4.23] 

2. CalCOFI 2- band Linear 

C=10 A (a0+al*R) 

R=log(Rrs490/Rrs5 55) 
a-fO.444,-2.431] 

3. CalCOFI 2- band Cubic 

C= lO^aO+a 1 *R+a2*R A 2+a3 *R A 3) 

R=log(Rrs490/Rrs555) 
a=[0.450 > -2.860,0.996, -0.36741 

4. Morel 1 

C=10 A (a0+al*R) 

R=Log(Rrs443/Rrs5 5 5) 
a=f0.2492,-1.7681 

5. Morel 2 

C=EXP(a0+al*R) 

R=Ln(Rrs490/Rrs5 5 5 ) 
a=ri. 077835,-2.542605] 

6. Morel 3 

C= 1 0 A (a0+a 1 *R+a2 *R A 2+a3 *R A 3) 

R=Log(Rrs443/Rrs5 55) 
a=f0.20766,-l. 82876,0.75885, -0.739791 

7. Power law model 

C=EXP(a0+al*R) 

R=Ln(Lwn44 3/L wn5 5 5) 
a=f-0.0177,-1.2331 

8. Carder Ratio 

Log(C)=a0+(a 1 +a2 *log(R)) *log(R) 

R=Rrs4 90/Rrs5 5 5 
a=f0.0469, -1.0129, -0.7416] 

9. SeaWiFS 

C= 1 0 A (a0-a 1 *R+a2 *R A 2+a3 *R A 3 )+a4 

R=log(Rrs490/Rrs5 55) 
a=[0.2974,-2.2429, 0.8358, -0.0077, -0.0929] 
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GLOSSARY 

CHORS 

Center for Hydro-Optics and Remote 

ACE 

Aerosol Characterization Experiment 


Sensing (San Diego State University) 

ADEOS 

Advanced Earth Observation Satellite 

CICESE 

Centro de Investigation Cientifica y de 


(Japan) 


Edcacion Superior de Ensenada 

AERONET 

Aerosol Robotic Network 


(Mexico) 

AM-1 

Not a acronym, used to designate the 

CIMEL 

The name of a sun photometer manu- 


morning platform of EOS 


facturer 

AOP 

Apparent Optical Properties 

CNES 

Centre National d’Edudes Spatialle 

AOT 

Aerosol Optical Thickness 

CONICIT 

Consejo Nacional de Investigaciones 

APV 

Autonomous Profiling Vehicle 


Cientifica y Technologicas (Venezuela) 

ARGOS 

Not an acronym, but the name given to 

CRAM 

Conditional Relaxation Analysis Method 


the data collection and location system 

CTD 

Conductivity-T emperature-Depth 


on the NOAA Operational Satellites. 

CZCS 

Coastal Zone Color Scanner 

ASCII 

American Standard Code for Information 

DAAC 

Distributed Active Archive Center 


Interchange 

DLR 

Deutshe Forshungsanstalt f r Luft-und 

AVHRR 

Advanced Very High Resolution 


Raumfarhrt (German Aerospace Center) 


Radiometer 

DMF 

Dimethylformamide 

AVIRIS 

Advanced Visible and Infrared Imaging 

DOC 

Dissolved Organic Carbon 


Spectrometer 

DoD 

Department of Defense 

BATS 

Bermuda Atlantic Time-series Study 

DOE 

Department of Energy 

BBOP 

Bermuda Bio-Optics Profiler 

ECOHAB 

Ecology of Harmful Algal Blooms 

BBSR 

Bermuda Biological Station for Research 

EEZ 

Exclusive Economic Zone 

BNL 

Brookhaven National Laboratory 

EORC 

Earth Observation Research Center 

BTM 

Bermuda Test Mooring 

EOS 

Earth Observing System 

Cal/Val 

Calibration and Validation 

FFP 

Firm-Fixed Price 

CalCOFI 

California Cooperative Oceanic 

FOV 

Field of View 


Fisheries Investigation 

ftp 

File transfer protocol 

CALVAL 

Calibration Validation 

FWHM 

Full Width Half Maximum 

CARICO 

Carbon Retention in a Colored Ocean 

GAC 

Global Area Coverage, coarse resolution 

Case-1 

Water whose reflectance is determined 


satellite data with a nominal ground re- 


by absorption. 


solution at nadir of approximately 4 Km 

Case-2 

Water whose reflectance is significantly 

GB 

Gigabyte, or about one billion bytes 


influenced by scattering. 

GF/F 

Not an acronym, but a specific type of 

CCD 

Charge-Coupled Device 


glass fiber manufactured by Whatman. 

CDOM 

Chromophoric Dissolved Organic Matter 

GLI 

Global Imager 

CHN 

Carbon, Hydrogen and Nitrogen 

GoCal 

Gulf of California 



GPS 

Global Positioning System 
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GSFC 

Goddard Space Flight Center 

MERIS 

Medium Resolution Imaging Spec- 

HIVE 

High-Latitude Intercomparison and 


trometer 


Validation Experiment 

MFRSSR 

Marine Fast Rotating Shadow-band 

HOBI 

Hydro-Optics, Biology and 


Spectral Radiometer 


Instrumentation 

MI SR 

Multi-angle Imaging SpectroRadiometer 

HOT 

Hawaii Ocean Time series 

MOBY 

Marine Optical Buoy 

HPLC 

High Performance Liquid 

MOD APS 

Modular Ocean Data and Power System 


Chromatography 

MODIS 

Moderate Resolution Imaging 

HQ 

Headquarters 


Spectroradiometer 

HRPT 

High Resolution Picture Transmission 

MOS 

Modular Optoelectronic Scanner 

ICESS 

Institute for Computational Earth 

MSI 12 

Multi-Sensor level- IB to level-2 code 


Science System 

MTPE 

Mission to Planet Earth 

IDL 

Interactive Data Language 

NASA 

National Aeronautics and Space 

INDOEX 

Indian Ocean Experiment 


Administration 

IOCCG 

International Ocean Color Coordinating 

NASDA 

National Space Development Agency 


Group 


of Japan 

IOP 

Inherent Optical Properties 

NAVOCEANO Naval Oceanographic Office 

IR 

Infrared 

NDBC 

National Data Buoy Center 

IRS 

Indian Remote Sensing Satellite 

NIMBUS 

Not an acronym, but a series of NASA 

ISCCP 

International Satellite Cloud 


experimental weather satellites contain- 


Climatology Project 


ing a wide variety of atmosphere, ice, 

ISPO 

Interim SIMBIOS Project 0 Tice 


and ocean sensors. 

JGOFS 

Joint Global Ocean Flux Study 

NIST 

National Institute of Standards and 

JPL 

Jet Propulsion Laboratory 


Technology 

LAC 

Local Area Coverage, tine resolution 

NOAA 

National Oceanic and Atmospheric Ad- 


satellite data with a nominal ground re- 


ministration 


solution at nadir of approximately lKm 

NODC 

National Oceanographic Data Center 

LIDAR 

Light Detection and Ranging Instrument 

NOPP 

NIMBUS Observation Processing Sys- 

LED 

Light Emitting Dicode 


tem 

LOA 

Laboratorie d’Optique Atmospherique 

NRA 

NASA Research Announcement 

MARMAP 

Marine Resources Monitoring, Assess- 

NSF 

National Science Foundation 


ment and Prediction 

NSPO 

National Space Program Office 

MBARI 

Monterey Bay Aquarium Research In- 

NTOU 

National Taiwan Ocean University 


stitute 

OCI 

Ocean Color Imager 

MB 

Megabyte, or about one million bytes 

OCTS 

Ocean Color and Temperature Sensor 

MER 

Multispectral Environmental Radiometer 

ONR 

Office of Naval Research 



ORCA 

Optical Research Consortium 
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of the Arctic 

SeaDAS 

SeaWiFs Data Analysis System 

osc 

Orbital Sciences Corporation 

SeaWiFS 

Sea-viewing Wide Field-of-view Sensor 

osu 

Oregon State University 

SIMBAD 

The name of a sun photometer 

PAR 

Photosynthetically Available Radiation 

SIMBIOS 

Sensor Intercomparison and Merger for 

PHILLS 

Hyperspectral Imager for Low Light 


Biological and Interdisciplinary Oceanic 


Spectroscopy 


Studies 

PI 

Principal Investigator 

SIO 

Scripps Institution of Oceanography 

PM-1 

Not a acronym, used to designate the 

SGI 

Silicon Graphics, Inc. 


afternoon platform of EOS 

SIRREX 

SeaWiFS Intercalibration Round-Robin 

POLDER 

Polarization Detecting Environmental 


Experiment 


Radiometer (France) 

SLOWDROP Slow Descendent Rate Optical Profiler 

PRR 

Profiling Reflectance Radiometer 

SMSR 

SeaWiFS Multispectral Surface 

R&D 

Research and Development 


Reference 

R/V 

Research Vessel 

SOOP 

Ship of Opportunity Program 

ROCSAT 

Republic of China (Taiwan) Satellite 

SOW 

Statement of Work 


platform for the OCI sensor. 

SPMR 

SeaWiFS Profiling Multi-channel 

RR 

Round-Robin 


Radiometer 

SAB 

South Atlantic Bight 

SQM 

SeaWiFS Quality Monitor 

SBE 

Sea-Bird Electronics 

ST 

Science Team 

SCOR 

Scientific Committee on Oceanic Re- 

TB 

Terabytes 


search 

TOA 

Top of the Atmosphere 

SDDC 

Science Data Distribution Center 

TOMS 

Total Ozone Mapping Spectrometer 

SeaOPS 

SeaWiFS Optical Profiling System 

TOTO 

Tongue of the Ocean 

SeaBAM 

SeaWiFS Bio-optical Algorithms Mini- 

TSM 

Total Suspended Matter 


workshop 

TSRB-II 

Tethered Spectral Radiometer Buoy - II 

SeaBASS 

SeaWiFS Bio-optical Archive and Stor- 
age System 

XBT 

Expendable Bathythermograph 
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